Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 426,157 G→T 93.8% E33* (GAG→TAG)  yajD → HNH nuclease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273426,1570GT93.8% 507.4 / 16.0 212E33* (GAG→TAG) yajDHNH nuclease family protein
Reads supporting (aligned to +/- strand):  ref base G (14/2);  new base T (97/98);  total (112/100)
Fisher's exact test for biased strand distribution p-value = 3.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

ATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGAT  >  CP009273/426013‑426300
                                                                                                                                                |                                                                                                                                               
aTCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgt                                                                                                                                              <  1:2941039/148‑1 (MQ=255)
 tCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAggt                                                                                                                                             >  2:513876/1‑146 (MQ=255)
   ttGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgtttg                                                                                                                                           <  1:1613437/148‑1 (MQ=255)
    tGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgtttgt                                                                                                                                          <  2:3427317/148‑1 (MQ=255)
    tGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgtttgt                                                                                                                                          >  2:1027093/1‑148 (MQ=255)
     gggAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgtttgtt                                                                                                                                         <  1:3231003/148‑1 (MQ=255)
      ggAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgtttgttt                                                                                                                                        >  1:1090627/1‑148 (MQ=255)
      ggAAACTGAGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAACCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAttttgttt                                                                                                                                        <  1:78190/148‑1 (MQ=255)
            tGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCa                                                                                                                                  >  2:3163434/1‑148 (MQ=255)
                cgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCt                                                                                                                              <  2:3095522/148‑1 (MQ=255)
                     cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTg                                                                                                                         >  1:1133002/1‑148 (MQ=255)
                     cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTg                                                                                                                         >  1:408026/1‑148 (MQ=255)
                       cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGaa                                                                                                                       >  2:961327/1‑148 (MQ=255)
                       cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGaa                                                                                                                       >  2:1353800/1‑148 (MQ=255)
                       cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGaa                                                                                                                       >  2:1077579/1‑148 (MQ=255)
                          cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACtt                                                                                                                    >  1:974934/1‑148 (MQ=255)
                          cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACtt                                                                                                                    >  1:2373535/1‑148 (MQ=255)
                           aTTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTa                                                                                                                   <  1:1642815/148‑1 (MQ=255)
                           aTTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTa                                                                                                                   <  1:812508/148‑1 (MQ=255)
                              tCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCg                                                                                                                <  1:3253654/148‑1 (MQ=255)
                                aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGtt                                                                                                              >  2:1038555/1‑148 (MQ=255)
                                 tCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGGTAGTTTGTTTTTTCCAACCTGCGTGAACTTACCGTTc                                                                                                             >  1:114878/1‑148 (MQ=255)
                                 tCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTc                                                                                                             <  2:2313490/148‑1 (MQ=255)
                                    aGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTcacc                                                                                                          >  2:2600126/1‑148 (MQ=255)
                                     ggTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTcacca                                                                                                         >  1:1631328/1‑148 (MQ=255)
                                       tAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACcaca                                                                                                       >  2:1335991/1‑148 (MQ=255)
                                        aaCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAt                                                                                                      <  1:434107/148‑1 (MQ=255)
                                            tttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTCACCACATTGAt                                                                                                  >  1:2430948/1‑148 (MQ=255)
                                             ttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATc                                                                                                 <  2:134186/148‑1 (MQ=255)
                                                 tGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGa                                                                                             >  2:3202970/1‑148 (MQ=255)
                                                 tGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGa                                                                                             <  1:2904089/148‑1 (MQ=255)
                                                 tGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGa                                                                                             >  1:3216749/1‑148 (MQ=255)
                                                    cTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCa                                                                                          >  2:3136355/1‑148 (MQ=255)
                                                    cTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCa                                                                                          <  1:2659436/148‑1 (MQ=255)
                                                     tATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCAt                                                                                         >  2:1880501/1‑148 (MQ=255)
                                                      atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATa                                                                                        >  2:1759779/1‑148 (MQ=255)
                                                      atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATa                                                                                        <  1:2007287/148‑1 (MQ=255)
                                                      atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTAACCACATTGATCACGACCATa                                                                                        >  1:1656618/1‑148 (MQ=255)
                                                       tcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTCTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATAc                                                                                       >  2:2153188/1‑148 (MQ=255)
                                                       tcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATAc                                                                                       >  2:2804765/1‑148 (MQ=255)
                                                           cccAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTAACCACAATGATCACGaacaaaacaaa                                                                                   >  1:970473/1‑147 (MQ=255)
                                                            ccAAAAAACTATGCGCGGTTATAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATa                                                                                  <  2:1374115/148‑1 (MQ=255)
                                                             cAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATaa                                                                                 <  1:2323407/148‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGa                                                                           <  2:529416/148‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGa                                                                           <  1:641450/148‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAATGATCACGACCATAACAATAAACCGGa                                                                           >  1:1029499/1‑148 (MQ=255)
                                                                     tATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAg                                                                         <  1:1198273/148‑1 (MQ=255)
                                                                       tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTTAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGAt                                                                       <  2:3297850/148‑1 (MQ=255)
                                                                       tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACGTACCGGTCACCACATTGATCACGACCATACCCATAAACCGGAACAt                                                                       >  1:2793275/1‑148 (MQ=255)
                                                                         cgcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATgg                                                                     <  1:1749767/148‑1 (MQ=255)
                                                                         cgcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATgg                                                                     <  2:1090627/148‑1 (MQ=255)
                                                                         cgcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGAACATACCAATAACCCGGAAGATgg                                                                     <  2:2654038/148‑1 (MQ=255)
                                                                          gcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGc                                                                    <  1:3027468/148‑1 (MQ=255)
                                                                           cgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCa                                                                   >  2:924631/1‑148 (MQ=255)
                                                                           cgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCa                                                                   >  2:2724037/1‑148 (MQ=255)
                                                                              ttAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTa                                                                >  2:2231818/1‑148 (MQ=255)
                                                                               tAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTaa                                                               <  1:1061865/148‑1 (MQ=255)
                                                                                aGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAAc                                                              <  1:3564917/148‑1 (MQ=255)
                                                                                 gAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACt                                                             >  2:588720/1‑148 (MQ=255)
                                                                                    aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTggg                                                          >  2:3028338/1‑148 (MQ=255)
                                                                                    aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAATGATCACGACCATACCAATAACCCGGAAGATGGCAGTTACgggg                                                          >  1:3149977/1‑148 (MQ=255)
                                                                                    aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCAGCGTGAACATACCGTTAACCACATTGATAAAGACCATACCAATAAACCCGAAGATGGCAGTAACTggg                                                          >  1:1368859/1‑148 (MQ=255)
                                                                                      cGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGGTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTgggca                                                        >  1:2523677/1‑146 (MQ=255)
                                                                                       ggCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAt                                                       <  2:2405184/148‑1 (MQ=255)
                                                                                         cTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCCCGAGTTTGTTTATTAAAAACTGCGTGAACTTACCGTTCAACAAATATATCACGACCAAACCAAAAACACAGAAGAAGGCAGGAACTGGGAAatg                                                     >  2:3304508/1‑148 (MQ=255)
                                                                                            tCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAttgttg                                                  >  2:1546010/1‑148 (MQ=255)
                                                                                             cgcgAAAAAGCATTAAAACTCTATCCGTGGGTCTGCGGTCGCTGGTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATtgttgt                                                 <  1:546438/148‑1 (MQ=255)
                                                                                             cgcgAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATtgttgt                                                 >  2:675584/1‑148 (MQ=255)
                                                                                               cgAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTgtgt                                               >  2:1937227/1‑148 (MQ=255)
                                                                                                 aaaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtct                                             <  1:1027093/148‑1 (MQ=255)
                                                                                                 aaaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtct                                             <  1:1978523/148‑1 (MQ=255)
                                                                                                  aaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtctc                                            >  1:2572930/1‑148 (MQ=255)
                                                                                                  aaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtctc                                            <  1:1308432/148‑1 (MQ=255)
                                                                                                   aaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTctct                                           >  2:1150774/1‑148 (MQ=255)
                                                                                                   aaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGGTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTctct                                           <  2:1192828/148‑1 (MQ=255)
                                                                                                         ctAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGcc                                     >  1:579675/2‑148 (MQ=255)
                                                                                                              aaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAAACAGCGTGAACTTAACGTTCACCACATTGATAACGACCATACCAATAAAACGGAAGATTGCAGTAACTGGGAATTGTTGTGTGTCTATTTCCACGAt                                >  1:229983/1‑148 (MQ=255)
                                                                                                               aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                               >  2:948855/1‑148 (MQ=255)
                                                                                                               aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                               <  1:3439697/148‑1 (MQ=255)
                                                                                                               aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                               <  1:352362/148‑1 (MQ=255)
                                                                                                                tCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTAACCACATTTATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTTCCACGATCa                              >  2:2113820/1‑148 (MQ=255)
                                                                                                                   aTCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGa                           <  1:1353767/148‑1 (MQ=255)
                                                                                                                   aTCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGa                           <  2:452358/148‑1 (MQ=255)
                                                                                                                    tCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCAGCGTGAACATAACGTTCACCACAATGATAACGACCCAAACAAAAAACCGGAACATTGCAGTAACTGGGAATTGTTGTGTGTATATTTCCACCATCATGAg                          >  1:692792/1‑148 (MQ=255)
                                                                                                                         gggTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTc                     >  2:1378199/1‑148 (MQ=255)
                                                                                                                            tCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaa                  <  1:3474116/148‑1 (MQ=255)
                                                                                                                                cGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat              <  1:2643305/148‑1 (MQ=255)
                                                                                                                                cGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat              <  1:2885957/148‑1 (MQ=255)
                                                                                                                                cGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat              <  2:1629855/148‑1 (MQ=255)
                                                                                                                                cGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat              <  2:1631443/148‑1 (MQ=255)
                                                                                                                                cGGTCGCTGTTCCCGCTAGTTTGTTGATTCCAACCGGCGTGAACTTACCGTTCACCACAGTGATCACGACCAGACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatag              <  1:662359/148‑2 (MQ=255)
                                                                                                                                 ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata             <  2:229983/148‑1 (MQ=255)
                                                                                                                                   tCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc           >  2:2420629/1‑148 (MQ=255)
                                                                                                                                      cTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGaa        >  1:2495551/1‑148 (MQ=255)
                                                                                                                                        gTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGc      >  2:1526855/1‑148 (MQ=255)
                                                                                                                                         ttCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCg     >  2:3400750/1‑148 (MQ=255)
                                                                                                                                          tCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCgg    >  1:481221/1‑148 (MQ=255)
                                                                                                                                           cccGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGa   <  2:2077708/148‑1 (MQ=255)
                                                                                                                                            ccGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGAt  <  2:3216797/148‑1 (MQ=255)
                                                                                                                                            ccGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGAt  <  2:3216749/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGAT  >  CP009273/426013‑426300

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