Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 786,632 A→G T334A (ACC→GCC)  aroG → 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,6320AG98.8% 304.7 / ‑4.5 82T334A (ACC→GCC) aroG3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/47);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGG  >  NC_000913/786495‑786768
                                                                                                                                         |                                                                                                                                        
agATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGcc                                                                                                                                        >  2:240364/1‑140 (MQ=255)
 gATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCTATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCg                                                                                                                                       <  2:652785/140‑1 (MQ=255)
 gATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCg                                                                                                                                       <  2:945555/140‑1 (MQ=255)
 gATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCg                                                                                                                                       <  2:759812/140‑1 (MQ=255)
 gATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCg                                                                                                                                       <  2:528161/140‑1 (MQ=255)
  aTTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATTACCGATGCCTGCATCGGCTGGGAAGATGCCgc                                                                                                                                      >  2:713968/1‑139 (MQ=255)
    tGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATg                                                                                                                                    <  1:535197/140‑1 (MQ=255)
    tGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATg                                                                                                                                    <  2:437710/140‑1 (MQ=255)
    tGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATg                                                                                                                                    <  2:910738/140‑1 (MQ=255)
      ccGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGTCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGct                                                                                                                                  <  1:1375281/140‑1 (MQ=255)
       cGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGctc                                                                                                                                 <  1:1211373/140‑1 (MQ=255)
                aaaGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTc                                                                                                                        >  1:314428/1‑140 (MQ=255)
                  aGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCaa                                                                                                                      <  2:238139/140‑1 (MQ=255)
                    gCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACt                                                                                                                    >  1:48141/1‑140 (MQ=255)
                      cATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACAGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTgg                                                                                                                  >  2:207991/1‑140 (MQ=255)
                       attattGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGc                                                                                                                 <  1:1348768/140‑1 (MQ=255)
                        ttattGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCg                                                                                                                <  2:1417386/140‑1 (MQ=255)
                          attGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGaa                                                                                                              <  2:383667/140‑1 (MQ=255)
                            tGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATg                                                                                                            <  2:929338/140‑1 (MQ=255)
                                   aTGGTGGAAAGCCATCTGGTGGAAGGCAATCGGAGCCTCGAGAGCGGGGGGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGGATGCAGTaaa                                                                                                     >  2:840690/1‑140 (MQ=255)
                                   aTGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGaaaa                                                                                                     >  2:1310195/1‑140 (MQ=255)
                                   aTGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTaaa                                                                                                     >  2:731579/1‑140 (MQ=255)
                                   aTGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTaaa                                                                                                     >  1:524429/1‑140 (MQ=255)
                                           aaGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTc                                                                                             <  1:417350/140‑1 (MQ=255)
                                             gCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTcgc                                                                                           >  1:1472310/1‑140 (MQ=255)
                                             gCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTcgc                                                                                           <  2:1486657/140‑1 (MQ=255)
                                              ccATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTcgcg                                                                                          >  2:254229/1‑140 (MQ=255)
                                              ccATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAAGGCAGTAAAAGCGCGTcgcg                                                                                          >  1:1139472/1‑140 (MQ=255)
                                                aTCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCggg                                                                                        <  2:979243/140‑1 (MQ=255)
                                                aTCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCggg                                                                                        >  1:1163545/1‑140 (MQ=255)
                                                    ggtggAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAg                                                                                    <  2:653093/140‑1 (MQ=255)
                                                      tggAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGt                                                                                  <  2:748416/140‑1 (MQ=255)
                                                             cAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTg                                                                           >  1:634724/1‑140 (MQ=255)
                                                             cAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTg                                                                           >  1:346130/1‑140 (MQ=255)
                                                             cAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGCTTAATTg                                                                           >  1:705438/1‑140 (MQ=255)
                                                               aTCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTc                                                                         >  1:1067157/1‑140 (MQ=255)
                                                                tCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCg                                                                        >  2:149676/1‑140 (MQ=255)
                                                                  agagCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGa                                                                      <  2:1412896/140‑1 (MQ=255)
                                                                   gagCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGAt                                                                     <  2:314428/140‑1 (MQ=255)
                                                                    agCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATg                                                                    <  2:524429/140‑1 (MQ=255)
                                                                     gCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATgc                                                                   <  1:207991/140‑1 (MQ=255)
                                                                      ccTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGa                                                                      >  2:1281324/1‑136 (MQ=255)
                                                                      ccTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGa                                                                      <  1:1281324/136‑1 (MQ=255)
                                                                        tCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGcc                                                                <  1:240364/140‑1 (MQ=255)
                                                                         cGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCg                                                               <  2:1021347/140‑1 (MQ=255)
                                                                         cGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCg                                                               <  2:1056094/140‑1 (MQ=255)
                                                                          gagaGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGt                                                              <  2:1082102/140‑1 (MQ=255)
                                                                          gagaGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGt                                                              <  2:1451025/140‑1 (MQ=255)
                                                                          gagaGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGt                                                              <  2:1435347/140‑1 (MQ=255)
                                                                              gCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCaga                                                          <  2:1378703/140‑1 (MQ=255)
                                                                                  ggAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTgg                                                      >  1:386014/1‑140 (MQ=255)
                                                                                   gAGCCGCTGGCCTACGGTAAGAGCATTACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGc                                                     <  1:713968/140‑1 (MQ=255)
                                                                                    aGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCg                                                    <  2:848106/140‑1 (MQ=255)
                                                                                    aGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCg                                                    <  2:1246101/140‑1 (MQ=255)
                                                                                      ccGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTa                                                  <  2:1430224/140‑1 (MQ=255)
                                                                                          tGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCg                                              <  1:149676/140‑1 (MQ=255)
                                                                                           ggCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGa                                             <  1:254229/140‑1 (MQ=255)
                                                                                           ggCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGa                                             <  2:1163545/140‑1 (MQ=255)
                                                                                               tACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGaa                                         >  1:638205/1‑140 (MQ=255)
                                                                                                aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAAt                                        <  2:300673/140‑1 (MQ=255)
                                                                                                aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAAt                                        <  2:1342780/140‑1 (MQ=255)
                                                                                                aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAAt                                        <  1:683816/140‑1 (MQ=255)
                                                                                                aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAAt                                        <  1:520097/140‑1 (MQ=255)
                                                                                                    tAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcacc                                    >  2:883549/1‑140 (MQ=255)
                                                                                                     aaGAGCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcacca                                   <  2:820601/140‑1 (MQ=255)
                                                                                                      agagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccac                                  <  2:498104/140‑1 (MQ=255)
                                                                                                      agagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccac                                  >  2:259797/1‑140 (MQ=255)
                                                                                                      agagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccac                                  >  2:235610/1‑140 (MQ=255)
                                                                                                      agagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccac                                  >  1:1242725/1‑140 (MQ=255)
                                                                                                      agagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGGTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccac                                  >  2:931491/1‑140 (MQ=255)
                                                                                                       gagCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACcaca                                 >  1:1123182/1‑140 (MQ=255)
                                                                                                        agCATCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAg                                <  2:48141/140‑1 (MQ=255)
                                                                                                           aTCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGcc                             >  2:1510890/1‑140 (MQ=255)
                                                                                                           aTCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGcc                             >  2:748844/1‑140 (MQ=255)
                                                                                                            tCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCt                            >  1:102016/1‑140 (MQ=255)
                                                                                                            tCACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCCCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCt                            >  1:1201952/1‑140 (MQ=255)
                                                                                                             cACCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTg                           >  2:164071/1‑140 (MQ=255)
                                                                                                              aCCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGa                          >  1:1043497/1‑140 (MQ=255)
                                                                                                              aCCGATGCCTGCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGa                          >  1:395318/1‑140 (MQ=255)
                                                                                                                        gCATCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAAAGCAGTAAAGGCGCGTCGCGGGAAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCgcgc                >  2:1431532/1‑140 (MQ=255)
                                                                                                                          aTCGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAg              >  2:679657/1‑140 (MQ=255)
                                                                                                                            cGGCTGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCg            <  2:1043497/140‑1 (MQ=255)
                                                                                                                                tGGGAAGATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGc        <  1:1310195/140‑1 (MQ=255)
                                                                                                                                      gATGCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAgg  <  2:346130/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
AGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGG  >  NC_000913/786495‑786768

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: