Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,065,158 T→A 12.1% S47T (TCA→ACA)  flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,065,1580TA12.1% 747.0 / 73.9 246S47T (TCA→ACA) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (93/123);  new base A (11/19);  total (104/142)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

GGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATAACGAG  >  CP009273/2065010‑2065307
                                                                                                                                                    |                                                                                                                                                     
  aaaaGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAgtcacgtcac                                                                                                                                                    >  2:269148/1‑150 (MQ=255)
    aaGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTc                                                                                                                                                  >  1:1714654/1‑150 (MQ=255)
     aGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTcc                                                                                                                                                 <  1:3942105/150‑1 (MQ=255)
     aGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTcc                                                                                                                                                 <  1:1681630/150‑1 (MQ=255)
       cTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCg                                                                                                                                               <  2:2744997/150‑1 (MQ=255)
        tgatgaAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCgg                                                                                                                                              >  1:689150/1‑150 (MQ=255)
         gatgaAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGt                                                                                                                                             <  1:1211650/150‑1 (MQ=255)
            gaAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgct                                                                                                                                          >  2:2972957/1‑150 (MQ=255)
            gaAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgct                                                                                                                                          >  2:691908/1‑150 (MQ=255)
             aaaCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctg                                                                                                                                         <  2:1141115/150‑1 (MQ=255)
             aaaCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctg                                                                                                                                         <  2:3719213/150‑1 (MQ=255)
              aaCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctgg                                                                                                                                        <  2:1678030/150‑1 (MQ=255)
              aaCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTgctgg                                                                                                                                        <  1:662123/150‑1 (MQ=255)
               aCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctggc                                                                                                                                       >  2:4423099/1‑150 (MQ=255)
               aCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctggc                                                                                                                                       >  2:4425162/1‑150 (MQ=255)
                 gACATCTGAAGACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTgctggctg                                                                                                                                     <  2:311833/150‑1 (MQ=255)
                   cATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGgctgct                                                                                                                                   <  1:2835319/150‑1 (MQ=255)
                    aTCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGctgctg                                                                                                                                  >  1:1246958/1‑150 (MQ=255)
                     tCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGa                                                                                                                                 >  2:330840/1‑150 (MQ=255)
                       tGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACa                                                                                                                               >  1:89462/1‑150 (MQ=255)
                         aaTACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATc                                                                                                                             <  1:4225403/150‑1 (MQ=255)
                               tGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGTCATCGTtgtg                                                                                                                       <  1:3234124/150‑1 (MQ=255)
                                gCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGc                                                                                                                      >  2:3039211/1‑150 (MQ=255)
                                 cTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCa                                                                                                                     >  2:4502488/1‑150 (MQ=255)
                                     aGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCg                                                                                                                 >  1:588986/1‑150 (MQ=255)
                                     aGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCg                                                                                                                 >  2:2742681/1‑150 (MQ=255)
                                      ggCTGGTATGAAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCgg                                                                                                                <  2:1225138/150‑1 (MQ=255)
                                        cTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCggca                                                                                                              <  2:3717327/150‑3 (MQ=255)
                                          ggTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGgaga                                                                                                            >  2:4072081/1‑150 (MQ=255)
                                          ggTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGgaga                                                                                                            >  1:137337/1‑150 (MQ=255)
                                           gTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGAGCACGGGGTAAACGTGGCGGTGTGGCGGTAGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGaa                                                                                                           >  1:2931124/1‑150 (MQ=255)
                                             aTGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCTGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAAc                                                                                                         >  1:4312724/1‑150 (MQ=255)
                                             aTGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAAc                                                                                                         <  2:1960104/150‑1 (MQ=255)
                                              tGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAAcc                                                                                                        <  1:3070330/150‑1 (MQ=255)
                                              tGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAAcc                                                                                                        <  1:3070316/150‑1 (MQ=255)
                                              tGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAAcc                                                                                                        <  1:3418073/150‑1 (MQ=255)
                                                 aaTCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTg                                                                                                     >  1:4376039/1‑150 (MQ=255)
                                                    cacaTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAc                                                                                                  <  2:2550786/150‑1 (MQ=255)
                                                    cacaTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAc                                                                                                  <  2:704150/150‑1 (MQ=255)
                                                     acaTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcg                                                                                                 <  2:3950163/150‑1 (MQ=255)
                                                     acaTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcg                                                                                                 <  2:922125/150‑1 (MQ=255)
                                                      caTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcgg                                                                                                >  1:887873/1‑150 (MQ=255)
                                                      caTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcgg                                                                                                >  1:1412309/1‑150 (MQ=255)
                                                       aTGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAGCGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcggc                                                                                               >  2:2122795/1‑150 (MQ=255)
                                                         gACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAAcggcgg                                                                                             >  1:4050389/1‑150 (MQ=255)
                                                            gggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAac                                                                                          <  1:790555/150‑1 (MQ=255)
                                                            gggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAac                                                                                          <  1:913479/150‑1 (MQ=255)
                                                            gggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAac                                                                                          <  1:3989908/150‑1 (MQ=255)
                                                            gggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAac                                                                                          <  1:3637461/150‑1 (MQ=255)
                                                             ggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAaca                                                                                         <  1:916990/150‑1 (MQ=255)
                                                             ggCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGATGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAaca                                                                                         <  1:2524616/150‑1 (MQ=255)
                                                               cgcTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACt                                                                                       <  1:542616/150‑1 (MQ=255)
                                                               cgcTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACt                                                                                       <  1:264670/150‑1 (MQ=255)
                                                               cgcTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACt                                                                                       <  1:259627/150‑1 (MQ=255)
                                                                 cTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTgg                                                                                     >  1:2105546/1‑150 (MQ=255)
                                                                  tttCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGc                                                                                    <  2:3334390/150‑1 (MQ=255)
                                                                      gTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAAt                                                                                <  1:3338751/150‑1 (MQ=255)
                                                                           tGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGa                                                                           >  2:22062/1‑150 (MQ=255)
                                                                               tCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAAc                                                                       <  1:3628248/150‑1 (MQ=255)
                                                                               tCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAAc                                                                       <  1:793696/150‑1 (MQ=255)
                                                                               tCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAAc                                                                       <  1:623540/150‑1 (MQ=255)
                                                                               tCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAAc                                                                       <  1:3076080/150‑1 (MQ=255)
                                                                                 cGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCa                                                                     >  1:1523160/1‑150 (MQ=255)
                                                                                  gAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAg                                                                    <  1:14613/150‑1 (MQ=255)
                                                                                     cTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGAtt                                                                 <  1:956146/150‑1 (MQ=255)
                                                                                         cccGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCt                                                             <  1:3517224/150‑1 (MQ=255)
                                                                                         cccGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCt                                                             >  1:1036379/1‑150 (MQ=255)
                                                                                             cgcACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCgg                                                         >  2:3768307/1‑150 (MQ=255)
                                                                                              gcACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGt                                                        >  2:2799454/1‑150 (MQ=255)
                                                                                               cACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTa                                                       >  2:3623460/1‑150 (MQ=255)
                                                                                               cACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTa                                                       <  2:3236942/150‑1 (MQ=255)
                                                                                               cACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTa                                                       <  1:205035/150‑1 (MQ=255)
                                                                                                aCGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTAc                                                      <  2:603130/150‑1 (MQ=255)
                                                                                                 cGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACg                                                     <  1:1834052/150‑1 (MQ=255)
                                                                                                      tAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCaa                                                <  2:20633/150‑1 (MQ=255)
                                                                                                       aaaCGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAAc                                               >  1:379037/1‑150 (MQ=255)
                                                                                                       aaaCGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAAc                                               <  2:683887/150‑1 (MQ=255)
                                                                                                        aaCGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGGTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACg                                              <  2:1531350/150‑1 (MQ=255)
                                                                                                         aCGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAATgg                                             <  1:1059966/150‑1 (MQ=255)
                                                                                                           gtggcggtgtggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGaa                                           >  2:3193352/1‑150 (MQ=255)
                                                                                                             ggcggtgtggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACATCCAGATTGTCTTCGGTACGACCAACGGAATg                                         <  2:3529742/150‑1 (MQ=255)
                                                                                                             ggcggtgtggcggtTGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATg                                         >  2:2640503/1‑150 (MQ=255)
                                                                                                                ggtgtggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGAcc                                      <  1:2496074/150‑1 (MQ=255)
                                                                                                                  tgtggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCAt                                    <  2:3399200/150‑1 (MQ=255)
                                                                                                                    tggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCa                                  <  1:526993/150‑1 (MQ=255)
                                                                                                                    tggcggtTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCa                                  <  2:1862411/150‑1 (MQ=255)
                                                                                                                          ttGCACTGTCTCTTGCCGCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCg                            >  2:2226757/1‑150 (MQ=255)
                                                                                                                           tGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCgg                           >  2:2405716/1‑150 (MQ=255)
                                                                                                                            gCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCggg                          <  1:1658205/150‑1 (MQ=255)
                                                                                                                            gCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCggg                          <  1:292713/150‑1 (MQ=255)
                                                                                                                                 gTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGa                     >  2:1754727/1‑150 (MQ=255)
                                                                                                                                    tctTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTa                  <  1:2158517/150‑1 (MQ=255)
                                                                                                                                           gCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCg           <  2:4172529/150‑1 (MQ=255)
                                                                                                                                           gCAGTCACGACACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCg           >  1:3695151/1‑150 (MQ=255)
                                                                                                                                            cAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCgg          <  1:2791826/150‑1 (MQ=255)
                                                                                                                                               tcacgtcacTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATa       <  1:4535806/150‑1 (MQ=255)
                                                                                                                                                cacgtcacTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATaa      >  1:2969144/1‑150 (MQ=255)
                                                                                                                                                  cgtcacTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATAACg    <  1:330840/150‑1 (MQ=255)
                                                                                                                                                   gtcacTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATAACGa   >  1:4077364/1‑150 (MQ=255)
                                                                                                                                                    tcacTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATAACGAg  <  1:2684011/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTGTGGCGGTTGCACTGTCTCTTGCCGCAGTCACGTCACTCCCGGTGCTGGCTGCTGACATCGTTGTGCACCCGGGAGAAACCGTGAACGGCGGAACACTGGCAAATCATGACAACCAGATTGTCTTCGGTACGACCAACGGAATGACCATCAGTACCGGGCTGGAGTATGGGCCGGATAACGAG  >  CP009273/2065010‑2065307

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: