Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 87,619 G→A 5.8% R34Q (CGA→CAA)  ftsL → membrane bound cell division protein at septum containing leucine zipper motif

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927387,6190GA 5.8% 1319.6 / 42.4 361R34Q (CGA→CAA) ftsLmembrane bound cell division protein at septum containing leucine zipper motif
Reads supporting (aligned to +/- strand):  ref base G (181/159);  new base A (12/9);  total (193/168)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 361 total aligned reads displayed.

CCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAATGGCGC  >  CP009273/87472‑87764
                                                                                                                                                   |                                                                                                                                                 
ccTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGa                                                                                                                                                  <  2:3151105/149‑1 (MQ=255)
 cTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAt                                                                                                                                                 <  1:1030478/149‑1 (MQ=255)
  tCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAtt                                                                                                                                                <  1:4892741/149‑1 (MQ=255)
   cGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAttt                                                                                                                                               <  1:4170053/149‑1 (MQ=255)
   cGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAttt                                                                                                                                               <  1:1664043/149‑1 (MQ=255)
    gTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTg                                                                                                                                              <  1:3093553/149‑1 (MQ=255)
    gTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTg                                                                                                                                              <  1:3839237/149‑1 (MQ=255)
       cccGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGa                                                                                                                                           <  2:4133949/149‑1 (MQ=255)
        ccGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGaa                                                                                                                                          <  2:1830388/149‑1 (MQ=255)
          gTAGTTCAGTTCTGCTTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGc                                                                                                                                        <  2:2951062/149‑1 (MQ=255)
          gTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGc                                                                                                                                        >  2:4689329/1‑149 (MQ=255)
             gttcagttcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGc                                                                                                                                     >  2:3279152/1‑149 (MQ=255)
                 agttcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                                                                                                                                 <  1:2315981/149‑1 (MQ=255)
                 agttcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                                                                                                                                 >  1:3067836/1‑149 (MQ=255)
                  gttcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCActc                                                                                                                                <  2:689148/149‑1 (MQ=255)
                      tGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGcc                                                                                                                            <  2:2642735/149‑1 (MQ=255)
                      tGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGcc                                                                                                                            <  1:2709069/149‑1 (MQ=255)
                      tGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGcc                                                                                                                            <  1:536781/149‑1 (MQ=255)
                       gcttATTGCAGAGATGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCt                                                                                                                           <  2:1471817/146‑1 (MQ=255)
                         gTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGCAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGt                                                                                                                         >  2:3341017/1‑149 (MQ=255)
                         gTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGt                                                                                                                         >  2:4829946/1‑149 (MQ=255)
                           aTTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTc                                                                                                                       <  1:74700/149‑1 (MQ=255)
                           aTTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTc                                                                                                                       <  1:1119230/149‑1 (MQ=255)
                            ttGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCa                                                                                                                      <  2:3861220/149‑1 (MQ=255)
                                agagagGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTg                                                                                                                  <  2:537239/149‑1 (MQ=255)
                                agagagGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTGGGAAGCTGCCACTCTGCCTGTTCATTTg                                                                                                                  >  2:865796/1‑149 (MQ=255)
                                   gagGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCat                                                                                                               <  1:4077215/149‑1 (MQ=255)
                                   gagGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCat                                                                                                               <  1:459725/149‑1 (MQ=255)
                                    agGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCatt                                                                                                              >  1:4730308/1‑149 (MQ=255)
                                     ggACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCatta                                                                                                             >  1:3073114/1‑149 (MQ=255)
                                        cGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTAttt                                                                                                          <  2:2174907/149‑1 (MQ=255)
                                           aTGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGa                                                                                                       <  1:956093/149‑1 (MQ=255)
                                             gCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcg                                                                                                     >  1:832448/1‑149 (MQ=255)
                                             gCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcg                                                                                                     >  1:1994635/1‑149 (MQ=255)
                                             gCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcg                                                                                                     >  2:4788113/1‑149 (MQ=255)
                                                tGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcggcg                                                                                                  >  1:1417553/1‑149 (MQ=255)
                                                 gATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcggcgg                                                                                                 >  1:3007801/1‑149 (MQ=255)
                                                  aTCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGt                                                                                                <  1:287732/149‑1 (MQ=255)
                                                    cagcagAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGa                                                                                              <  1:4339692/149‑1 (MQ=255)
                                                    cagcagAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGa                                                                                              <  2:2483768/149‑1 (MQ=255)
                                                     agcagAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGTCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGAc                                                                                             >  2:5063378/1‑149 (MQ=255)
                                                           gTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGt                                                                                       >  2:2590103/1‑149 (MQ=255)
                                                           gTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGt                                                                                       >  2:1670475/1‑149 (MQ=255)
                                                            tGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTa                                                                                      >  2:1180541/1‑149 (MQ=255)
                                                               cAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAAcc                                                                                   >  1:4410984/1‑149 (MQ=255)
                                                                 gAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCAc                                                                                 >  2:3622604/1‑149 (MQ=255)
                                                                  aaGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACg                                                                                >  2:3828910/1‑149 (MQ=255)
                                                                   aGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACgg                                                                               >  2:4506487/1‑149 (MQ=255)
                                                                    gCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGgc                                                                              >  1:3274659/1‑149 (MQ=255)
                                                                     ctctAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGgcg                                                                             >  2:603401/1‑149 (MQ=255)
                                                                      tctAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGgcgc                                                                            <  1:3477662/149‑1 (MQ=255)
                                                                         aaGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCAcc                                                                         <  2:1974559/149‑1 (MQ=255)
                                                                          aGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCa                                                                        >  1:358222/1‑149 (MQ=255)
                                                                             aaaGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATAc                                                                     >  2:4644287/1‑149 (MQ=255)
                                                                              aaGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATAcc                                                                    >  2:4936997/1‑149 (MQ=255)
                                                                                gTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCg                                                                  >  1:2085178/1‑149 (MQ=255)
                                                                                 ttAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGt                                                                 >  1:9422/1‑149 (MQ=255)
                                                                                  tAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGtt                                                                >  2:1013578/1‑149 (MQ=255)
                                                                                   aaaGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGttt                                                               >  1:3495085/1‑149 (MQ=255)
                                                                                   aaaGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGttt                                                               >  1:5036874/1‑149 (MQ=255)
                                                                                    aaGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTa                                                              <  1:1974722/149‑1 (MQ=255)
                                                                                      ggATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACt                                                            >  2:2838679/1‑149 (MQ=255)
                                                                                       gATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTg                                                           >  2:3795954/1‑149 (MQ=255)
                                                                                       gATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTg                                                           <  2:3670327/149‑1 (MQ=255)
                                                                                        aTCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGa                                                          >  2:319618/1‑149 (MQ=255)
                                                                                          cGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGAcc                                                        >  1:3411262/1‑149 (MQ=255)
                                                                                          cGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGAcc                                                        >  2:2885951/1‑149 (MQ=255)
                                                                                          cGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGAcc                                                        <  1:1302334/149‑1 (MQ=255)
                                                                                            aTGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGc                                                      <  1:376553/149‑1 (MQ=255)
                                                                                              gggAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTc                                                    >  2:4396852/1‑149 (MQ=255)
                                                                                                     cACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGa                                             <  2:4216775/149‑1 (MQ=255)
                                                                                                      aCGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGaa                                            <  2:4347698/149‑1 (MQ=255)
                                                                                                      aCGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGaa                                            >  1:2503358/1‑149 (MQ=255)
                                                                                                         aGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGaacaa                                         <  2:3298816/149‑1 (MQ=255)
                                                                                                         aGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGaacaa                                         >  1:3188856/1‑149 (MQ=255)
                                                                                                             ccATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTgg                                     <  2:3389273/149‑1 (MQ=255)
                                                                                                              cATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGt                                    <  2:2041062/149‑1 (MQ=255)
                                                                                                               aTGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTg                                   <  2:2367049/149‑1 (MQ=255)
                                                                                                                tGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCAATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGc                                  >  1:2570092/1‑149 (MQ=255)
                                                                                                                 gCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCt                                 >  1:657534/1‑149 (MQ=255)
                                                                                                                       ccTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCg                           <  2:4286638/149‑1 (MQ=255)
                                                                                                                          ggTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCgaga                        >  2:1291821/1‑149 (MQ=255)
                                                                                                                          ggTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCgaga                        >  1:559581/1‑149 (MQ=255)
                                                                                                                             gtTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGc                     >  2:3877035/1‑149 (MQ=255)
                                                                                                                              ttATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCt                    <  2:1461631/149‑1 (MQ=255)
                                                                                                                               tATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCtt                   <  1:1468837/149‑1 (MQ=255)
                                                                                                                               tATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCtt                   >  1:2673076/1‑149 (MQ=255)
                                                                                                                                aTCGGTGACGATCTTTTGCGATGTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCttt                  <  2:3943313/149‑1 (MQ=255)
                                                                                                                                   ggTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGa               <  1:603401/149‑1 (MQ=255)
                                                                                                                                   ggTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGa               <  1:3321370/149‑1 (MQ=255)
                                                                                                                                    gTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGAc              >  1:1522295/1‑149 (MQ=255)
                                                                                                                                    gTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGAc              >  1:1753136/1‑149 (MQ=255)
                                                                                                                                      gACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACAt            >  2:1588970/1‑149 (MQ=255)
                                                                                                                                       aCGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACAtt           >  2:3977517/1‑149 (MQ=255)
                                                                                                                                          aTCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGaa        <  1:4498311/149‑1 (MQ=255)
                                                                                                                                          aTCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGaa        >  2:547275/1‑149 (MQ=255)
                                                                                                                                           tCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAAt       <  1:1085246/149‑1 (MQ=255)
                                                                                                                                            cTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAATg      <  2:3128573/149‑1 (MQ=255)
                                                                                                                                            cTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAATg      >  2:4055579/1‑149 (MQ=255)
                                                                                                                                             ttttGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAATgg     >  2:4184499/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAACTGGTGCTGGAGCGAGATGCTTTAGACATTGAATGGCGC  >  CP009273/87472‑87764

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: