Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,585,213 G→T 7.5% F152L (TTC→TTA yjjM ← yjjMN operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,585,2130GT 7.5% 1137.1 / 40.1 325F152L (TTC→TTAyjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (128/172);  new base T (10/15);  total (138/187)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 326 total aligned reads displayed.

TGTTGTCACGCAGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTG  >  CP009273/4585065‑4585362
                                                                                                                                                    |                                                                                                                                                     
  ttgtCACGCAGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAAt                                                                                                                                                    <  2:3088566/150‑1 (MQ=255)
         gCAGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCa                                                                                                                                             <  1:2473181/150‑1 (MQ=255)
           aGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATg                                                                                                                                           >  2:4680623/1‑150 (MQ=255)
           aGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATg                                                                                                                                           <  1:2611575/150‑1 (MQ=255)
             tAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCt                                                                                                                                         >  2:3535803/1‑150 (MQ=255)
             tAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCt                                                                                                                                         <  2:593661/150‑1 (MQ=255)
              aGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTc                                                                                                                                        >  2:3625800/1‑150 (MQ=255)
                  gATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACt                                                                                                                                    <  2:4246808/150‑1 (MQ=255)
                  gATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACt                                                                                                                                    >  1:3271960/1‑150 (MQ=255)
                  gATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACt                                                                                                                                    >  1:3224135/1‑150 (MQ=255)
                       cGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCcg                                                                                                                               >  1:3587226/1‑150 (MQ=255)
                              gCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCtttt                                                                                                                        <  1:686666/150‑1 (MQ=255)
                              gCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATATGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCtttt                                                                                                                        <  1:699883/150‑1 (MQ=255)
                               cATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTc                                                                                                                       >  1:1231423/1‑150 (MQ=255)
                               cATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTc                                                                                                                       <  1:674203/150‑1 (MQ=255)
                                     tttCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCtgtg                                                                                                                 >  2:4977392/1‑150 (MQ=255)
                                      ttCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGa                                                                                                                <  1:1854630/150‑1 (MQ=255)
                                      ttCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGa                                                                                                                <  1:1014906/150‑1 (MQ=255)
                                      ttCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGa                                                                                                                <  1:1985417/150‑1 (MQ=255)
                                        cGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAAt                                                                                                              <  1:1468174/150‑1 (MQ=255)
                                         gCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATa                                                                                                             <  2:5271268/150‑1 (MQ=255)
                                          ccTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATaa                                                                                                            <  2:2861726/150‑1 (MQ=255)
                                           cTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAg                                                                                                           <  2:2570793/150‑1 (MQ=255)
                                             gCAGCCAGCGGGGATCATGATCAGGCAGGTTGNGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGan                                                                                                         <  1:5057067/150‑2 (MQ=255)
                                              cagccagcGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGAtt                                                                                                        <  1:2863226/150‑1 (MQ=255)
                                               agccagcGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATta                                                                                                       >  2:3407886/1‑150 (MQ=255)
                                               agccagcGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATta                                                                                                       <  2:4728351/150‑1 (MQ=255)
                                                 ccagcGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATtata                                                                                                     <  2:228904/150‑1 (MQ=255)
                                                   agcGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAAc                                                                                                   <  1:2192946/150‑1 (MQ=255)
                                                      ggggATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGc                                                                                                <  1:1687175/150‑1 (MQ=255)
                                                      ggggATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGc                                                                                                <  1:1926380/150‑1 (MQ=255)
                                                      ggggATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGc                                                                                                >  2:712121/1‑150 (MQ=255)
                                                         gATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGa                                                                                             >  2:3234885/1‑150 (MQ=255)
                                                          aTCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGaa                                                                                            >  1:2302293/1‑150 (MQ=255)
                                                          aTCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGaa                                                                                            >  1:364349/1‑150 (MQ=255)
                                                           tCATGATCTGGCTGGTTGAGAAAGTGCTGATGTTAATGTGTTTCGAGCTTTTCGCGTAGCTCAAACAGCTGCTGGGCGTAGTTTTGATCGAATTGCTTCTTGCTCCACTGGCCGCGTTTTTCGCTGTTAATAATATTATAACGGCCGAAt                                                                                           <  1:2337155/150‑1 (MQ=255)
                                                            cATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAAtt                                                                                          >  1:1658306/1‑150 (MQ=255)
                                                             aTGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAAttt                                                                                         <  2:4091678/150‑1 (MQ=255)
                                                             aTGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAAttt                                                                                         >  2:4379220/1‑150 (MQ=255)
                                                              tGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAAtttt                                                                                        >  1:4178131/1‑150 (MQ=255)
                                                                aTCAGGCAGGTTGAGAAAGTGCTGATGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTaa                                                                                      <  1:4957998/150‑1 (MQ=255)
                                                                   aggcaggTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAg                                                                                   >  2:1758224/1‑150 (MQ=255)
                                                                    ggcaggTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGa                                                                                  <  1:3979205/150‑1 (MQ=255)
                                                                    ggcaggTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGATTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGa                                                                                  <  1:75421/150‑1 (MQ=255)
                                                                       aggTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGatat                                                                               <  2:5294238/150‑1 (MQ=255)
                                                                       aggTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGatat                                                                               <  2:4600887/150‑1 (MQ=255)
                                                                          ttGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATtct                                                                            <  1:1746020/150‑1 (MQ=255)
                                                                          ttGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACGGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATtct                                                                            <  2:2009070/150‑1 (MQ=255)
                                                                           tGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATtatc                                                                           >  1:3805898/1‑150 (MQ=255)
                                                                            gagaAAGTGCTGATGTGAATGTGTTTCGAGCTTTTGGGGTAGCTCAAACAGCTGCTCGGCGTTTTATTGATCGATTTGATTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATATGATTATAACGGCCGAATTTTAAAAGATATTctct                                                                          <  1:2741552/150‑1 (MQ=255)
                                                                            gagaAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTctct                                                                          >  2:1631434/1‑150 (MQ=255)
                                                                            gagaAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTctct                                                                          <  2:5167640/150‑1 (MQ=255)
                                                                             agaAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCtt                                                                         >  1:4018271/1‑150 (MQ=255)
                                                                              gaAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTta                                                                        <  1:1499765/150‑1 (MQ=255)
                                                                               aaaGTGCTGATGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATATGATTATAACGGCCGAATTTTAAAAGATATTCTCTtac                                                                       <  1:5082966/150‑1 (MQ=255)
                                                                               aaaGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTtac                                                                       <  1:4291181/150‑1 (MQ=255)
                                                                                  gTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTtactac                                                                    >  2:2394954/1‑150 (MQ=255)
                                                                                     cTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAAc                                                                 >  1:2063321/1‑150 (MQ=255)
                                                                                     cTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGACGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAAc                                                                 >  2:2821539/1‑150 (MQ=255)
                                                                                      ttaaGTTAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAAcc                                                                <  1:552262/148‑1 (MQ=255)
                                                                                      tGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAAcc                                                                >  2:1094943/1‑150 (MQ=255)
                                                                                         aGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCggg                                                             >  1:3303411/1‑150 (MQ=255)
                                                                                           tGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATACCGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCt                                                           >  2:224727/1‑150 (MQ=255)
                                                                                           tGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCgggat                                                           <  1:4198793/150‑3 (MQ=255)
                                                                                           tGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCt                                                           >  2:263160/1‑150 (MQ=255)
                                                                                             aaTGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGa                                                         <  2:139299/150‑1 (MQ=255)
                                                                                              aTGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGac                                                        <  2:379976/150‑1 (MQ=255)
                                                                                                 tgtTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACAcc                                                     <  1:3016251/150‑1 (MQ=255)
                                                                                                  gtTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACAccc                                                    >  1:3570273/1‑150 (MQ=255)
                                                                                                  gtTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACAccc                                                    <  2:796059/150‑1 (MQ=255)
                                                                                                  gtTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACAccc                                                    >  2:5083836/1‑150 (MQ=255)
                                                                                                        aGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGcc                                              <  2:4873862/150‑1 (MQ=255)
                                                                                                         gCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGccc                                             <  2:1231423/150‑1 (MQ=255)
                                                                                                          cATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCg                                            >  1:1084647/1‑150 (MQ=255)
                                                                                                          cATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCg                                            >  2:2941265/1‑150 (MQ=255)
                                                                                                            tttCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAg                                          >  1:5215087/1‑150 (MQ=255)
                                                                                                               cgcgTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCa                                       >  1:4768736/1‑150 (MQ=255)
                                                                                                                gcgTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAg                                      >  1:1304092/1‑150 (MQ=255)
                                                                                                                 cgTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGt                                     >  2:4512915/1‑150 (MQ=255)
                                                                                                                  gTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGTCCGCGATTTTCGCTGTGAATAATATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGtt                                    <  1:4489832/150‑1 (MQ=255)
                                                                                                                  gTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGtt                                    <  2:3021631/150‑1 (MQ=255)
                                                                                                                  gTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGtt                                    >  2:2118655/1‑150 (MQ=255)
                                                                                                                    aGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGc                                  >  1:894310/1‑150 (MQ=255)
                                                                                                                         aaaCAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTc                             >  2:891048/1‑150 (MQ=255)
                                                                                                                               ctgctCGGCGTAGGATTGATCTAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGtaaa                       <  1:4794145/150‑1 (MQ=255)
                                                                                                                                 gctCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACg                     <  1:4673964/150‑1 (MQ=255)
                                                                                                                                 gctCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACg                     <  1:1796285/150‑1 (MQ=255)
                                                                                                                                  ctCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGt                    <  2:2510411/150‑1 (MQ=255)
                                                                                                                                   tCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGtt                   <  2:4178106/150‑1 (MQ=255)
                                                                                                                                    cGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATCACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGttt                  <  1:1516227/150‑1 (MQ=255)
                                                                                                                                    cGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGttt                  <  2:1169861/150‑1 (MQ=255)
                                                                                                                                     ggCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTc                 >  1:4454265/1‑150 (MQ=255)
                                                                                                                                      gCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGTATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTcc                <  2:4041661/150‑1 (MQ=255)
                                                                                                                                      gCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTcc                >  1:2694566/1‑150 (MQ=255)
                                                                                                                                       cGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTccc               >  1:618577/1‑150 (MQ=255)
                                                                                                                                         tAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTCCCCt             <  2:4795340/150‑1 (MQ=255)
                                                                                                                                         tAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTCCCCt             <  1:3480248/150‑1 (MQ=255)
                                                                                                                                               tGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTCCCCTGGGCGt       <  1:1266408/150‑1 (MQ=255)
                                                                                                                                                  tCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTCCCCTGGGCGTAAc    <  1:4977392/150‑1 (MQ=255)
                                                                                                                                                   cGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGTAAACGTTTCCCCTGGGCGTAACt   <  2:188289/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTTGTCACGCAGTAGCCGATGACGTTCGAGCATGGTTTTCGCCTGCAGCCAGCGGGGATCATGATCAGGCAGGTTGAGAAAGTGCTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTG  >  CP009273/4585065‑4585362

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: