Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,067,754 T→A 70.0% H189L (CAC→CTC)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,067,7540TA70.0% 313.5 / 172.6 199H189L (CAC→CTC) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (38/21);  new base A (85/54);  total (124/75)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

TTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTA  >  CP009273/1067609‑1067897
                                                                                                                                                 |                                                                                                                                               
ttGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGgaga                                                                                                                                               <  1:755816/148‑1 (MQ=255)
 tGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGagag                                                                                                                                              <  2:1649555/148‑1 (MQ=255)
        gTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTc                                                                                                                                       <  2:602492/148‑1 (MQ=255)
         ttAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCt                                                                                                                                      >  2:2451971/1‑148 (MQ=255)
         ttAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCt                                                                                                                                      <  1:3092280/148‑1 (MQ=255)
         ttAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCt                                                                                                                                      <  2:247580/148‑1 (MQ=255)
              caacTTAAGCGATATGATCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGACGCTGACCCAGCCAAAAAATGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCAGGAGAGTTGCGTCATCAAg                                                                                                                                 >  2:593745/4‑148 (MQ=255)
                ccTTAAGCGATATGAGCACAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGca                                                                                                                               <  2:2087182/148‑1 (MQ=255)
                ccTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGca                                                                                                                               >  2:410857/1‑148 (MQ=255)
                  ttAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcacc                                                                                                                             >  2:199114/1‑148 (MQ=255)
                   tAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcacca                                                                                                                            >  1:2854098/1‑148 (MQ=255)
                   tAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGcacca                                                                                                                            <  1:343527/148‑1 (MQ=255)
                   tAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGcacca                                                                                                                            >  1:885256/1‑148 (MQ=255)
                    aaGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcaccac                                                                                                                           >  2:1866015/1‑148 (MQ=255)
                      gCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGc                                                                                                                         >  1:2927286/1‑148 (MQ=255)
                         atatGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGa                                                                                                                      >  2:1982445/1‑148 (MQ=255)
                          tatGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCCAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGaa                                                                                                                     >  2:3343783/1‑148 (MQ=255)
                          tatGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGaa                                                                                                                     >  2:302747/1‑148 (MQ=255)
                           atGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGaaa                                                                                                                    >  2:494653/1‑148 (MQ=255)
                             gAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATa                                                                                                                  >  1:1550184/1‑148 (MQ=255)
                             gAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATa                                                                                                                  >  1:2854019/1‑148 (MQ=255)
                             gAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATa                                                                                                                  >  1:2099815/1‑148 (MQ=255)
                              aGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATAc                                                                                                                 <  2:2165491/148‑1 (MQ=255)
                               gCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACt                                                                                                                <  2:2840596/148‑1 (MQ=255)
                               gCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACt                                                                                                                <  2:2841299/148‑1 (MQ=255)
                               gCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACt                                                                                                                >  1:961756/1‑148 (MQ=255)
                                 acaGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTcc                                                                                                              >  1:1682407/3‑148 (MQ=255)
                                  aaGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCa                                                                                                             <  1:2464801/148‑1 (MQ=255)
                                   aGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAg                                                                                                            >  1:1770055/1‑148 (MQ=255)
                                        gTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGaaaaa                                                                                                       <  2:2450247/148‑1 (MQ=255)
                                        gTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGaaaaa                                                                                                       <  1:3297321/148‑1 (MQ=255)
                                             agaAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCg                                                                                                  <  2:960369/148‑1 (MQ=255)
                                                   cgcgTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGt                                                                                            >  2:2323962/1‑148 (MQ=255)
                                                       tCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCg                                                                                        <  2:2540544/148‑1 (MQ=255)
                                                          cAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCg                                                                                     >  2:3369221/1‑148 (MQ=255)
                                                            gAATGTTTCGACGTCGCTGACCCAGCCAAACAAGGTGTCGATATTGCACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAt                                                                                   <  2:259643/148‑1 (MQ=255)
                                                             aaTGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAtt                                                                                  >  1:966932/1‑148 (MQ=255)
                                                             aaTGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAtt                                                                                  >  2:1574098/1‑148 (MQ=255)
                                                             aaTGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAtt                                                                                  >  2:575284/1‑148 (MQ=255)
                                                              aTGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTc                                                                                 >  1:1092298/1‑148 (MQ=255)
                                                                gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGa                                                                               <  1:302747/148‑1 (MQ=255)
                                                                gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGa                                                                               <  1:3407603/148‑1 (MQ=255)
                                                                gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGa                                                                               <  1:412697/148‑1 (MQ=255)
                                                                      aCGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGa                                                                         >  2:963539/1‑148 (MQ=255)
                                                                        gTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACg                                                                       <  2:3442896/148‑1 (MQ=255)
                                                                         tCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTACAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGt                                                                      >  1:3456393/1‑148 (MQ=255)
                                                                          cGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGtt                                                                     >  1:204909/1‑148 (MQ=255)
                                                                            cTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTgg                                                                   <  1:532908/148‑1 (MQ=255)
                                                                             tGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGt                                                                  <  2:1396537/148‑1 (MQ=255)
                                                                              gACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTACGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTa                                                                 <  2:853548/148‑1 (MQ=255)
                                                                                 ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCg                                                              <  2:676136/148‑1 (MQ=255)
                                                                                    gCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATg                                                           >  1:1870125/1‑148 (MQ=255)
                                                                                    gCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATg                                                           >  2:2311243/1‑148 (MQ=255)
                                                                                     ccAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGc                                                          >  2:3096356/1‑148 (MQ=255)
                                                                                        aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCgg                                                       >  2:3217402/1‑148 (MQ=255)
                                                                                        aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCgg                                                       >  2:863905/1‑148 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                      >  2:1264047/1‑148 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                      >  1:2108375/1‑148 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                      >  1:1473779/1‑148 (MQ=255)
                                                                                          aaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTg                                                     >  2:3399526/1‑148 (MQ=255)
                                                                                          aaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTg                                                     >  2:690592/1‑148 (MQ=255)
                                                                                            aGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaa                                                   >  2:1938575/1‑148 (MQ=255)
                                                                                            aGGTTTCGAGATTGCACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaa                                                   <  2:2968616/148‑1 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  1:1399859/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  1:3150790/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  1:1950410/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  <  1:1651865/148‑1 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  1:372150/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  1:372038/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGACGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                  >  2:768116/1‑148 (MQ=255)
                                                                                                ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAcc                                               >  2:1421524/1‑148 (MQ=255)
                                                                                                ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAcc                                               >  2:3051228/1‑148 (MQ=255)
                                                                                                ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAcc                                               >  1:1117864/1‑148 (MQ=255)
                                                                                                  cGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAg                                             <  1:667161/148‑1 (MQ=255)
                                                                                                   gATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGa                                            <  2:2045203/148‑1 (MQ=255)
                                                                                                        tGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                       <  1:2510292/148‑1 (MQ=255)
                                                                                                         gAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGt                                      >  1:2766375/1‑148 (MQ=255)
                                                                                                          aacaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGTGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                     >  1:2366326/1‑148 (MQ=255)
                                                                                                            caacGCTGCTTTCTTCGAAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGACGTAGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATTCCGGTGAAAACCAGATGGCGTAtt                                   <  2:42093/148‑1 (MQ=255)
                                                                                                             aacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTc                                  <  1:248064/148‑1 (MQ=255)
                                                                                                              acGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTcc                                 >  2:1067312/1‑148 (MQ=255)
                                                                                                              acGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTcc                                 >  1:1162672/1‑148 (MQ=255)
                                                                                                                 cgGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGATGAAAACCAGATGGCGTATTCcgcg                              >  2:2290756/1‑148 (MQ=255)
                                                                                                                   gCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGc                            >  2:962350/1‑148 (MQ=255)
                                                                                                                    cTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCt                           >  2:1633558/1‑148 (MQ=255)
                                                                                                                      ttCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTgc                         >  2:1364093/1‑148 (MQ=255)
                                                                                                                      ttCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTgc                         >  1:492114/1‑148 (MQ=255)
                                                                                                                      ttCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTgc                         >  1:3272374/1‑148 (MQ=255)
                                                                                                                          gcgcAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaa                     >  2:2490182/1‑148 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                    >  2:1964921/1‑148 (MQ=255)
                                                                                                                              aaaTTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                 <  1:806837/148‑1 (MQ=255)
                                                                                                                              aaaTTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                 <  2:96895/148‑1 (MQ=255)
                                                                                                                                aTTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCt               >  2:1912211/1‑148 (MQ=255)
                                                                                                                                aTTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCt               >  1:2174722/1‑148 (MQ=255)
                                                                                                                                   tCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTc            <  1:2663518/148‑1 (MQ=255)
                                                                                                                                    cGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTcc           >  1:1702067/1‑148 (MQ=255)
                                                                                                                                         ccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcgg      <  1:410857/148‑1 (MQ=255)
                                                                                                                                           gCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcggcg    >  1:2269626/1‑148 (MQ=255)
                                                                                                                                             cTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTa  <  1:727584/148‑1 (MQ=255)
                                                                                                                                             cTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTa  <  1:3377777/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTA  >  CP009273/1067609‑1067897

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: