Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,670 A→T 5.6% intergenic (+18/‑415) 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,6700AT 5.6% 794.2 / 31.5 285intergenic (+18/‑415)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (139/125);  new base T (11/5);  total (151/134)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 297 total aligned reads displayed.

GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781521‑781818
                                                                                                                                                     |                                                                                                                                                    
gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                                      >  2:2813298/1‑150 (MQ=255)
gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                                      >  1:971864/1‑150 (MQ=255)
gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTAGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                                      >  2:638859/1‑150 (MQ=255)
cacctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                                      <  2:587070‑M1/94‑1 (MQ=255)
  gTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                    <  1:4517378/150‑1 (MQ=255)
   ctcccttactaggtgggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCATTTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                   >  1:5691761‑M1/54‑150 (MQ=255)
   ctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                   <  1:4007818‑M1/97‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                  <  1:50802‑M1/98‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCATGCACGACCCACCAATGTAAAAAAGCGCCCAAAAgg                                                                                                                                                  <  1:4395891‑M1/98‑1 (MQ=255)
     aTCCTGCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                 <  1:1640690/150‑1 (MQ=255)
      tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                >  1:1167558/1‑150 (MQ=255)
      tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                >  1:1972309/1‑150 (MQ=255)
       ccTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                               <  2:4850467/150‑1 (MQ=255)
       ccTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                               >  1:2952275/1‑150 (MQ=255)
        ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAggcgcg                                                                                                                                              <  1:710498‑M1/102‑2 (MQ=255)
        ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                              <  2:2319824‑M1/102‑1 (MQ=255)
        ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                              <  2:2369126‑M1/102‑1 (MQ=255)
           cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTAAAGGCGCttgt                                                                                                                                           >  1:3263003/1‑148 (MQ=255)
            aaggtgggggtcggcggttcgatcccgtcatcacccaccacttgGGTTCTTTAGCTCAGTTGGTAGCGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                          >  1:4930058‑M1/45‑150 (MQ=255)
            aaggagggggtcggcgtttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                          >  1:2556188‑M1/45‑150 (MQ=255)
            aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAACAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                          >  1:4331377/1‑150 (MQ=255)
            aCGACCCACCATCCTGAATGATTAAGGCAGCATAACCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                          <  1:2016349/150‑1 (MQ=255)
             cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                                         <  2:3466920/150‑1 (MQ=255)
             cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                         <  1:1025150/150‑1 (MQ=255)
             aggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                         >  1:3669317‑M1/44‑150 (MQ=255)
             aggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                         >  2:4030966‑M1/44‑150 (MQ=255)
              gccccACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                        <  1:3996340/148‑1 (MQ=255)
                  caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                    >  2:3649375/1‑150 (MQ=255)
                     caTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTg                                                                                                                                 <  1:4147805/150‑1 (MQ=255)
                       tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                                                                                                                               >  1:5573264/1‑150 (MQ=255)
                        ccTGAATCAATAAGGCAGCATAAGCCCGGCAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                              <  2:5733073/150‑1 (MQ=14)
                                ctgaaGGCCGCCTAAGCCCGCCCGGGGTCGTTAGGTGTGTTGGGAGTGGTGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                      <  1:1941209/147‑1 (MQ=11)
                                aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                      <  1:3860426/150‑1 (MQ=255)
                                aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                      >  1:2143114/1‑150 (MQ=255)
                                     gtcttcaTATTCCCGCAAGGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTATTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                 <  1:934333/145‑1 (MQ=255)
                                     gtcagcaTAAACCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                 <  1:529457/148‑1 (MQ=255)
                                      gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                                <  2:1704241/150‑1 (MQ=255)
                                       cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTGTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGTTATCTGCGATACTAAAAGATTCGa                                                                                                               >  2:5404608/1‑150 (MQ=255)
                                        agccTAATCCCCCAAGGGGTCGTTAGCTCAGTCGGTAGAGCAGTCGACTCTGCATCCAGTGGTCGCAGGTACGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                              <  2:4210753/150‑1 (MQ=14)
                                           ctAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCt                                                                                                           <  2:5381165/149‑1 (MQ=255)
                                            tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCATTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                          <  1:5720230/150‑1 (MQ=255)
                                                 ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcagca                                                                                                     <  2:4475587/150‑1 (MQ=38)
                                                 ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATACAGCACCACCCACCAATGTAAAAAATCGCCCTAAAGGCGCTTTTTTTCTATCTGCGATAAAAAAAGATTAGAACCAgcagca                                                                                                     >  1:232585/1‑150 (MQ=11)
                                                  cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                                    >  2:1385963/1‑150 (MQ=37)
                                                  cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACACAAAGATTCGAACcagcagcag                                                                                                    >  1:3747949/1‑150 (MQ=255)
                                                   gCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                   <  1:2737120/150‑1 (MQ=18)
                                                        gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                              <  2:2581609/150‑1 (MQ=32)
                                                        gggTCGTTAGCTCAGGTGGGCGCGCAGTTGACGGTGCGTCCCGTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                              <  1:585510/150‑1 (MQ=11)
                                                         ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                             <  1:3059080/150‑1 (MQ=32)
                                                          gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                            <  1:1546422/150‑1 (MQ=17)
                                                           tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                           <  2:4719465/150‑1 (MQ=32)
                                                             gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGTGCGATTCTCAAAGATTCGATCCTGCAGCTGGTttgagttga                                                                                         >  2:1910529/1‑150 (MQ=255)
                                                             gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTttgagttga                                                                                         >  2:1084395/1‑150 (MQ=17)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                        <  1:698645/150‑1 (MQ=32)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                        <  1:1251278/150‑1 (MQ=32)
                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCGAAAGGCGCTTTTTTGCTATCTGCGATACTCCAGGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                   >  1:5448254/1‑150 (MQ=255)
                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                   >  2:1590276/1‑150 (MQ=32)
                                                                     aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGGGCGATACTTAAAGTTTCGTAGGTGCAGCTGGTTTGATTTTAgctgagct                                                                                 >  2:1455731/1‑144 (MQ=11)
                                                                      gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                                >  1:3436714/1‑150 (MQ=32)
                                                                       ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                               >  1:175745/1‑150 (MQ=32)
                                                                        tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                              <  1:2396850/150‑1 (MQ=32)
                                                                          gTAGAGCAGTTGACTTTTAATCAAGTGGTCGCCGGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTATTGCTATCTGCGCTACTCCAAGCTTCGAACCAGCCGCAGGGTTGAGGTGAGCGCAGCGaaccc                                                                            >  2:1416701/1‑147 (MQ=11)
                                                                             gagCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTTCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                                                                         >  1:3535398/1‑150 (MQ=17)
                                                                               gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                       <  2:4430709/150‑1 (MQ=32)
                                                                                     gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                 <  2:765618/150‑1 (MQ=17)
                                                                                       cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                               <  2:2185180/150‑1 (MQ=17)
                                                                                       cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                               <  2:3326332/150‑1 (MQ=32)
                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  2:177739/1‑150 (MQ=32)
                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  1:1377411/1‑150 (MQ=32)
                                                                                             aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                         <  2:496559/150‑1 (MQ=32)
                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                        >  1:767480/1‑150 (MQ=32)
                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                     <  2:5573264/150‑1 (MQ=32)
                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                     >  1:2851194/1‑150 (MQ=17)
                                                                                                    ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                  >  1:2447724/1‑150 (MQ=32)
                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                <  2:2945513/150‑1 (MQ=17)
                                                                                                             ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                               <  2:947530/80‑1 (MQ=34)
                                                                                                                gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGgcgagcga                                      >  1:2252365/1‑150 (MQ=17)
                                                                                                                     cTGCACGCCCCACCCAGGTACAACAGCGCCCTAACGGCGCTTTTTGGCTATCTGCGATACTCAAAGATTCGAACCCGCAGCCGGTTTGAGTTGAGCGCAGCGAACCACCGGCGCCGCTCGCGGCGACGGCCCGACGGGCGagcgaagcga                                 <  2:4276750/150‑1 (MQ=11)
                                                                                                                      ggCACGCCCCACCAAGGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTTTCTGCGATACTCAAAGATTCGACCCTGCAGCCGGTTTGAGTTGCGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                <  2:5662886/149‑1 (MQ=11)
                                                                                                                             cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                         >  1:5231432/1‑150 (MQ=17)
                                                                                                                             cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                         >  1:4078199/1‑150 (MQ=32)
                                                                                                                             cccACCAATGTAAAAAAGCGCCCTAAAGCCGTTTTTTTGCTATCTGCTATACTCAACGATTCCAACCGGCAGCTGGTTTGCGTTGCGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                         <  1:3387148/150‑1 (MQ=11)
                                                                                                                              ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                        >  1:4568853/1‑150 (MQ=32)
                                                                                                                               caccaATGAAAAAAAGCCCCCTACCGGCGTTTTTTTGCTTTCTGCTATACACAAAGATCCGAACCTGAAGCTGGTTTGAGTGGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                       <  1:5353458/150‑1 (MQ=11)
                                                                                                                                cccccggTAAAACAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACCGGAGGCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                      <  2:5147670/144‑1 (MQ=11)
                                                                                                                                accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                      >  1:3159413/1‑150 (MQ=32)
                                                                                                                                   aaggTAAAAAAGCGCCCTCCCGGCGATTTTTTGCTACATGCGATACAAAAATATTCGACCCTCCGGCTGTTTTGAGTTGAGCCCACCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                   <  2:12116/147‑1 (MQ=11)
                                                                                                                                    aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAGGGGCGAGCGAAGCGAGTCATCCTGCACGAc                  >  2:2198181/1‑150 (MQ=17)
                                                                                                                                     tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                 >  2:2037397/1‑150 (MQ=32)
                                                                                                                                        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)
                                                                                                                                         aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc             >  2:4347228/1‑150 (MQ=12)
                                                                                                                                         aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc             >  2:3066165/1‑150 (MQ=32)
                                                                                                                                          aaaaGCGCCTTTAGGGCGCTTTTTTGATATATTCGATAATCAAAGATTCTAACCTGCAGCAGGTTTTAGTTTAGCGCAGCGAAACAACGGAGCCGCACGCGGCGACGGCCCGAAGGGCGAGAGAATCGAGTCCACCATCAAGACccacca            >  1:2047793/1‑150 (MQ=11)
                                                                                                                                                gcCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa      <  1:5969793/150‑1 (MQ=17)
                                                                                                                                                  ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa    >  1:4466591/1‑150 (MQ=32)
                                                                                                                                                   cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa   >  2:3781305/1‑150 (MQ=32)
                                                                                                                                                    tAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa  >  2:3839952/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                    
GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781521‑781818

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