Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 329,887 G→T 6.6% I86I (ATC→ATA yahB ← putative DNA‑bindng transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273329,8870GT 6.6% 669.7 / 12.9 194I86I (ATC→ATAyahBputative DNA‑bindng transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (95/80);  new base T (6/13);  total (101/93)
Fisher's exact test for biased strand distribution p-value = 8.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

CCGTTATACACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAG  >  CP009273/329738‑330034
                                                                                                                                                     |                                                                                                                                                   
ccGTTATACACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCg                                                                                                                                                     >  1:600664/1‑150 (MQ=255)
 cGTTATACACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGa                                                                                                                                                    >  2:3133467/1‑150 (MQ=255)
   ttATACACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGAtt                                                                                                                                                  >  2:4234899/1‑150 (MQ=255)
    tataCACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATACAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTc                                                                                                                                                 <  2:4054302/150‑1 (MQ=255)
     ataCACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCt                                                                                                                                                <  1:4283117/150‑1 (MQ=255)
      taCACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTg                                                                                                                                               <  1:4472164/150‑1 (MQ=255)
       acacTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGc                                                                                                                                              <  2:52400/150‑1 (MQ=255)
       acacTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAAATGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGc                                                                                                                                              <  1:1330625/150‑1 (MQ=255)
         acTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGt                                                                                                                                            >  2:1442550/1‑150 (MQ=255)
         acTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATGCTGCGt                                                                                                                                            >  2:2916515/1‑150 (MQ=255)
          cTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGtt                                                                                                                                           <  1:2991126/150‑1 (MQ=255)
          cTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGtt                                                                                                                                           <  1:2141199/150‑1 (MQ=255)
            tCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTg                                                                                                                                         <  2:3474424/150‑1 (MQ=255)
                tggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATgg                                                                                                                                     >  2:3904175/1‑150 (MQ=255)
                tggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATgg                                                                                                                                     <  1:3055971/150‑1 (MQ=255)
                tggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATgg                                                                                                                                     >  2:194290/1‑150 (MQ=255)
                tggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATgg                                                                                                                                     <  2:3960437/150‑1 (MQ=255)
                 ggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGc                                                                                                                                    <  1:2073263/150‑1 (MQ=255)
                 ggtgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGc                                                                                                                                    >  1:3623113/1‑150 (MQ=255)
                   tgACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCAt                                                                                                                                  <  2:2050213/150‑1 (MQ=255)
                    gACGGTAATCGGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATa                                                                                                                                 <  2:4002608/150‑1 (MQ=255)
                     aCGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATaa                                                                                                                                <  1:1694117/150‑1 (MQ=255)
                        gTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATaaaaa                                                                                                                             >  2:1832298/1‑150 (MQ=255)
                          aaTCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATc                                                                                                                           <  1:3738274/150‑1 (MQ=255)
                              tGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTc                                                                                                                       >  1:3021461/1‑150 (MQ=255)
                               ggCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTca                                                                                                                      <  1:1282153/150‑1 (MQ=255)
                                gCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTcag                                                                                                                     <  2:3346905/150‑1 (MQ=255)
                                  aGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTcagca                                                                                                                   <  1:4717025/150‑1 (MQ=255)
                                      atgAAACTTCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCATACGCGCTTCTATGCCCTGCTAAATAGTATCGATTCTGCGTTTGATGTCATAAAAATCATTCAGCAGGTc                                                                                                               <  1:256146/148‑1 (MQ=255)
                                      aGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCATACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTc                                                                                                               <  1:662025/150‑1 (MQ=255)
                                      aGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTc                                                                                                               <  2:2003826/150‑1 (MQ=255)
                                         aaaCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGt                                                                                                            <  2:1478057/150‑1 (MQ=255)
                                         aaaCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGt                                                                                                            >  2:3559156/1‑150 (MQ=255)
                                         aaaCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCATATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGt                                                                                                            >  1:3290769/1‑150 (MQ=255)
                                           aCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAg                                                                                                          <  1:3396585/150‑1 (MQ=255)
                                           aCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAg                                                                                                          <  1:80369/150‑1 (MQ=255)
                                           aCTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAg                                                                                                          <  1:2320992/150‑1 (MQ=255)
                                              gCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCtt                                                                                                       >  1:1620280/1‑150 (MQ=255)
                                              gCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCtt                                                                                                       <  2:1090146/150‑1 (MQ=255)
                                                tttttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTg                                                                                                     >  2:2620688/1‑150 (MQ=255)
                                                 ttttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGc                                                                                                    >  1:2126247/1‑150 (MQ=255)
                                                   ttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGTGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGAgg                                                                                                  <  2:2354399/150‑1 (MQ=255)
                                                   ttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCgg                                                                                                  >  1:1563013/1‑150 (MQ=255)
                                                     cAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCTATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGaa                                                                                                <  1:1846759/150‑1 (MQ=255)
                                                     cAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCTATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGaa                                                                                                <  1:374035/150‑1 (MQ=255)
                                                         aCCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGaaaaaa                                                                                            <  2:3257844/150‑1 (MQ=255)
                                                              tAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACg                                                                                       <  2:2969244/150‑1 (MQ=255)
                                                              tAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACg                                                                                       >  1:3302579/1‑150 (MQ=255)
                                                               aaCAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGc                                                                                      >  2:1015050/1‑150 (MQ=255)
                                                                aCAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCTAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGcc                                                                                     <  1:4426269/150‑1 (MQ=255)
                                                                 cAAGCGTGCGGTGTGTTTGGGCGTATAAAGCTGCTGATTAATACAAATGCGCACACGCGATTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCTGGTCGGTAGCTTTGCGGAAAAAATAACGccc                                                                                    <  1:2790823/150‑1 (MQ=255)
                                                                 cAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGccc                                                                                    >  1:107567/1‑150 (MQ=255)
                                                                  aaGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCTATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCg                                                                                   >  2:2017946/1‑150 (MQ=255)
                                                                       tGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTc                                                                              >  2:4133904/1‑150 (MQ=255)
                                                                       tGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTc                                                                              >  1:2774932/1‑150 (MQ=255)
                                                                          ggTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGt                                                                           >  2:4697911/1‑150 (MQ=255)
                                                                           gtgtgtTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGtt                                                                          >  2:4201086/1‑150 (MQ=255)
                                                                            tgtgtTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTa                                                                         >  2:2293716/1‑150 (MQ=255)
                                                                                  tGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCaa                                                                   >  1:4739743/1‑150 (MQ=255)
                                                                                   gggCGTATAAAGCAGATTATTAATACAAATGCGCACACGCGATTATATGCCCTGCGAAATATTATCGATTCTGCTTTTGTTGGCTTAAAAATCATTCCGCAGGTCGTTTGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAt                                                                  <  2:1849939/150‑1 (MQ=255)
                                                                                    ggCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCTATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCaaag                                                                 >  1:3910667/1‑148 (MQ=255)
                                                                                    ggCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAtg                                                                 >  1:4318205/1‑150 (MQ=255)
                                                                                       gTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAtgctg                                                              >  2:2483377/1‑150 (MQ=255)
                                                                                          taAATCAGATTATAAAAAAAAATGCGCACCCGCGCTTCTATGCCCTGCTAAATAGTATCTATTCTGCTTTTTATGGAAAAAAAATCATTCAGCTGGTCGTTAGCTTTGCGTAAAAAAAAACGCCCGGACTCCGTTAACTCAATGCTGCgg                                                           <  2:3010999/150‑1 (MQ=255)
                                                                                          taAAGCCGCTTTTTAATACAAAGGCGCACACGCGCTTCTTTGCCCTGCGAAATTTTATCGATTCTGGGTTTGATGGAATAAAAATCATCCAGCTGGTCGGTAGCTTTGCGTAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCgg                                                           <  2:621510/150‑1 (MQ=255)
                                                                                           aaaGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCggg                                                          >  1:458909/1‑150 (MQ=255)
                                                                                           aaaGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCggg                                                          >  1:2851567/1‑150 (MQ=255)
                                                                                            acgcagcTTATTAATACAAAGGCGCACCCGCGATTCTAGGCCCTGCTAATTTGTATCTTTTCTGCTTTTTATGGAATAAAAATCATTCAGCAGGGCGGTAGCTTTGGTTAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGt                                                         <  2:2126721/148‑1 (MQ=255)
                                                                                            aaGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGt                                                         >  1:385759/1‑150 (MQ=255)
                                                                                                agcTTATAAATACAACTGCGCACCCGGGATTCTAGGCCCTGCTAATTTGTATCTTTTCGGCTTTTGATGGCATAAAAATCTTTCCGCTGGTGGGTTGTTTTGGGTAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTg                                                     <  1:2263627/150‑1 (MQ=255)
                                                                                                agcTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTg                                                     >  1:3142492/1‑150 (MQ=255)
                                                                                                 gcTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGc                                                    >  2:4015040/1‑150 (MQ=255)
                                                                                                  cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCg                                                   >  1:3327226/1‑150 (MQ=255)
                                                                                                  cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCg                                                   >  1:3327223/1‑150 (MQ=255)
                                                                                                   tGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGa                                                  >  2:143025/1‑150 (MQ=255)
                                                                                                    gATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAg                                                 >  2:4458614/1‑150 (MQ=255)
                                                                                                       tAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGa                                              >  1:4101741/1‑150 (MQ=255)
                                                                                                       tAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGa                                              >  2:1333813/1‑150 (MQ=255)
                                                                                                          tACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACa                                           >  2:1122395/1‑150 (MQ=255)
                                                                                                          tACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACa                                           >  2:3527025/1‑150 (MQ=255)
                                                                                                          caCAACTGCGCACACGCGATTATAGGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCCGGTCGGTAGCTTTGCGTAAAAAATCACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACa                                           <  2:1633803/149‑1 (MQ=255)
                                                                                                            cAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGc                                         >  1:4161197/1‑150 (MQ=255)
                                                                                                                tGCGCACCCGCGGTTATTTGCCCTGCTAATTTTTATCGTTTGTGCTTTTGTTGGAATAAAAATCATTCAGCTGGTCGTTAGCTTGGCGTAAAAAATCACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGcacca                                     <  2:3615289/150‑1 (MQ=255)
                                                                                                                   gcACACGCGATTATATGCCCTGCGAAATATTATCTTTTCTGGTTTTGTTGGAATAAAAATCATTCAGCTGGTCGGTAGATTTGCGTAAAAAATCACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACAt                                  <  1:1132128/150‑1 (MQ=255)
                                                                                                                    cacaCGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATc                                 >  1:3507527/1‑150 (MQ=255)
                                                                                                                       acGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAg                              >  2:374818/1‑150 (MQ=255)
                                                                                                                       acGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAg                              >  2:2765795/1‑150 (MQ=255)
                                                                                                                       acGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAg                              >  1:554631/1‑150 (MQ=255)
                                                                                                                           gcTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGcccc                          >  1:1987969/1‑150 (MQ=255)
                                                                                                                            cTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCg                         >  1:634850/1‑150 (MQ=255)
                                                                                                                                tttGCCCTGCGAAATTGTATCTTTTCTGCGTTTGATGGCATAAAAATCATTCAGCTGGTCGGTAGCTTTGCGTAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGtctc                     <  2:2616446/148‑1 (MQ=255)
                                                                                                                                 aTGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTcc                    >  2:1822812/1‑150 (MQ=255)
                                                                                                                                   gCCCTGCTAAATAGTATCTTTTCTGCGTTTGATGGAATAAAAATCATTCCGCTGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCAt                  <  2:600634/150‑1 (MQ=255)
                                                                                                                                    cccTGCGAAATAGTATCGATTCTGCGTTTGTTGGCATAAAAATCATTCAGCAGGTCGGTTGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATa                 <  2:377096/150‑1 (MQ=255)
                                                                                                                                          gAATTAGTATCGATTCTGCTTTTGATGGAATAAAAATCATTCAGCTGGTGGGTAGCTTTGGGTAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTa           <  2:3361272/150‑1 (MQ=255)
                                                                                                                                          gAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTa           >  2:4327220/1‑150 (MQ=255)
                                                                                                                                          gAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTa           >  2:4784144/1‑150 (MQ=255)
                                                                                                                                          gAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTa           >  2:669917/1‑150 (MQ=255)
                                                                                                                                                 tATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGgtgt    >  1:2492347/1‑150 (MQ=255)
                                                                                                                                                   tCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAg  >  1:3899947/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTTATACACTTCGGTGGTGACGGTAATCTGGCAGGTAGGAAACTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAG  >  CP009273/329738‑330034

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: