Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,277,786 G→A 100% D823N (GAC→AAC)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,277,7860GA98.7% 510.3 / ‑2.0 157D823N (GAC→AAC) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (101/54);  total (103/54)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-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 157 total aligned reads displayed.

CCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAA  >  CP009273/1277641‑1277931
                                                                                                                                                 |                                                                                                                                                  
ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAAc                                                                                                                                                 <  1:732334/148‑1 (MQ=255)
  gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTg                                                                                                                                               >  1:1337467/1‑148 (MQ=255)
       tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGaca                                                                                                                                          >  2:102384/1‑146 (MQ=255)
       tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGaaa                                                                                                                                          >  2:102454/1‑148 (MQ=255)
           gCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTa                                                                                                                                      >  2:1265000/1‑148 (MQ=255)
            cACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTAc                                                                                                                                     >  2:1131510/1‑148 (MQ=255)
              ccTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACaa                                                                                                                                   <  2:1084523/148‑1 (MQ=255)
              ccTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACaa                                                                                                                                   >  2:2466839/1‑148 (MQ=255)
                 gTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGc                                                                                                                                >  1:2362709/1‑148 (MQ=255)
                 gTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGc                                                                                                                                >  1:1821207/1‑148 (MQ=255)
                 gTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGc                                                                                                                                >  2:1983112/1‑148 (MQ=255)
                 gTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGc                                                                                                                                >  2:51319/1‑148 (MQ=255)
                  tctcTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGcc                                                                                                                               <  1:873064/148‑1 (MQ=255)
                   ctctATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCa                                                                                                                              >  2:856111/1‑148 (MQ=255)
                       aTTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATcgc                                                                                                                          <  1:3029118/148‑1 (MQ=255)
                        ttCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATcgcg                                                                                                                         >  2:1308347/1‑148 (MQ=255)
                         tCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGACATCTACAAAGCCATCGCga                                                                                                                        >  2:1563999/1‑148 (MQ=255)
                                tcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAAtt                                                                                                                 >  2:3165355/1‑148 (MQ=255)
                                 catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtc                                                                                                                >  2:3203242/1‑148 (MQ=255)
                                 catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtc                                                                                                                >  2:1052246/1‑148 (MQ=255)
                                  atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct                                                                                                               <  2:2518349/148‑1 (MQ=255)
                                    tttGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTcc                                                                                                             >  1:2124994/1‑148 (MQ=255)
                                      tGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGa                                                                                                           <  1:2788568/148‑1 (MQ=255)
                                        ccGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg                                                                                                         <  1:2711954/148‑1 (MQ=255)
                                         cGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgt                                                                                                        <  1:1492559/148‑1 (MQ=255)
                                           actgCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGACCCGAAAGCAACCTGGGAAAGATAAAAACCCAACCCCACCAAATAATCCGAAgtgt                                                                                                      >  2:2843144/4‑148 (MQ=17)
                                             ggCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                                                                                                    <  1:3621329/148‑1 (MQ=255)
                                                    tGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGcc                                                                                             >  2:2332769/1‑148 (MQ=255)
                                                       tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                                                                                          >  1:2561742/1‑148 (MQ=255)
                                                       tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                                                                                          >  2:1717412/1‑148 (MQ=255)
                                                        aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt                                                                                         >  1:382381/1‑148 (MQ=255)
                                                         cGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg                                                                                        >  1:1000542/1‑148 (MQ=255)
                                                             aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                                                                                    <  1:3144949/148‑1 (MQ=255)
                                                             aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                                                                                    >  1:3251627/1‑148 (MQ=255)
                                                             aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                                                                                    >  1:3434070/1‑148 (MQ=255)
                                                                gacgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAg                                                                                 >  2:2871078/1‑148 (MQ=255)
                                                                   gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGTCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                                                                              <  2:545690/148‑1 (MQ=255)
                                                                   gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                                                                              <  1:2284205/148‑1 (MQ=255)
                                                                    acATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAAc                                                                             <  2:3330878/148‑1 (MQ=255)
                                                                     cATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAAcc                                                                            >  2:408078/1‑148 (MQ=255)
                                                                      aTGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTTCGTCGGCCATCTGGGTAAAGAAACCg                                                                           <  2:3494277/148‑1 (MQ=255)
                                                                      aTGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTTCGTCGGCCATCTGGGTAAAGAAACCg                                                                           <  2:3494296/148‑1 (MQ=255)
                                                                        gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                                                                         <  2:2179110/148‑1 (MQ=255)
                                                                        gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                                                                         >  2:2457407/1‑148 (MQ=255)
                                                                        gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                                                                         >  1:2365740/1‑148 (MQ=255)
                                                                        gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                                                                         >  1:79336/1‑148 (MQ=255)
                                                                          aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACAt                                                                       >  2:1385666/1‑148 (MQ=255)
                                                                          aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACAt                                                                       <  2:130216/148‑1 (MQ=255)
                                                                          aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACAGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACAt                                                                       >  1:2226578/1‑148 (MQ=255)
                                                                           tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATc                                                                      >  1:1213469/1‑148 (MQ=255)
                                                                                 ggATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACg                                                                <  2:2324224/148‑1 (MQ=255)
                                                                                  gATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGc                                                               <  2:2439117/148‑1 (MQ=255)
                                                                                   atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCTGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCt                                                              >  1:31709/1‑148 (MQ=255)
                                                                                    tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTg                                                             <  1:1390818/148‑1 (MQ=255)
                                                                                       gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCt                                                          <  1:2487413/148‑1 (MQ=255)
                                                                                       gCATCCGTTTATTAACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAAATGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGAGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCt                                                          >  2:3600942/1‑148 (MQ=255)
                                                                                         aTCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTAt                                                        <  1:1901451/148‑1 (MQ=255)
                                                                                          tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAGGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATc                                                       >  1:566475/1‑148 (MQ=255)
                                                                                          tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATc                                                       >  2:3645122/1‑148 (MQ=255)
                                                                                           ccGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATcc                                                      >  2:1190754/1‑148 (MQ=255)
                                                                                           ccGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATcc                                                      >  2:2458757/1‑148 (MQ=255)
                                                                                              tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGc                                                   >  2:2922574/1‑148 (MQ=255)
                                                                                              tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGc                                                   >  1:1355730/1‑148 (MQ=255)
                                                                                               ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCa                                                  >  1:487646/1‑148 (MQ=255)
                                                                                               ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCa                                                  <  1:1983112/148‑1 (MQ=255)
                                                                                               ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCa                                                  <  1:51319/148‑1 (MQ=255)
                                                                                               ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCa                                                  >  2:2576645/1‑148 (MQ=255)
                                                                                                 aTTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATg                                                >  1:2899258/1‑148 (MQ=255)
                                                                                                  ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGa                                               <  2:2692306/148‑1 (MQ=255)
                                                                                                      cccGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGActct                                           >  1:791286/1‑148 (MQ=255)
                                                                                                       ccGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTg                                          >  2:1679881/1‑148 (MQ=255)
                                                                                                       ccGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTg                                          >  1:2538311/1‑148 (MQ=255)
                                                                                                       ccGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTg                                          >  1:2538776/1‑148 (MQ=255)
                                                                                                         gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGcc                                        >  1:3551054/1‑148 (MQ=255)
                                                                                                         gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGcc                                        >  1:3551070/1‑148 (MQ=255)
                                                                                                          cTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCg                                       <  1:1731463/148‑1 (MQ=255)
                                                                                                                 cggcggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACt                                <  1:3529604/148‑1 (MQ=255)
                                                                                                                 cggcggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACt                                <  1:2051117/148‑1 (MQ=255)
                                                                                                                   gcggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTgg                              <  1:1029871/148‑1 (MQ=255)
                                                                                                                   gcggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTgg                              <  1:1131510/148‑1 (MQ=255)
                                                                                                                    cggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGgc                             >  2:3451287/1‑148 (MQ=255)
                                                                                                                    cggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGgc                             >  2:3451304/1‑148 (MQ=255)
                                                                                                                           tCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTCAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCg                      >  2:1950570/1‑148 (MQ=255)
                                                                                                                           tCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCg                      >  2:3616575/1‑148 (MQ=255)
                                                                                                                            ccGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGc                     >  1:2604777/1‑148 (MQ=255)
                                                                                                                                ccTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGa                 <  1:269719/148‑1 (MQ=255)
                                                                                                                                  tGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtg               <  2:1821207/148‑1 (MQ=255)
                                                                                                                                    ggAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtgtg             >  1:2162353/1‑148 (MQ=255)
                                                                                                                                    ggAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtgtg             >  1:2162861/1‑148 (MQ=255)
                                                                                                                                        gCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAg         <  1:1296532/148‑1 (MQ=255)
                                                                                                                                        gCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAg         >  2:2325780/1‑148 (MQ=255)
                                                                                                                                        gCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCAGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAg         >  2:1904553/1‑148 (MQ=255)
                                                                                                                                         cGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGa        >  1:1131869/1‑148 (MQ=255)
                                                                                                                                           aaaaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACt      >  2:479175/1‑148 (MQ=255)
                                                                                                                                           aaaaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACt      >  2:1534810/1‑148 (MQ=255)
                                                                                                                                             aaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTgg    >  2:394415/1‑148 (MQ=255)
                                                                                                                                             aaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTgg    >  2:1944502/1‑148 (MQ=255)
                                                                                                                                             aaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTgg    >  2:567682/1‑148 (MQ=255)
                                                                                                                                              aGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGa   >  2:1796073/1‑148 (MQ=255)
                                                                                                                                               gCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGaa  >  2:2976600/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAA  >  CP009273/1277641‑1277931

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: