Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 81,102 G→A S286L (TCG→TTG)  leuB ← 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091381,1020GA95.6% 119.6 / ‑3.8 46S286L (TCG→TTG) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (28/15);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGATCTCCCG  >  NC_000913/80856‑81343
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cTTAGCCATGATTACACCCCTTCTGCTACATGGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGAGTTCAAAGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCCGCAGTGACAGCGAAAGGAATTGTGCAATCGGGTTTGCGGTGGTTTTGCCT‑GGGGTATCTGGGGCCAAgcc                                                                                                                                                                                                                                               >  2:287214/1‑251 (MQ=255)
             acacCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCgggtct                                                                                                                                                                                                                                  >  2:481201/1‑247 (MQ=255)
               acCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGGTGGCGGTGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCAt                                                                                                                                                                                                                                 >  2:178279/1‑250 (MQ=255)
                               agagggCAATGATACCGCCCATTCCATCGGTACTACCGGCGGCAGCGCCGCGGGCTAAATCCCCGGTGCGAATGCCTTCTTCAAAGGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGc                                                                                                                                                                                                                <  1:208948/248‑1 (MQ=255)
                               aGCGGGAAAGGATACCCCCCATTTCATCGGTCCTAACGGCGGCAGCGCCCGGGGCTAAATCCCCGGTGCGAATGCCTTCTTCAAATGCGCGGTTAAGGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTGGCCT‑GCGATATCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGc                                                                                                                                                                                                                <  1:481201/251‑1 (MQ=255)
                                   ggCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGt                                                                                                                                                                                                            >  1:311679/1‑251 (MQ=255)
                                    gCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGtt                                                                                                                                                                                                           >  1:128187/1‑251 (MQ=255)
                                              cGCCCATTTCATCGGTCCAACCGGCGCAAGCGCCACGGGCTAAACCCCCGGGCCGAAGGCTTTCTTCAAATCCGCGTTTATTGGCGCGTTCAATGGCGCAAGCCGCATCAGCGGCATCAAGGCGGTAACGCAGCAGCATTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGAGGTTGTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCg                                                                                                                                                                                                 <  1:28822/251‑1 (MQ=255)
                                                 ccATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGa                                                                                                                                                                                               >  1:459708/1‑250 (MQ=255)
                                                 ccATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGa                                                                                                                                                                                               >  1:166484/1‑250 (MQ=255)
                                                          cGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTGATGCGCGGGTCATGGCGCGTTCAATGGCGCAAGCCGCCTCATCGGCATCCAAGCTGTAACGCCGCAGCCGGGCCAGCGAAAGGATTTGTGCAATCGGGGTGGCGAAGTTTTTTCCT‑GCCAACTTTGGGACCAAGCAGCCCGCCCGTTCATCCAGTACAAAACCTTGAGCGTTCAGGGTGGGGGGAAGCAACacc                                                                                                                                                                                     >  2:358782/1‑249 (MQ=255)
                                                          cGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCgg                                                                                                                                                                                                                                      <  1:87060/202‑1 (MQ=255)
                                                          cGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCgg                                                                                                                                                                                                                                      >  2:87060/1‑202 (MQ=255)
                                                                 aaCGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGCATGCCTTCTTCTAATGCGCGGGTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGGCAAACCGCCCGCCGGGTCATACAGTCCAAAACCTTGCTCGTTCTGGCTGGCGGAAGGCAACATCCCCCTCg                                                                                                                                                                              >  2:330183/1‑251 (MQ=255)
                                                                              ccACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGCTCAt                                                                                                                                                                 >  1:181061/1‑251 (MQ=255)
                                                                              ccACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCAt                                                                                                                                                                 >  1:288276/1‑251 (MQ=255)
                                                                                  ggcctaatcccccGGTGCGAATTCCTTCTTCAACTGCGCGTTTAATGGCTCTTTCAATGGCGCAAGCCGCATCAGCGGCATCCAGGCTATAACGCAGCAGCAGTCCCAGCAAAAGGATTTGTTCAATCGGGTTGGCGATTTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTgcg                                                                                                                                                             <  1:251775/242‑1 (MQ=255)
                                                                                         acccccGGTGCGAATGCCTTCTTCAAATGCGCGTTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGt                                                                                                                                                      <  1:178279/249‑1 (MQ=255)
                                                                                              cGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcagac                                                                                                                                                 >  2:170158/1‑251 (MQ=255)
                                                                                                    gAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGATCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATg                                                                                                                                           >  1:118121/1‑251 (MQ=255)
                                                                                                         ccTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGAACTCGTCAGACAGAATGTCGcc                                                                                                                                      >  1:224899/1‑251 (MQ=255)
                                                                                                            tattcaaaTGGGCGGTTAATGGCGCGTAAAGGGGCGAAAGCGGCACCATCGGCATCCAGGCTGTACCGCAGCAGCAGGGCCGGCGAAAGGATTTGTGCAGTCGGGTTGGCGATGTTTTTGCCT‑GCGATAGCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCa                                                                                                                                     <  1:394071/243‑1 (MQ=255)
                                                                                                                   aTGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATAAAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACCTCCCCATCGAGCCAGTGATCATTGCGCACTCTTCAACCAGAATGTCGCCAAACAggatg                                                                                                                            >  2:46959/1‑248 (MQ=255)
                                                                                                                        cgGTTAAAGGCGCGTTCAATGGCGCAAGCCGAATCATCAGCATCCAGGCTGTACCGCAGCAGCAGCGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTCCCG‑GCGATATCTGTTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTACAGGCTGGCGGAAGACCACATCACCCTCTAGCCCGTGTTCATTGCGCACTCGTCCGACAGAATGTAGCCaacacgattggagc                                                                                                                        >  2:247393/1‑238 (MQ=255)
                                                                                                                           ttAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGGGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGcac                                                                                                                      >  1:11379/1‑249 (MQ=255)
                                                                                                                                   gcgtaaaagggCGAAAGCCGCATGACCGGCAGCCAGGGTGTAACACAGCAGCAGTGCCAGCGAAGGGTTTTGTGCAAGCGGGTTGGCGATGTTTTGGCCT‑GCGATAGCGGGTGCCACGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGt                                                                                                            <  2:118121/242‑1 (MQ=255)
                                                                                                                                       ccaagggCGCAACCCGCATCACCGGCATCCAGGCTGGAACGCAGCAGCAGTGCCAGCGAAAGGTTTTGTGCAATCGGGTTGCCGATGTTTTTGCCT‑GCGATATCTGGTCCCAAGCCCCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCaaa                                                                                                        <  1:252103/246‑1 (MQ=255)
                                                                                                                                           tGGCGAAAGCCGCATCATCGGCATCCAGGCTGTAACGAAGCAGCAGTGCCAGAGCACGGATTTGTGCACTCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGACGGCCACATCCCCATCGAGCCAGTGATCATTGCGCCCTCGTCAGACAGAATGTCGCCAAACAGGTTGGTGCACAGCCGCACGTcccactgt                                                                                                    >  2:228460/1‑244 (MQ=255)
                                                                                                                                            ggCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGATGGCGGAAGGCAACCTCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACCGGGTGGAGCACAGCAGGACGTCAAACtg                                                                                                     >  2:446825/1‑249 (MQ=255)
                                                                                                                                                  aGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGTTTTGTGCAATCGTGTGGGCGATTTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGAt                                                                                             <  1:295267/251‑1 (MQ=255)
                                                                                                                                                      ccatcaccgtcctCCAGGCTGTAACGCACCAGCAGTGCCAGCGAAAGTATTTGTGCATCCGGGTTGGCGATGTTTTTGCCT‑GCAATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCttt                                                                                         <  1:287214/239‑1 (MQ=255)
                                                                                                                                                                         tGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGCTATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGGGGCg                                                                      >  2:82275/1‑251 (MQ=255)
                                                                                                                                                                          gTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGt                                                                     >  1:380102/1‑251 (MQ=255)
                                                                                                                                                                          gTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGt                                                                     >  2:200303/1‑251 (MQ=255)
                                                                                                                                                                               gcggcagcagGGCCAGCGAAAGGATTTTTGCTATCGGGTTGGCG‑TTTTTTTTCCTCGCGATATCTGGTGCAACGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTTGCGTTGTCg                                                                <  1:323580/248‑1 (MQ=255)
                                                                                                                                                                                    agcagTCCCAGCGAAAGGAGTTGGGCAACCGGGTTGGCGATGTTTTTGCCT‑GCGATACCTGGTGCCATGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACGCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTa                                                           <  1:59517/251‑1 (MQ=255)
                                                                                                                                                                                                 aaaGGTTTGGTGAAATCGGGTTGGCGATGTTTTGGCCT‑GCGATAGCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGt                                              <  1:10054/251‑1 (MQ=255)
                                                                                                                                                                                                    ggttttgggaatgcggggtgggtaTGTTTTTGCCT‑GGGATATCGGGTGCCAGGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCg                                           <  2:128187/228‑1 (MQ=255)
                                                                                                                                                                                                         gtgtcAATCGGGTGGGGGTTGTTTGTCCCT‑GCGTTTTCTGGGGCCCACCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATc                                      <  1:450316/247‑1 (MQ=255)
                                                                                                                                                                                                             cAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCggg                                  >  1:8555/1‑251 (MQ=255)
                                                                                                                                                                                                              aaTCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGt                                 >  1:253187/1‑251 (MQ=255)
                                                                                                                                                                                                               atggggTGGGCGATTTTTTTGCCT‑GCGTTATCTGGTGCCGAGCCCCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTa                                <  2:288276/248‑1 (MQ=255)
                                                                                                                                                                                                                            tGTTTTTGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGAtc                   >  2:63283/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttttGCCT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCCGTTCGACATCCGGGTATTCCGTGGCGATCTCg                >  1:440189/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccT‑GCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCAGGGTATTCCGTGGCGATCTCGTTacc           >  1:267782/1‑249 (MQ=255)
                                                                                                                                                                                                                                              tCTGGTCCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGCTCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGATCTCCCg  >  1:440331/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCT‑GCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGATCTCCCG  >  NC_000913/80856‑81343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑