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.8% 835.3 / ‑0.2 240D823N (GAC→AAC) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (123/114);  total (123/117)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-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 240 total aligned reads displayed.

TCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGG  >  CP009273/1277640‑1277929
                                                                                                                                                  |                                                                                                                                               
tCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACt                                                                                                                                              >  2:2228602/1‑150 (MQ=255)
 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTg                                                                                                                                             >  1:4171472/1‑150 (MQ=255)
   gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTggg                                                                                                                                           <  1:2780033/150‑1 (MQ=255)
   gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTggg                                                                                                                                           <  1:1040905/150‑1 (MQ=255)
   gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGAGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTggg                                                                                                                                           <  1:2864236/150‑1 (MQ=255)
    tCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGa                                                                                                                                          <  1:542542/150‑1 (MQ=255)
    tCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGa                                                                                                                                          <  1:981041/150‑1 (MQ=255)
     cTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGaa                                                                                                                                         >  2:1885120/1‑150 (MQ=255)
      tGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGaaa                                                                                                                                        <  2:4508431/150‑1 (MQ=255)
        tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATc                                                                                                                                      >  2:2204380/1‑150 (MQ=255)
        tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATc                                                                                                                                      >  1:979113/1‑150 (MQ=255)
         cGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCt                                                                                                                                     >  2:2811725/1‑150 (MQ=255)
         cGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCAGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCt                                                                                                                                     >  1:4535197/1‑150 (MQ=255)
          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTa                                                                                                                                    >  2:209100/1‑150 (MQ=255)
          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTa                                                                                                                                    >  1:3503318/1‑150 (MQ=255)
          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTa                                                                                                                                    >  2:2672060/1‑150 (MQ=255)
             cACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACaa                                                                                                                                 <  1:408606/150‑1 (MQ=255)
              aCCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACaaa                                                                                                                                <  2:2114158/150‑1 (MQ=255)
              aCCTGTATCTTTTCCGACATCATTTTGCCGACGGCGACCTGTTACGAAAAAGACGACATGTATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACata                                                                                                                                <  2:610582/150‑3 (MQ=25)
               ccTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAg                                                                                                                               <  1:1411919/150‑1 (MQ=255)
               ccTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAg                                                                                                                               >  1:2840189/1‑150 (MQ=255)
                cTGTCTCTATTCCGACATCATATTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGAGGCGGTCGATCCGGCCTGGGAAGCGAAAAGAAACTGGGAAATCTACAAAGc                                                                                                                              >  1:2211303/1‑150 (MQ=255)
                   tctcTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCAt                                                                                                                           <  2:4232202/150‑1 (MQ=255)
                             gACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGaaa                                                                                                                 <  1:1193412/150‑1 (MQ=255)
                             gACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGaaa                                                                                                                 <  2:3567136/150‑1 (MQ=255)
                             gACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGaaa                                                                                                                 <  1:75674/150‑1 (MQ=255)
                                atcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtc                                                                                                              >  1:2291401/1‑150 (MQ=255)
                                 tcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct                                                                                                             <  2:3641098/150‑1 (MQ=255)
                                  catTATGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCAGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATtctc                                                                                                            >  1:3164468/1‑150 (MQ=255)
                                   atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTcc                                                                                                           <  2:157392/150‑1 (MQ=255)
                                   atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTcc                                                                                                           <  1:2687880/150‑1 (MQ=255)
                                   atTTTGACGACGGCGAACTTGTACGAAAAAGACGTCATTTATACTTCGGATATGCTTACGTTTATTCACCCGCTGTCTGCTGAGGTCGATCCGTACTGGGAAGCGAAAAGCAACTGGGAAATATTCAAATCATTCGCGAAGTAATTCTcc                                                                                                           <  2:560504/150‑1 (MQ=255)
                                       tGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg                                                                                                       >  2:4121889/1‑150 (MQ=255)
                                        gCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgt                                                                                                      <  1:215612/150‑1 (MQ=255)
                                          cGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt                                                                                                    <  1:2467474/150‑1 (MQ=255)
                                          cGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt                                                                                                    <  2:4641925/150‑1 (MQ=255)
                                          cGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt                                                                                                    <  2:3157149/150‑1 (MQ=255)
                                          cGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt                                                                                                    <  1:1257883/150‑1 (MQ=255)
                                            aCGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                                                                                                  >  2:485707/1‑150 (MQ=255)
                                              ggCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                                                                                                >  2:736028/1‑150 (MQ=255)
                                              ggCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                                                                                                >  2:3378227/1‑150 (MQ=255)
                                              ggCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                                                                                                >  1:2608729/1‑150 (MQ=255)
                                               gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACCAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTc                                                                                               <  2:4056338/150‑1 (MQ=255)
                                               gCGACCTGGTACGAAAAAGACGACATGAATAATTAGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGATGTGTGCGTc                                                                                               <  1:1123998/150‑1 (MQ=255)
                                                 gACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCgg                                                                                             >  1:4424358/1‑150 (MQ=255)
                                                  aCCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGc                                                                                            >  2:3380977/1‑150 (MQ=255)
                                                    cTGGTACGAACAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCa                                                                                          <  2:2406084/150‑1 (MQ=255)
                                                    cTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCa                                                                                          <  2:2406812/150‑1 (MQ=255)
                                                    cTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCa                                                                                          >  1:1735870/1‑150 (MQ=255)
                                                       gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt                                                                                       >  1:3079096/1‑150 (MQ=255)
                                                        tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg                                                                                      >  2:4404520/1‑150 (MQ=255)
                                                        tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg                                                                                      >  1:2271767/1‑150 (MQ=255)
                                                         ccGAAAAAGACGCCCGGAATCCTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg                                                                                     <  1:1618194/149‑1 (MQ=25)
                                                         aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg                                                                                     <  2:1500820/150‑1 (MQ=255)
                                                            aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                                                                                  >  2:392758/1‑150 (MQ=255)
                                                            aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                                                                                  <  2:595114/150‑1 (MQ=255)
                                                             aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaa                                                                                 >  1:3114446/1‑150 (MQ=255)
                                                             aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaa                                                                                 >  1:3114340/1‑150 (MQ=255)
                                                             aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaa                                                                                 <  2:644787/150‑1 (MQ=255)
                                                               aaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAg                                                                               <  2:3402290/150‑1 (MQ=255)
                                                                   cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAAc                                                                           <  2:2061031/150‑1 (MQ=255)
                                                                   cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAAc                                                                           >  1:3626744/1‑150 (MQ=255)
                                                                       aTGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                                                                       <  1:2015654/150‑1 (MQ=255)
                                                                        tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa                                                                      <  1:3264272/150‑1 (MQ=255)
                                                                          aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATc                                                                    >  1:4105772/1‑150 (MQ=255)
                                                                                tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACg                                                              >  1:4632342/1‑150 (MQ=255)
                                                                                tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACg                                                              <  1:3558205/150‑1 (MQ=255)
                                                                                    atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGc                                                          >  2:13492/1‑150 (MQ=255)
                                                                                    atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGc                                                          >  1:620075/1‑150 (MQ=255)
                                                                                    atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGc                                                          >  2:3889539/1‑150 (MQ=255)
                                                                                      atGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCt                                                        <  1:3007498/150‑1 (MQ=255)
                                                                                        gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTAt                                                      >  2:204278/1‑150 (MQ=255)
                                                                                           tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCa                                                   >  2:3298389/1‑150 (MQ=255)
                                                                                              gTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCa                                                >  2:2977135/1‑150 (MQ=255)
                                                                                                 tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGa                                             >  2:946026/1‑150 (MQ=255)
                                                                                                 tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGa                                             >  2:947082/1‑150 (MQ=255)
                                                                                                  aTTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGAc                                            >  2:1787568/1‑150 (MQ=255)
                                                                                                   ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAGAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGAct                                           <  2:4104114/150‑1 (MQ=255)
                                                                                                   ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGAct                                           <  2:4538285/150‑1 (MQ=255)
                                                                                                     cACCCGCTGTCTGCGGCGGACGAACCGGCCGGGGAAGCGAAAAGAAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTAAGAAGTGTGCGTCGTCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGActct                                         >  1:1902825/1‑150 (MQ=255)
                                                                                                      aCCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTg                                        <  1:1365209/150‑1 (MQ=255)
                                                                                                      aCCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTg                                        <  1:4236810/150‑1 (MQ=255)
                                                                                                         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCg                                     >  2:4164819/1‑150 (MQ=255)
                                                                                                            tGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTg                                  >  2:1595882/1‑150 (MQ=255)
                                                                                                                tGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACt                              <  1:2313865/150‑1 (MQ=255)
                                                                                                                      ggTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCa                        >  2:1200544/1‑150 (MQ=255)
                                                                                                                       gTCGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAg                       >  1:2017524/1‑150 (MQ=255)
                                                                                                                         cGATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGcc                     <  1:737506/150‑1 (MQ=255)
                                                                                                                          gATCCGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCg                    >  2:1730076/1‑150 (MQ=255)
                                                                                                                             ccGGCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTg                 <  2:2508283/150‑1 (MQ=255)
                                                                                                                               ggCCTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGa               >  1:3817979/1‑150 (MQ=255)
                                                                                                                                 ccTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtg             >  2:607237/1‑150 (MQ=255)
                                                                                                                                  cTGGGAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtgt            <  2:2991701/150‑1 (MQ=255)
                                                                                                                                     ggAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGaa         <  2:1014414/150‑1 (MQ=255)
                                                                                                                                      gAAGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGaaa        <  2:4515663/150‑1 (MQ=255)
                                                                                                                                       aaGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAg       >  2:4565402/1‑150 (MQ=255)
                                                                                                                                        aGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGACCTGGCGCAGCCGCTGGATGTGAAAGa      >  1:1291607/1‑150 (MQ=255)
                                                                                                                                        aGCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGa      >  2:3123489/1‑150 (MQ=255)
                                                                                                                                         gCGAAAAGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGAc     >  1:1141303/1‑150 (MQ=255)
                                                                                                                                            aaaaGCAACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTgg  >  2:2301525/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGG  >  CP009273/1277640‑1277929

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: ‑