Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,209,475 A→T 100% N387I (AAC→ATC)  PP_1948 → benzaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,209,4750AT99.5% 782.3 / ‑4.8 214N387I (AAC→ATC) PP_1948benzaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (117/96);  total (118/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 218 total aligned reads displayed.

GCCCCCTAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCAACAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCTCGC  >  NC_002947/2209328‑2209620
                                                                                                                                                   |                                                                                                                                                 
gCCCCCTAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGtcatca                                                                                                                                                 >  1:559378/1‑150 (MQ=255)
  ccccTAACGACCACGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCatcaac                                                                                                                                               <  2:4275296/150‑1 (MQ=255)
    ccTAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTc                                                                                                                                             <  2:3799247/150‑1 (MQ=255)
    ccTAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTc                                                                                                                                             >  2:243866/1‑150 (MQ=255)
      tAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGa                                                                                                                                           >  2:5150807/1‑150 (MQ=255)
      tAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGa                                                                                                                                           >  2:4818226/1‑150 (MQ=255)
         ccaCCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGc                                                                                                                                        <  2:5544765/148‑1 (MQ=255)
          gACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCg                                                                                                                                       >  2:1778605/1‑150 (MQ=255)
            ccAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGa                                                                                                                                     >  2:4408259/1‑150 (MQ=255)
                ggCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGc                                                                                                                                 >  1:60858/1‑150 (MQ=255)
                ggCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGc                                                                                                                                 >  1:374187/1‑150 (MQ=255)
                 gCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCt                                                                                                                                >  2:1490370/1‑150 (MQ=255)
                  cTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTc                                                                                                                               <  1:2482543/150‑1 (MQ=255)
                    cAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCgg                                                                                                                             >  1:1947786/1‑150 (MQ=255)
                      gTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGga                                                                                                                           >  2:532203/1‑150 (MQ=255)
                       tttGAGCAGGTTCACCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGgaa                                                                                                                          <  2:2093822/150‑1 (MQ=255)
                        ttGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCGCCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGgaag                                                                                                                         <  1:4922704/150‑1 (MQ=255)
                        ttGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGgaag                                                                                                                         <  1:2042480/150‑1 (MQ=255)
                            gCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAAt                                                                                                                     >  1:3104512/1‑150 (MQ=255)
                            gCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAAt                                                                                                                     >  1:1145994/1‑150 (MQ=255)
                               ggTTCATCGTTACAACGAGCGGGCCAGCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAAtttt                                                                                                                  <  2:3826179/150‑1 (MQ=255)
                                gTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTc                                                                                                                 <  1:1069094/150‑1 (MQ=255)
                                   cATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGa                                                                                                              <  2:3555892/150‑1 (MQ=255)
                                    atggTTTCATCGAGCTGGCGTTGGTGGAGGGGGGTGTGTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGAc                                                                                                             <  2:5516315/147‑1 (MQ=255)
                                     tCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGAcc                                                                                                            >  2:1628957/1‑150 (MQ=255)
                                     tCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGAcc                                                                                                            >  1:5453020/1‑150 (MQ=255)
                                        ttACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGt                                                                                                         >  1:1446666/1‑150 (MQ=255)
                                          aCAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTc                                                                                                       <  1:817816/150‑1 (MQ=255)
                                          aCAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTc                                                                                                       >  1:4562038/1‑150 (MQ=255)
                                            aaCGAGCTGGTCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTg                                                                                                     >  1:699207/1‑150 (MQ=255)
                                            aaCGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTg                                                                                                     >  1:1283843/1‑150 (MQ=255)
                                             aCGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCttg                                                                                                    >  2:3283448/1‑148 (MQ=255)
                                                  cTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTg                                                                                               >  2:177203/1‑150 (MQ=255)
                                                       cATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTcc                                                                                          <  2:4071050/150‑1 (MQ=255)
                                                       cATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTcc                                                                                          >  2:3985081/1‑150 (MQ=255)
                                                       cATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTcc                                                                                          >  1:5581581/1‑150 (MQ=255)
                                                       cATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTcc                                                                                          >  1:3009288/1‑150 (MQ=255)
                                                        aTCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCg                                                                                         >  1:409152/1‑150 (MQ=255)
                                                        aTCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCg                                                                                         >  1:319238/1‑150 (MQ=255)
                                                            aggaggGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCa                                                                                     >  2:2185646/1‑150 (MQ=255)
                                                             ggaggGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACATAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAAGTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCaa                                                                                    <  2:5331427/150‑1 (MQ=255)
                                                             ggaggGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGTACCCGTTCTGGTACTGATTCCGTTCaa                                                                                    <  2:2573321/150‑1 (MQ=255)
                                                             ggaggGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCaa                                                                                    <  2:5072309/150‑1 (MQ=255)
                                                             ggaggGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCaa                                                                                    <  1:3684509/150‑1 (MQ=255)
                                                                 ggggCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGAc                                                                                <  2:3158585/150‑1 (MQ=255)
                                                                 ggggCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGAc                                                                                <  2:115873/150‑1 (MQ=255)
                                                                      tCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGaaaa                                                                           >  2:3051708/1‑150 (MQ=255)
                                                                       cATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAAc                                                                          <  2:2557552/150‑1 (MQ=255)
                                                                       cATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAAc                                                                          <  2:2321531/150‑1 (MQ=255)
                                                                         ttGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGa                                                                        <  2:4996920/150‑1 (MQ=255)
                                                                         ttGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGa                                                                        >  1:4777163/1‑150 (MQ=255)
                                                                            cTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGc                                                                     <  2:237200/150‑1 (MQ=255)
                                                                             ttGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCg                                                                    >  2:2692397/1‑150 (MQ=255)
                                                                               ggTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGt                                                                  <  1:3865917/150‑1 (MQ=255)
                                                                                gTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTg                                                                 <  2:782555/150‑1 (MQ=255)
                                                                                 tCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTgg                                                                >  2:2594376/1‑150 (MQ=255)
                                                                                 tCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTgg                                                                <  2:664826/150‑1 (MQ=255)
                                                                                   ggTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGcc                                                              <  1:2309211/150‑1 (MQ=255)
                                                                                      ggggAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTc                                                                                                                               <  1:991654/82‑1 (MQ=255)
                                                                                         gAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATcgcg                                                        >  2:3586245/1‑150 (MQ=255)
                                                                                          aaGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGa                                                       <  2:695808/150‑1 (MQ=255)
                                                                                           aGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGaa                                                      <  2:5396737/150‑1 (MQ=255)
                                                                                              tCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGa                                                   >  1:2629454/1‑150 (MQ=255)
                                                                                                 aaGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTc                                                <  2:1416035/150‑1 (MQ=255)
                                                                                                  aGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCg                                               <  1:1778605/150‑1 (MQ=255)
                                                                                                  aGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCg                                               >  2:4511582/1‑150 (MQ=255)
                                                                                                  aGTGACAAAACCGGCTGGTACGTACGACCCACCATCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCg                                               <  2:987629/150‑1 (MQ=255)
                                                                                                   gTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCgg                                              >  2:1874173/1‑150 (MQ=255)
                                                                                                   gTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCgg                                              >  1:3430864/1‑150 (MQ=255)
                                                                                                    tGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGa                                             <  1:2278403/150‑1 (MQ=255)
                                                                                                       cAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTa                                          >  1:2359965/1‑150 (MQ=255)
                                                                                                        aaaaCCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTAc                                         >  1:1269210/1‑150 (MQ=255)
                                                                                                         aaaCCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTc                                                                                                                 <  2:274542/77‑1 (MQ=255)
                                                                                                         aaaCCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTc                                                                                                                 >  1:274542/1‑77 (MQ=255)
                                                                                                          aaCCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACgg                                       >  1:5533838/1‑150 (MQ=255)
                                                                                                           aCCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACggg                                      <  2:1096189/150‑1 (MQ=255)
                                                                                                                 tGGTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTcgcg                                >  1:1654855/1‑150 (MQ=255)
                                                                                                                  ggTACGTACGACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCgg                               >  1:1524176/1‑150 (MQ=255)
                                                                                                                     acgtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGTAGAAATTTTCGGACCCTTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTg                            <  1:2204486/150‑1 (MQ=255)
                                                                                                                     acgtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTg                            <  1:5554373/150‑1 (MQ=255)
                                                                                                                      cgtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTgg                           <  1:3546854/150‑1 (MQ=255)
                                                                                                                      cgtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTgg                           >  1:2951606/1‑150 (MQ=255)
                                                                                                                       gtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGc                          >  1:1315673/1‑150 (MQ=255)
                                                                                                                       gtacgACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGc                          <  1:853345/150‑1 (MQ=255)
                                                                                                                           gACCCACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTAt                      <  1:2813868/150‑1 (MQ=255)
                                                                                                                              ccACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGa                   >  1:3130234/1‑150 (MQ=255)
                                                                                                                               cACCGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGac                  >  2:134768/1‑150 (MQ=255)
                                                                                                                                  cGTCTACACCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGaccac               >  2:4554360/1‑150 (MQ=255)
                                                                                                                                       acacCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATc          <  2:726518/150‑1 (MQ=255)
                                                                                                                                       acacCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATc          >  1:1958670/1‑150 (MQ=255)
                                                                                                                                       acacCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATc          <  1:514625/150‑1 (MQ=255)
                                                                                                                                         acCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTTGCGTATGGACCACCGATCTg        <  2:374187/150‑1 (MQ=255)
                                                                                                                                         acCAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTg        <  1:464543/150‑1 (MQ=255)
                                                                                                                                          ccAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTgg       <  1:1702689/150‑1 (MQ=255)
                                                                                                                                          ccAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTgg       <  2:4942373/150‑1 (MQ=255)
                                                                                                                                          ccAACGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTgg       <  1:3283448/150‑1 (MQ=255)
                                                                                                                                            aaCGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCt     >  1:2523144/1‑150 (MQ=255)
                                                                                                                                            aaCGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCt     <  2:559378/150‑1 (MQ=255)
                                                                                                                                              cGTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCTcg   >  1:4135920/1‑150 (MQ=255)
                                                                                                                                               gTCATCAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCTcgc  <  2:2851226/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCCCCTAACGACCAAGGCTCAGTTTGAGCAGGTTCATCGTTACAACGAGCTGGCCATCCAGGAGGGGGCTCATTGCTTGGTCGGTGGGGAAGCTCCAAGTGACAAAACCGGCTGGTACGTACGACCCACCGTCTACACCAACGTCAACAACTCGATGCGGATTGCTCGGGAAGAAATTTTCGGACCCGTTCTGGTACTGATTCCGTTCAAGGACGAAAACGAGGCGGTGGCCATCGCGAATGACTCGGACTACGGGCTCGCGGCTGGCGTATGGACCACCGATCTGGCTCGC  >  NC_002947/2209328‑2209620

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