Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,088,850 C→A 6.6% R93S (CGT→AGT)  rhaR → transcriptional activator of rhaSR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,088,8500CA 6.6% 1008.9 / 29.1 264R93S (CGT→AGT) rhaRtranscriptional activator of rhaSR
Reads supporting (aligned to +/- strand):  ref base C (141/105);  new base A (6/12);  total (147/117)
Fisher's exact test for biased strand distribution p-value = 5.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

TTGTCTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATGTG  >  CP009273/4088701‑4088996
                                                                                                                                                     |                                                                                                                                                  
ttGTCTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGc                                                                                                                                                    >  1:276195/1‑150 (MQ=255)
 tGTCTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCg                                                                                                                                                   >  1:3027903/1‑150 (MQ=255)
               aTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATc                                                                                                                                     <  2:2531005/150‑1 (MQ=255)
                        aTGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGc                                                                                                                            >  1:3104633/1‑150 (MQ=255)
                         tGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCa                                                                                                                           <  2:1034271/150‑1 (MQ=255)
                           tACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAgg                                                                                                                         >  2:2033202/1‑150 (MQ=255)
                             cTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAgggg                                                                                                                       <  1:612515/150‑1 (MQ=255)
                               cAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGc                                                                                                                     <  2:1201187/150‑1 (MQ=255)
                                aaCGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                                    >  2:4499905/1‑150 (MQ=255)
                                   gATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAtt                                                                                                                 <  1:2543816/150‑1 (MQ=255)
                                    aTCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTc                                                                                                                >  1:4238313/1‑150 (MQ=255)
                                    aTCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTc                                                                                                                >  1:737662/1‑150 (MQ=255)
                                     tCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTcc                                                                                                               <  2:1976268/150‑1 (MQ=255)
                                      cGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCg                                                                                                              >  1:3166781/1‑150 (MQ=255)
                                       gCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCgg                                                                                                             >  1:1975057/1‑150 (MQ=255)
                                        cccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCggg                                                                                                            >  2:125062/1‑150 (MQ=255)
                                         ccTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGa                                                                                                           >  2:4340193/1‑150 (MQ=255)
                                             aTCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                       >  2:1773410/1‑150 (MQ=255)
                                             aTCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTa                                                                                                       >  1:3422531/1‑150 (MQ=255)
                                                 cATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGc                                                                                                   >  2:2791110/1‑150 (MQ=255)
                                                    tACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACCCTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAg                                                                                                <  2:2439267/150‑1 (MQ=255)
                                                    tACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAg                                                                                                >  1:635272/1‑150 (MQ=255)
                                                     aCCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAgc                                                                                               >  1:2811598/1‑150 (MQ=255)
                                                        cGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCa                                                                                            >  1:1888921/1‑150 (MQ=255)
                                                          tGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAgg                                                                                          <  2:542227/150‑1 (MQ=255)
                                                           ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                         <  1:3530521/150‑1 (MQ=255)
                                                             ctaTCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                       <  1:3810364/148‑1 (MQ=255)
                                                             cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                       >  1:1859844/1‑150 (MQ=255)
                                                             cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                       >  2:2263907/1‑150 (MQ=255)
                                                             cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                       >  2:2019276/1‑150 (MQ=255)
                                                              gATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCaa                                                                                      >  1:3256057/1‑150 (MQ=255)
                                                              gATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAAAATTATTTATTGCCCGGAGAGTCTGATTCTGAATCATGACTGGCAGTGGGCGAATCCGGGATGTAACGCCCGCGCAGGGcca                                                                                      >  2:3202462/1‑148 (MQ=255)
                                                               aTCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAc                                                                                     <  2:3775266/150‑1 (MQ=255)
                                                                tctcTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAcc                                                                                    <  1:2001833/150‑1 (MQ=255)
                                                                      ttACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTg                                                                              <  2:826284/150‑1 (MQ=255)
                                                                        aCATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgc                                                                            >  2:1351782/1‑150 (MQ=255)
                                                                         cATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                           <  1:431043/150‑1 (MQ=255)
                                                                         cATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                           <  1:4333/150‑1 (MQ=255)
                                                                         cATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                           >  2:4528882/1‑150 (MQ=255)
                                                                              aTGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAg                                                                      <  2:4102670/150‑1 (MQ=255)
                                                                                gCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGt                                                                    >  1:29280/1‑150 (MQ=255)
                                                                                gCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGt                                                                    >  1:2350710/1‑150 (MQ=255)
                                                                                 cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTa                                                                   >  1:469905/1‑150 (MQ=255)
                                                                                 cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTa                                                                   <  2:1791495/150‑1 (MQ=255)
                                                                                  tGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCgtaggtag                                                                  >  2:1346682/1‑150 (MQ=255)
                                                                                  tGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAg                                                                  <  1:166474/150‑1 (MQ=255)
                                                                                     cGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCAt                                                               <  2:2949700/150‑1 (MQ=255)
                                                                                      gATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATg                                                              >  2:4040194/1‑150 (MQ=255)
                                                                                      gATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCACGGAGCGTCTGAAGCTGAATCTTGACTGGCAGTGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCACCCACACTGGCGCTTAGGTAGCATg                                                              >  2:2803884/1‑150 (MQ=255)
                                                                                      gATAAACACTCATACGCTTCCGTTAACGATCTGGTTTAGCAGAATATTATTTATTGCCCGGAGCGACTGAAGCTGAATCTTGACTGTCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACAGGCGCTTAGGTAGCaag                                                              >  2:3614880/1‑148 (MQ=255)
                                                                                            cacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATg                                                        >  1:1480424/1‑150 (MQ=255)
                                                                                               tCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGgcg                                                     <  1:2424908/150‑1 (MQ=255)
                                                                                               tCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGgcg                                                     >  1:499995/1‑150 (MQ=255)
                                                                                                  tACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAg                                                  >  2:2871965/1‑150 (MQ=255)
                                                                                                   aCGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAgg                                                 >  1:4229629/1‑150 (MQ=255)
                                                                                                   aCGCGTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAgg                                                 <  2:2709350/150‑1 (MQ=255)
                                                                                                     gCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGgcg                                               >  2:793427/1‑150 (MQ=255)
                                                                                                          cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCa                                          >  2:1905207/1‑150 (MQ=255)
                                                                                                           gTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAg                                         >  2:1606678/1‑150 (MQ=255)
                                                                                                            ttAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAgg                                        <  2:2152853/150‑1 (MQ=255)
                                                                                                             tAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGt                                       >  2:3158512/1‑150 (MQ=255)
                                                                                                               aCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTa                                     >  1:2613877/1‑150 (MQ=255)
                                                                                                                cGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTAt                                    <  1:908188/150‑1 (MQ=255)
                                                                                                                 gATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                                   <  1:420198/150‑1 (MQ=255)
                                                                                                                 gATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                                   <  2:689856/150‑1 (MQ=255)
                                                                                                                  aTCTGGTTTTGCAGAATATTCTTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCg                                  <  1:3356623/150‑1 (MQ=255)
                                                                                                                   tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCgg                                 <  2:92003/150‑1 (MQ=255)
                                                                                                                      ggTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCa                              >  2:782590/1‑150 (MQ=255)
                                                                                                                      ggTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCa                              >  1:179066/1‑150 (MQ=255)
                                                                                                                      ggTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCa                              <  1:3836819/150‑1 (MQ=255)
                                                                                                                        ttttGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGc                            >  1:1046139/1‑150 (MQ=255)
                                                                                                                        ttttGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGc                            >  1:3674151/1‑150 (MQ=255)
                                                                                                                         tttGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCt                           >  2:3368760/1‑150 (MQ=255)
                                                                                                                          ttGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCtt                          <  2:1306349/150‑1 (MQ=255)
                                                                                                                           tGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                         <  1:4479145/150‑1 (MQ=255)
                                                                                                                           tGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                         <  1:3070680/150‑1 (MQ=255)
                                                                                                                            gCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGa                        >  2:3310921/1‑150 (MQ=255)
                                                                                                                            gCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGa                        <  2:482749/150‑1 (MQ=255)
                                                                                                                             cAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAg                       >  1:3736779/1‑150 (MQ=255)
                                                                                                                               gAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCa                     <  2:4215003/150‑1 (MQ=255)
                                                                                                                               gAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCa                     >  1:1381970/1‑150 (MQ=255)
                                                                                                                               gAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCa                     <  2:4465994/150‑1 (MQ=255)
                                                                                                                                aaTATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCAt                    >  1:3875524/1‑150 (MQ=255)
                                                                                                                                  tattatTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGa                  >  2:2065853/1‑150 (MQ=255)
                                                                                                                                     tatttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAag               >  1:3246189/1‑150 (MQ=255)
                                                                                                                                     tatttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAag               >  2:2611005/1‑150 (MQ=255)
                                                                                                                                      atttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGTATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagt              >  2:1111043/1‑150 (MQ=255)
                                                                                                                                        ttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtag            >  1:224243/1‑150 (MQ=255)
                                                                                                                                        ttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtag            >  2:2414258/1‑150 (MQ=255)
                                                                                                                                        ttattGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtag            >  2:4144700/1‑150 (MQ=255)
                                                                                                                                          attGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTc          >  2:310739/1‑150 (MQ=255)
                                                                                                                                          attGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTc          >  1:1125029/1‑150 (MQ=255)
                                                                                                                                            tGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAg        <  2:3680743/150‑1 (MQ=255)
                                                                                                                                              cccGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCa      >  2:641435/1‑150 (MQ=255)
                                                                                                                                               ccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAt     <  1:4289969/150‑1 (MQ=255)
                                                                                                                                                cGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtg    <  2:1939290/150‑1 (MQ=255)
                                                                                                                                                 ggAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtgt   <  1:2841291/150‑1 (MQ=255)
                                                                                                                                                 ggAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtgt   >  2:1530570/1‑150 (MQ=255)
                                                                                                                                                 ggAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtgt   >  2:3653449/1‑150 (MQ=255)
                                                                                                                                                  gAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtgtg  <  1:784595/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTGTCTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATGTG  >  CP009273/4088701‑4088996

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: