Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,124,268 A→T 5.4% F31Y (TTC→TAC)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,124,2680AT 5.4% 824.8 / 25.6 259F31Y (TTC→TAC) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base A (129/115);  new base T (10/4);  total (140/119)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

ATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCA  >  CP009273/2124119‑2124415
                                                                                                                                                     |                                                                                                                                                   
 tttACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTACGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGaa                                                                                                                                                    <  1:2994434/150‑1 (MQ=255)
    aCTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACt                                                                                                                                                 >  1:2845119/1‑150 (MQ=255)
      tGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACttt                                                                                                                                               <  2:2561893/150‑1 (MQ=255)
      tGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCGAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTACTCAAGTGACGATGGTGAGGGGGAAACttt                                                                                                                                               >  2:5083126/1‑150 (MQ=255)
       ggCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTc                                                                                                                                              <  2:3840873/150‑1 (MQ=255)
         cATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGc                                                                                                                                            >  2:482215/1‑150 (MQ=255)
          aTCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGcc                                                                                                                                           >  1:3464012/1‑150 (MQ=255)
              aGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGc                                                                                                                                       >  1:1339553/1‑150 (MQ=255)
               tcGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGATGATGAGTTCTAAAAACAGCAGTGACAAGTCAGGTAAAAAAGTGACGATGGTGATGGGGAAACTGTCGCCAAGCt                                                                                                                                      >  2:2945528/2‑150 (MQ=255)
                cGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTg                                                                                                                                     >  2:307744/1‑150 (MQ=255)
                cGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTg                                                                                                                                     >  1:2220625/1‑150 (MQ=255)
                cGACTGACCAAACAGATAGTTAAACGCCGTGAGTCACGCGCCAATCCCTAAGGCAATCATCAATTTTTTAAAGCTGATGCTTTCTAAAAAAAGCAGTAACTAGACAGGTAACAAAGTGACGATGGTGATGGGGAAACATTATCCAATCTg                                                                                                                                     >  2:156718/1‑150 (MQ=255)
                 gACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTgg                                                                                                                                    >  2:3304215/1‑150 (MQ=255)
                  aCTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTggg                                                                                                                                   <  1:4605673/150‑1 (MQ=255)
                   cTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGc                                                                                                                                  <  2:812738/150‑1 (MQ=255)
                      aCCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAt                                                                                                                               <  1:4471022/150‑1 (MQ=255)
                      aCCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAt                                                                                                                               <  2:4429508/150‑1 (MQ=255)
                       ccAACCAGATAGTGAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAtt                                                                                                                              <  1:1952114/150‑1 (MQ=255)
                       ccAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGAtt                                                                                                                              >  2:666327/1‑150 (MQ=255)
                        cAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGAttt                                                                                                                             >  2:3328739/1‑150 (MQ=255)
                        cAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAtta                                                                                                                             <  1:918358/150‑2 (MQ=255)
                          aaCAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTg                                                                                                                           >  2:3957458/1‑150 (MQ=255)
                          aaCAGAGAGTTAAACGCCGTGAGTCCCGCGCCAATCACTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTATAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTg                                                                                                                           <  1:1362074/150‑1 (MQ=255)
                            cAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGGAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGAc                                                                                                                         <  2:2645404/150‑1 (MQ=255)
                             aGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACg                                                                                                                        >  1:703915/1‑150 (MQ=255)
                                tAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTa                                                                                                                     <  2:3504252/150‑1 (MQ=255)
                                tAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTa                                                                                                                     <  2:1964269/150‑1 (MQ=255)
                                 aGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTaa                                                                                                                    <  2:1770256/150‑1 (MQ=255)
                                  gTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAAc                                                                                                                   <  2:2122404/150‑1 (MQ=255)
                                   ttAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAAcc                                                                                                                  <  1:2519977/150‑1 (MQ=255)
                                    tAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACca                                                                                                                 >  2:1353234/1‑150 (MQ=255)
                                         gCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGa                                                                                                            <  1:1013444/150‑1 (MQ=255)
                                         gCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGa                                                                                                            >  1:84688/1‑150 (MQ=255)
                                         gCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGa                                                                                                            >  1:5021326/1‑150 (MQ=255)
                                          ccGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGAt                                                                                                           <  2:2455277/150‑1 (MQ=255)
                                           cGTGAGTCCCGCGCCAATCCATAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATa                                                                                                          <  2:4492887/150‑1 (MQ=255)
                                            gTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATan                                                                                                         <  1:4620469/150‑2 (MQ=255)
                                            gTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATaa                                                                                                         >  1:2929362/1‑150 (MQ=255)
                                             tGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATaaa                                                                                                        >  1:504014/1‑150 (MQ=255)
                                                 tCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAg                                                                                                    >  2:4734451/1‑150 (MQ=255)
                                                    cgcgCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCgg                                                                                                 >  2:5129545/1‑150 (MQ=255)
                                                      cgcCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGgca                                                                                               <  1:1173408/150‑1 (MQ=255)
                                                       gcCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGAGGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGgcag                                                                                              <  2:4363435/150‑1 (MQ=255)
                                                         cAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGgcagca                                                                                            >  1:376008/1‑150 (MQ=255)
                                                          aaTCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCaa                                                                                           >  1:3334213/1‑150 (MQ=255)
                                                           aTCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAAc                                                                                          >  2:4140805/1‑150 (MQ=255)
                                                            tCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACa                                                                                         <  1:4529/150‑1 (MQ=255)
                                                              ccTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAgg                                                                                       >  2:3416254/1‑150 (MQ=255)
                                                              ccTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAgg                                                                                       <  1:4178750/150‑1 (MQ=255)
                                                               cTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGt                                                                                      >  2:626538/1‑150 (MQ=255)
                                                                tAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTa                                                                                     <  1:3603473/150‑1 (MQ=255)
                                                                    gCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTa                                                                                 >  1:675300/1‑150 (MQ=255)
                                                                    gCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTa                                                                                 >  1:4941409/1‑150 (MQ=255)
                                                                     cAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTAc                                                                                <  1:1697385/150‑1 (MQ=255)
                                                                     cAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAGTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTAc                                                                                >  1:1784745/1‑150 (MQ=255)
                                                                      aaTCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACa                                                                               <  1:3632846/150‑1 (MQ=255)
                                                                        tcatcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGa                                                                             >  1:1571492/1‑150 (MQ=255)
                                                                        tcatcaATTTTTTAAAGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGTGGAAACTTTCTCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGGAGCTACAGa                                                                             >  1:4817371/1‑150 (MQ=255)
                                                                         catcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGAt                                                                            <  1:4038914/150‑1 (MQ=255)
                                                                           tcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTc                                                                          >  2:525062/1‑150 (MQ=255)
                                                                                    tAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAg                                                                 <  1:3903764/150‑1 (MQ=255)
                                                                                     aaCGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGt                                                                >  1:1475377/1‑150 (MQ=255)
                                                                                        gCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGa                                                             >  2:2064292/1‑150 (MQ=255)
                                                                                        gCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGa                                                             >  2:459237/1‑150 (MQ=255)
                                                                                            aTGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa                                                         <  1:3380221/150‑1 (MQ=255)
                                                                                            aTGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa                                                         <  1:415212/150‑1 (MQ=255)
                                                                                                gTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATActcc                                                     <  1:908650/150‑1 (MQ=255)
                                                                                                gTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATAAGCGGCAGCAACAGTTAGCTACAGATTCTGATAGAAGTTGACATActcc                                                     <  1:2022301/150‑1 (MQ=255)
                                                                                                        aaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                             >  1:1893578/1‑150 (MQ=255)
                                                                                                         aCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATc                                            >  1:3697616/1‑150 (MQ=255)
                                                                                                             cagTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt                                        >  1:2808475/1‑150 (MQ=255)
                                                                                                                 aaCAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTcc                                    <  2:2932207/150‑1 (MQ=255)
                                                                                                                 aaCAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTcc                                    >  1:3975768/1‑150 (MQ=255)
                                                                                                                     aGACAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTg                                >  2:583782/1‑150 (MQ=255)
                                                                                                                       aCAGGTAAAAAAGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTa                              <  2:1858736/150‑1 (MQ=255)
                                                                                                                       aCAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTa                              >  2:214665/1‑150 (MQ=255)
                                                                                                                          ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggc                           <  1:2715984/150‑1 (MQ=255)
                                                                                                                           gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCGGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                          >  2:2836902/1‑150 (MQ=255)
                                                                                                                           gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                          <  2:5033717/150‑1 (MQ=255)
                                                                                                                            tAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                         >  1:49047/1‑150 (MQ=255)
                                                                                                                             aaaaaaGTGACGATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcggc                        >  2:484883/1‑150 (MQ=255)
                                                                                                                             aaaaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcggc                        >  2:3652221/1‑150 (MQ=255)
                                                                                                                             aaaaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcggc                        >  1:3837062/1‑150 (MQ=255)
                                                                                                                                aaaGTGTCGATGGTGTTGGGGTAACTTTCGCGAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcg                     <  2:2509860/150‑1 (MQ=255)
                                                                                                                                     gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGc                >  2:3106063/1‑150 (MQ=255)
                                                                                                                                     gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGc                <  1:830122/150‑1 (MQ=255)
                                                                                                                                      aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCg               >  1:2362092/1‑150 (MQ=255)
                                                                                                                                       cGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGt              <  2:3463276/150‑1 (MQ=255)
                                                                                                                                        taTGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGtt             >  2:721538/2‑150 (MQ=255)
                                                                                                                                        gATGGTGATGGGGTAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGtt             >  2:3369121/1‑150 (MQ=255)
                                                                                                                                        gATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGtt             >  1:2367323/1‑150 (MQ=255)
                                                                                                                                          tGGTGATGTGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTgg           <  1:1416920/150‑1 (MQ=255)
                                                                                                                                            gTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCt         <  2:1069723/150‑1 (MQ=255)
                                                                                                                                            gTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCt         >  2:4074363/1‑150 (MQ=255)
                                                                                                                                             tGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTg        <  1:626562/150‑1 (MQ=255)
                                                                                                                                             tGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTg        <  2:3181263/150‑1 (MQ=255)
                                                                                                                                             tGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTg        <  2:4277074/150‑1 (MQ=255)
                                                                                                                                                 ggggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACt    >  1:4319476/1‑150 (MQ=255)
                                                                                                                                                  gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTc   <  2:2765014/150‑1 (MQ=255)
                                                                                                                                                  gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTc   >  2:2173182/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCA  >  CP009273/2124119‑2124415

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