Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,040,800 C→A 14.1% T794T (ACC→ACA yeeJ → putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,040,8000CA14.1% 663.3 / 63.4 207T794T (ACC→ACAyeeJputative adhesin
Reads supporting (aligned to +/- strand):  ref base C (85/92);  new base A (12/18);  total (97/110)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

GCCAGCAATAATGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGATCT  >  CP009273/2040651‑2040949
                                                                                                                                                     |                                                                                                                                                     
gCCAGCAATAATGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACAcc                                                                                                                                                       <  2:2605528/150‑1 (MQ=255)
 ccAGCAATAATGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCg                                                                                                                                                      <  2:1963807/150‑1 (MQ=255)
  caggtataatGGTGTGCTCGCCATTGAGAATGCAGCAAACACCGTCTCGTTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGTATCTGCAACTTCCTTCAACTATCAAAACACCGc                                                                                                                                                     <  2:2501308/145‑1 (MQ=255)
    gCAATAATGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCaa                                                                                                                                                   <  2:422421/150‑1 (MQ=255)
        taatGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAAc                                                                                                                                               >  1:2031637/1‑150 (MQ=255)
               gtgCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGtt                                                                                                                                        <  2:2707371/150‑1 (MQ=255)
                tgCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTa                                                                                                                                       <  2:1341416/150‑1 (MQ=255)
                tgCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTa                                                                                                                                       >  1:2046257/1‑150 (MQ=255)
                tgCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTa                                                                                                                                       >  2:3383764/1‑150 (MQ=255)
                    cGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATgg                                                                                                                                   <  1:2810826/150‑1 (MQ=255)
                     gCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGt                                                                                                                                  <  1:2693661/150‑1 (MQ=255)
                     gCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGt                                                                                                                                  <  1:2619404/150‑1 (MQ=255)
                      ccAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTc                                                                                                                                 <  1:2286168/150‑1 (MQ=255)
                          tGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGc                                                                                                                             <  2:3641131/150‑1 (MQ=255)
                          tGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGc                                                                                                                             <  1:1320974/150‑1 (MQ=255)
                            agaATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGa                                                                                                                           >  1:4216839/1‑150 (MQ=255)
                              aaTGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACt                                                                                                                         <  1:32116/150‑1 (MQ=255)
                              aaTGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACt                                                                                                                         >  1:3919814/1‑150 (MQ=255)
                               aTGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACtt                                                                                                                        >  1:786845/1‑150 (MQ=255)
                               aTGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAAACAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACtt                                                                                                                        <  2:3750466/150‑1 (MQ=255)
                                   agcaAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGa                                                                                                                    >  1:4040374/1‑150 (MQ=255)
                                     caAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATc                                                                                                                  <  1:4173206/150‑1 (MQ=255)
                                      aaaCACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCCATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCt                                                                                                                 <  2:1365378/150‑1 (MQ=255)
                                       aaCACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCAGACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTg                                                                                                                <  2:587930/150‑1 (MQ=255)
                                        acacCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGa                                                                                                               >  1:2942039/1‑150 (MQ=255)
                                         cacCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGaa                                                                                                              >  2:2618571/1‑150 (MQ=255)
                                         cacCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGaa                                                                                                              >  1:776196/1‑150 (MQ=255)
                                          acCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGaaa                                                                                                             >  2:3762864/1‑150 (MQ=255)
                                                  ggTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTaa                                                                                                     >  1:3033844/1‑150 (MQ=255)
                                                  ggTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTaa                                                                                                     >  2:149922/1‑150 (MQ=255)
                                                       aTGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAgg                                                                                                <  1:609234/150‑1 (MQ=255)
                                                         gTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGaa                                                                                              <  1:1960135/150‑1 (MQ=255)
                                                             cTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGaca                                                                                          <  1:1558748/150‑1 (MQ=255)
                                                             cTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGaca                                                                                          >  1:1264745/1‑150 (MQ=255)
                                                             cTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGaca                                                                                          <  1:648408/150‑1 (MQ=255)
                                                                   aaggaagCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACgg                                                                                    >  1:3163186/1‑150 (MQ=255)
                                                                    aggaagCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGt                                                                                   >  1:1239531/1‑150 (MQ=255)
                                                                       aagCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGa                                                                                <  2:2879637/150‑1 (MQ=255)
                                                                       aagCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGa                                                                                >  1:3617477/1‑150 (MQ=255)
                                                                        agCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGaa                                                                               >  2:2565240/1‑150 (MQ=255)
                                                                              ccAATCAATGATCAGACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCAcc                                                                         <  1:2707638/150‑1 (MQ=255)
                                                                                aatcaatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCt                                                                       >  1:347505/1‑150 (MQ=255)
                                                                                aatcaatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCt                                                                       >  1:1494725/1‑150 (MQ=255)
                                                                                 atcaatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCtt                                                                      >  1:4030098/1‑150 (MQ=255)
                                                                                 atcaatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCtt                                                                      <  2:344672/150‑1 (MQ=255)
                                                                                   caatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGACGACAACACGGTTGAAGTCACCCTTGa                                                                    >  2:3067031/1‑150 (MQ=255)
                                                                                   caatGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGa                                                                    <  1:1502123/150‑1 (MQ=255)
                                                                                     atGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGaaa                                                                  >  1:1167298/1‑150 (MQ=255)
                                                                                        aTCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATg                                                               <  2:760645/150‑1 (MQ=255)
                                                                                        aTCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATg                                                               >  2:1540666/1‑150 (MQ=255)
                                                                                          cATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGc                                                             >  1:2646690/1‑150 (MQ=255)
                                                                                           agaCCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCg                                                            <  1:2889623/148‑1 (MQ=255)
                                                                                            tACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGt                                                           <  2:1571533/150‑1 (MQ=255)
                                                                                                 tCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGaaac                                                      <  1:4226889/150‑1 (MQ=255)
                                                                                                 tCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGaaac                                                      <  1:2206652/150‑1 (MQ=255)
                                                                                                   aCGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGaaacaa                                                    <  2:2046257/150‑1 (MQ=255)
                                                                                                     gTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGaaacaaac                                                  <  2:817803/150‑1 (MQ=255)
                                                                                                          cGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTaa                                             <  2:2370/150‑1 (MQ=255)
                                                                                                           ggTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAAt                                            <  1:864240/150‑1 (MQ=255)
                                                                                                            gTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATc                                           <  1:3383764/150‑1 (MQ=255)
                                                                                                            gTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATc                                           <  2:1234806/150‑1 (MQ=255)
                                                                                                               ttAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTc                                        <  1:384629/150‑1 (MQ=255)
                                                                                                                tAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCa                                       >  2:4212741/1‑150 (MQ=255)
                                                                                                                tAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCa                                       >  2:4213167/1‑150 (MQ=255)
                                                                                                                tAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCa                                       >  1:3396561/1‑150 (MQ=255)
                                                                                                                  aGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCAGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGGTGAAGTCACCCTTGAAAATGGCGTGAAACACACGTTACTCGTCAGt                                     >  2:142487/1‑150 (MQ=255)
                                                                                                                     ggATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGtttt                                  >  2:1916333/1‑150 (MQ=255)
                                                                                                                     ggATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGtttt                                  >  2:4334485/1‑150 (MQ=255)
                                                                                                                      gATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTg                                 >  1:2759051/1‑150 (MQ=255)
                                                                                                                       aTCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGt                                <  2:1596569/150‑1 (MQ=255)
                                                                                                                         cGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCg                              >  1:1021088/1‑150 (MQ=255)
                                                                                                                          ggCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCgg                             <  2:505730/150‑1 (MQ=255)
                                                                                                                          ggCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCgg                             >  2:1740494/1‑150 (MQ=255)
                                                                                                                          ggCAACTTCCTACAACAATCAAAACACCGCAACAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCgg                             <  1:2368503/150‑1 (MQ=255)
                                                                                                                            cAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCg                           >  2:2102459/1‑150 (MQ=255)
                                                                                                                             aaCTTCCTTCAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGa                          <  2:1049917/150‑1 (MQ=255)
                                                                                                                              aCTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGAc                         <  2:1101493/150‑1 (MQ=255)
                                                                                                                              aCTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGAc                         >  1:3520871/1‑150 (MQ=255)
                                                                                                                               cttccttcAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACt                        >  1:1115787/1‑150 (MQ=255)
                                                                                                                               cttccttcAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACt                        <  2:2512963/150‑1 (MQ=255)
                                                                                                                                   cttcAACAATCAAAACACAGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAg                    >  2:4060470/1‑150 (MQ=255)
                                                                                                                                    ttcAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGt                   <  1:1606836/150‑1 (MQ=255)
                                                                                                                                     acaacaaTCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTa                  <  2:3399393/149‑1 (MQ=255)
                                                                                                                                        acaaTCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTg               >  1:3813801/1‑150 (MQ=255)
                                                                                                                                        acaaTCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTg               >  2:2071959/1‑150 (MQ=255)
                                                                                                                                           aTCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTCCTgcgc            >  1:2982616/1‑150 (MQ=255)
                                                                                                                                           aTCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTgcgc            >  1:372726/1‑150 (MQ=255)
                                                                                                                                            tCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCa           <  1:3354484/150‑1 (MQ=255)
                                                                                                                                              aaaaCACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAgg         >  1:2044891/1‑150 (MQ=255)
                                                                                                                                               aaaCACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGt        <  1:50859/150‑1 (MQ=255)
                                                                                                                                                aaCACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGtt       <  1:4204224/150‑1 (MQ=255)
                                                                                                                                                 acacCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTg      <  1:663668/150‑1 (MQ=255)
                                                                                                                                                  cacCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGa     <  1:4101374/150‑1 (MQ=255)
                                                                                                                                                  cacCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGa     >  2:1425176/1‑150 (MQ=255)
                                                                                                                                                  cacCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGa     <  1:1911452/150‑1 (MQ=255)
                                                                                                                                                   acCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGAt    >  1:595143/1‑150 (MQ=255)
                                                                                                                                                   acCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGAt    <  1:1759600/150‑1 (MQ=255)
                                                                                                                                                    ccGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGATc   <  2:3975630/150‑1 (MQ=255)
                                                                                                                                                    ccGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGATc   >  1:2947174/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGATCt  <  2:2276163/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGCAATAATGGTGTGCTCGCCAATGAGAATGCAGCAAACACCGTCTCGGTCAATGTCGCTGATGAAGGAAGCAACCCAATCAATGATCATACCGTCACGTTTGCGGTATTAAGCGGATCGGCAACTTCCTTCAACAATCAAAACACCGCAAAAACGGATGTTAATGGTCTGGCGACTTTTGATCTGAAAAGTAGTAAGCAGGAAGACAACACGGTTGAAGTCACCCTTGAAAATGGCGTGAAACAAACGTTAATCGTCAGTTTTGTCGGCGACTCGAGTACTGCGCAGGTTGATCT  >  CP009273/2040651‑2040949

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: