Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,488,966 G→C 100% A533P (GCC→CCC)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,488,9660GC97.1% 1077.4 / 19.4 273A533P (GCC→CCC) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base C (148/117);  total (151/122)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 278 total aligned reads displayed.

AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488819‑3489112
                                                                                                                                                   |                                                                                                                                                  
aacCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGccc                                                                                                                                                  <  1:2433667/150‑1 (MQ=255)
 acCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcct                                                                                                                                                 <  2:679413/150‑1 (MQ=255)
       gCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCctgct                                                                                                                                           <  1:4148994/150‑1 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  2:3860585/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:2813455/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:232083/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:1802697/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:4098450/1‑150 (MQ=255)
         gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  1:1966740/1‑150 (MQ=255)
          aaCACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGa                                                                                                                                        >  1:3630844/1‑150 (MQ=255)
          aaCACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGa                                                                                                                                        >  2:4469591/1‑150 (MQ=255)
           acacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGaa                                                                                                                                       <  1:402972/150‑1 (MQ=255)
           acacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGATTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGaa                                                                                                                                       >  1:3579209/1‑150 (MQ=255)
            cacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAAc                                                                                                                                      <  2:477866/150‑1 (MQ=255)
              cTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGa                                                                                                                                    >  2:20960/1‑150 (MQ=255)
               tttGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGaa                                                                                                                                   >  1:3684262/1‑150 (MQ=255)
                  gCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAc                                                                                                                                <  1:1612013/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAGGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  1:2514851/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  2:713246/150‑1 (MQ=255)
                   cTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACa                                                                                                                               <  1:3332563/150‑1 (MQ=255)
                    ttATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAt                                                                                                                              >  1:1506632/1‑150 (MQ=255)
                       ttCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCt                                                                                                                           <  1:833990/150‑1 (MQ=255)
                         cgcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGa                                                                                                                         <  1:919713/150‑1 (MQ=255)
                           cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       <  2:4547386/150‑1 (MQ=255)
                           cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       <  1:1072901/150‑1 (MQ=255)
                            gTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGc                                                                                                                      >  2:342167/1‑150 (MQ=255)
                            gTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGc                                                                                                                      <  1:3208146/150‑1 (MQ=255)
                                 gAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGaa                                                                                                                 >  2:1217043/1‑150 (MQ=255)
                                  aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                <  1:3172348/150‑1 (MQ=255)
                                  aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                <  2:826861/150‑1 (MQ=255)
                                      gTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAc                                                                                                            <  1:837533/150‑1 (MQ=255)
                                      gTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATAc                                                                                                            >  1:447020/1‑150 (MQ=255)
                                       ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                           <  1:1729099/150‑1 (MQ=255)
                                       ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                           >  1:1203468/1‑150 (MQ=255)
                                        tCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTc                                                                                                          <  1:523654/150‑1 (MQ=255)
                                           aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGc                                                                                                       >  1:4449840/1‑150 (MQ=255)
                                            tttGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCt                                                                                                      <  1:2066892/150‑1 (MQ=255)
                                                      gTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCt                                                                                            >  1:3098682/1‑150 (MQ=255)
                                                         gAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAAc                                                                                         >  1:726614/1‑150 (MQ=255)
                                                           aGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGa                                                                                       >  1:86180/1‑150 (MQ=255)
                                                            ggCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGAc                                                                                      >  2:2648039/1‑150 (MQ=255)
                                                              cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                    >  1:1863662/1‑150 (MQ=255)
                                                              cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                    >  2:312015/1‑150 (MQ=255)
                                                                 tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                 >  2:1666229/1‑150 (MQ=255)
                                                                 tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                 <  2:4405814/150‑1 (MQ=255)
                                                                  aaaaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTc                                                                                >  2:2127927/1‑150 (MQ=255)
                                                                    aaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCt                                                                              <  2:3675124/150‑1 (MQ=255)
                                                                     aCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTc                                                                             >  2:1364810/1‑150 (MQ=255)
                                                                       cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                           >  1:3747224/1‑150 (MQ=255)
                                                                        ttCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCt                                                                          <  2:1823363/150‑1 (MQ=255)
                                                                          cGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTaa                                                                        >  2:2278307/1‑150 (MQ=255)
                                                                            aagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACa                                                                      <  1:4365764/150‑1 (MQ=255)
                                                                             agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                     >  2:731549/1‑150 (MQ=255)
                                                                                    ctgGCGAAACACGGCATAGGCTACGGTTGTGTAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                              <  2:2435473/150‑1 (MQ=255)
                                                                                    ctgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                              <  2:1752238/150‑1 (MQ=255)
                                                                                         gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACgg                                                         >  2:1578391/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  2:2685618/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  1:4061513/1‑150 (MQ=255)
                                                                                          aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  1:2176407/1‑150 (MQ=255)
                                                                                           aaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCa                                                       >  2:168466/1‑150 (MQ=255)
                                                                                            acacGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCAc                                                      <  2:4126391/150‑1 (MQ=255)
                                                                                                ggCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTAc                                                  >  2:1946021/1‑150 (MQ=255)
                                                                                                 gCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACt                                                 <  1:521543/150‑1 (MQ=255)
                                                                                                 gCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACt                                                 >  2:2917846/1‑150 (MQ=255)
                                                                                                   aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                               >  2:3942265/1‑150 (MQ=255)
                                                                                                      ggCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTg                                            <  2:2822220/150‑1 (MQ=255)
                                                                                                          aCGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATcc                                        >  2:3753216/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  1:2135776/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  2:1435729/1‑150 (MQ=255)
                                                                                                               tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   >  1:3659618/1‑150 (MQ=255)
                                                                                                                gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGtt                                  <  1:2218873/150‑1 (MQ=255)
                                                                                                                gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGtt                                  <  1:3864080/150‑1 (MQ=255)
                                                                                                                 tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                 <  1:1476109/150‑1 (MQ=255)
                                                                                                                 tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                 >  1:101173/1‑150 (MQ=255)
                                                                                                                  gAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtct                                <  2:748187/150‑1 (MQ=255)
                                                                                                                    aGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATCCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTcc                              >  2:4050618/1‑150 (MQ=255)
                                                                                                                     gTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCg                             <  2:4529943/150‑1 (MQ=255)
                                                                                                                       ttGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCggg                           <  2:18938/150‑1 (MQ=255)
                                                                                                                          tAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCgg                        >  1:1384141/1‑150 (MQ=255)
                                                                                                                          tAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCgg                        >  1:3059146/1‑150 (MQ=255)
                                                                                                                           aaaCCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGt                       >  1:1824304/1‑150 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:1090881/150‑1 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:859259/150‑1 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     >  2:3034191/1‑150 (MQ=255)
                                                                                                                             accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACGTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     <  2:4268086/150‑1 (MQ=255)
                                                                                                                                 accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                 <  2:1431907/150‑1 (MQ=255)
                                                                                                                                 accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                 <  2:147448/150‑1 (MQ=255)
                                                                                                                                  ccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCa                >  2:1284641/1‑150 (MQ=255)
                                                                                                                                   cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc               <  2:1966740/150‑1 (MQ=255)
                                                                                                                                    gTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcc              >  2:4528594/1‑150 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGACTACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            >  2:2690137/1‑150 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:3359866/150‑1 (MQ=255)
                                                                                                                                      cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:3630844/150‑1 (MQ=255)
                                                                                                                                        gTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCgg          >  2:261273/1‑150 (MQ=255)
                                                                                                                                            gctgctGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAgg      <  1:1257343/150‑1 (MQ=255)
                                                                                                                                            gctgctGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAgg      <  2:3592084/150‑1 (MQ=255)
                                                                                                                                              tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc    <  2:3101088/150‑1 (MQ=255)
                                                                                                                                              tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc    >  1:4224619/1‑150 (MQ=255)
                                                                                                                                               gctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCt   >  2:1596207/1‑150 (MQ=255)
                                                                                                                                                ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg  >  2:898033/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488819‑3489112

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: