Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 5.8% intergenic (+20/‑413) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6720AG 5.8% 784.9 / 30.4 285intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (140/113);  new base G (11/5);  total (151/134)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

Only 100 of 290 total aligned reads displayed.

GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCG  >  NC_000913/781523‑781821
                                                                                                                                                     |                                                                                                                                                     
cctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGaaa                                                                                                                                                       <  1:4088193‑M1/96‑1 (MQ=255)
 ctcccttactaggtgggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCATTTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                      >  1:5691761‑M1/54‑150 (MQ=255)
 ctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                      >  1:2290530‑M1/54‑150 (MQ=255)
  tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                     <  1:50802‑M1/98‑1 (MQ=255)
  cATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                     <  1:2480484/150‑1 (MQ=255)
   cccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  2:5744974‑M1/99‑1 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  1:1640690/150‑1 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    >  1:2231390/1‑150 (MQ=255)
    tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   >  1:1167558/1‑150 (MQ=255)
    tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   >  1:1972309/1‑150 (MQ=255)
    ccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  1:2628345‑M1/100‑1 (MQ=255)
     ccTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAATGGGTCGTTAGCTCAGTTGGTAGTGCATTTGACTTTTATTCTTTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                                  <  1:4262361/150‑1 (MQ=21)
     ccTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                                  >  1:2952275/1‑150 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  1:2021427‑M1/102‑1 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  2:2319824‑M1/102‑1 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  2:2369126‑M1/102‑1 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCt                                                                                                                                                 <  1:4970108‑M1/102‑1 (MQ=255)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                              <  2:5639313/150‑1 (MQ=255)
          aCGACCCACCATCCTGAATGATTAAGGCAGCATAACCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                             <  1:2016349/150‑1 (MQ=255)
           cggagggggtcgggggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCAATGGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                                            >  1:3030158‑M1/44‑150 (MQ=255)
           cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                                            <  2:3466920/150‑1 (MQ=255)
           cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                            <  1:1025150/150‑1 (MQ=255)
           aggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                            >  2:4030966‑M1/44‑150 (MQ=255)
            gACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                           <  1:3848423/150‑1 (MQ=255)
            gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                           <  1:3996379/150‑1 (MQ=255)
            gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAATTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                           <  1:4532626/150‑1 (MQ=255)
             ccccACCAGCCTGACGGATTAAGGCCGCATAAGCCCGCAAGGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCGGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCttttttac                                                                                                                                          <  1:2352883/149‑3 (MQ=17)
             aCCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGc                                                                                                                                          >  1:95062/1‑150 (MQ=255)
             aCCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGc                                                                                                                                          >  2:160929/1‑150 (MQ=255)
              cccACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCt                                                                                                                                         >  1:4475735/1‑150 (MQ=255)
                ccccaTCGTGAATGATTAAGGCAGCTTTATCCCGCATGGGGTCGGTAGCTCCGTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTan                                                                                                                                       <  1:3240684/148‑2 (MQ=17)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                       >  1:765618/1‑150 (MQ=255)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                       >  2:3649375/1‑150 (MQ=255)
                     tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                                                                                                                                  >  1:5573264/1‑150 (MQ=255)
                      ccTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                                 <  2:1889888/150‑1 (MQ=255)
                         gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATAc                                                                                                                              >  1:3648298/1‑150 (MQ=255)
                            tGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTATTTTGCTAGCGGCGATACTcc                                                                                                                           >  2:4421129/1‑149 (MQ=255)
                              ctgaaGGCCGCCTAAGCCCGCCCGGGGTCGTTAGGTGTGTTGGGAGTGGTGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                         <  1:1941209/147‑1 (MQ=11)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                         >  1:327621/1‑150 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                         <  2:1487114/150‑1 (MQ=255)
                                   gtcagcaTAAACCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                    <  1:529457/148‑1 (MQ=255)
                                     cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                                  <  1:2803512/150‑1 (MQ=255)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGTGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                                 <  2:4211236/150‑1 (MQ=255)
                                          tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCATTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                             <  1:5720230/150‑1 (MQ=255)
                                            aTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                                           >  1:342778/1‑150 (MQ=255)
                                               ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcagca                                                                                                        <  2:4475587/150‑1 (MQ=38)
                                               ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTgcagca                                                                                                        >  1:3710979/1‑150 (MQ=38)
                                                cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                                       >  2:1385963/1‑150 (MQ=37)
                                                   aaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCTGCAGGtt                                                                                                    <  1:2436200/150‑1 (MQ=21)
                                                      gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                                 <  2:2581609/150‑1 (MQ=32)
                                                      gggTCGTTAGCTCAGGTGGGCGCGCAGTTGACGGTGCGTCCCGTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                                 <  1:585510/150‑1 (MQ=11)
                                                        gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                               <  1:1546422/150‑1 (MQ=17)
                                                           gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                            <  1:4443895/150‑1 (MQ=32)
                                                           gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTttgagttga                                                                                            >  2:1084395/1‑150 (MQ=17)
                                                            ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                           <  2:1740224/150‑1 (MQ=32)
                                                            ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                           <  1:1251278/150‑1 (MQ=32)
                                                               gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCGGGTTTGAGTTGAgcgc                                                                                        >  2:4166023/1‑150 (MQ=17)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCGAAAGGCGCTTTTTTGCTATCTGCGATACTCCAGGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                      >  1:5448254/1‑150 (MQ=255)
                                                                   aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGGGCGATACTTAAAGTTTCGTAGGTGCAGCTGGTTTGATTTTAgctgagct                                                                                    >  2:1455731/1‑144 (MQ=11)
                                                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                                  >  1:175745/1‑150 (MQ=32)
                                                                        gTAGAGCAGTTGACTTTTAATCAAGTGGTCGCCGGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTATTGCTATCTGCGCTACTCCAAGCTTCGAACCAGCCGCAGGGTTGAGGTGAGCGCAGCGaaccc                                                                               >  2:1416701/1‑147 (MQ=11)
                                                                             gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                          <  2:4291222/150‑1 (MQ=32)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                       <  2:5976834/150‑1 (MQ=32)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                    <  2:765618/150‑1 (MQ=17)
                                                                                     cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                                  <  2:3326332/150‑1 (MQ=32)
                                                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                 >  1:3156862/1‑150 (MQ=32)
                                                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                 >  1:1377411/1‑150 (MQ=32)
                                                                                           aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                            <  2:496559/150‑1 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           <  2:1972309/150‑1 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           >  1:785559/1‑150 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGATTAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           >  2:2776911/1‑150 (MQ=12)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                         >  2:4759563/1‑150 (MQ=32)
                                                                                                    tCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                   <  2:2945513/150‑1 (MQ=17)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                                  <  2:947530/80‑1 (MQ=34)
                                                                                                                    tGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                   <  2:5662951/150‑1 (MQ=17)
                                                                                                                    tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                   >  2:5152345/1‑150 (MQ=32)
                                                                                                                         gACCCACCTATGTAATAACGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAACGATTCGAACCTGCAGCAGGTTTGAGTTGTGCGCAGCGATACACCGGAGCCGGTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                              <  2:5404390/150‑1 (MQ=11)
                                                                                                                         gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                              >  1:4019858/1‑150 (MQ=32)
                                                                                                                           cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                            >  1:5231432/1‑150 (MQ=17)
                                                                                                                             caccaATGAAAAAAAGCCCCCTACCGGCGTTTTTTTGCTTTCTGCTATACACAAAGATCCGAACCTGAAGCTGGTTTGAGTGGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                          <  1:5353458/150‑1 (MQ=11)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                         >  1:3992594/1‑150 (MQ=32)
                                                                                                                                 aaggTAAAAAAGCGCCCTCCCGGCGATTTTTTGCTACATGCGATACAAAAATATTCGACCCTCCGGCTGTTTTGAGTTGAGCCCACCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                      <  2:12116/147‑1 (MQ=11)
                                                                                                                                  aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAGGGGCGAGCGAAGCGAGTCATCCTGCACGAc                     >  2:2198181/1‑150 (MQ=17)
                                                                                                                                      aacaaaGCGCCCTACCGGCGCTTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCCGTTTTGTGTTGCGCGCACCTAAACACCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                 <  1:1410795/147‑1 (MQ=11)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                 >  2:3646403/1‑150 (MQ=12)
                                                                                                                                      aaaaaaGCGCCTTTAGGCCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                 >  2:4616549/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGCGCCGTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCCAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccccc                 >  1:4664626/1‑148 (MQ=255)
                                                                                                                                       aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                >  2:3066165/1‑150 (MQ=32)
                                                                                                                                          aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt             >  1:631173/1‑150 (MQ=32)
                                                                                                                                            gcgcCCTAATGGCTCTTTTTTGCTATCTGCGATAATCAAAGATTCTAACCTGCTGCTGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGt           <  1:4365304/150‑1 (MQ=11)
                                                                                                                                              gcCCTAATGGCTCTTTTTTGCTATCTGCGATACTCAAAGATTCTAACCTGAAGCCGGTTTGTGTTGAGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa         <  2:2172719/150‑1 (MQ=11)
                                                                                                                                              gcCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa         <  1:5969793/150‑1 (MQ=17)
                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa      >  1:700751/1‑150 (MQ=32)
                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa      >  2:3781305/1‑150 (MQ=32)
                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGAAACTCAAAGATTCGAACCTTCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa      >  1:3534970/1‑150 (MQ=17)
                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     >  2:3839952/1‑150 (MQ=32)
                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     >  2:5498730/1‑150 (MQ=32)
                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     >  1:2091229/1‑150 (MQ=32)
                                                                                                                                                    aaGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgc   >  2:4940457/1‑150 (MQ=32)
                                                                                                                                                    aaGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgc   >  2:5065622/1‑150 (MQ=17)
                                                                                                                                                     |                                                                                                                                                     
GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCG  >  NC_000913/781523‑781821

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