Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,107,261 C→A 11.9% D151Y (GAT→TAT)  pppA ← bifunctional prepilin leader peptidase/ methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,107,2610CA11.9% 885.6 / 57.0 268D151Y (GAT→TAT) pppAbifunctional prepilin leader peptidase/ methylase
Reads supporting (aligned to +/- strand):  ref base C (123/112);  new base A (10/22);  total (133/135)
Fisher's exact test for biased strand distribution p-value = 3.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

ACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAAC  >  CP009273/3107112‑3107410
                                                                                                                                                     |                                                                                                                                                     
aCGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATc                                                                                                                                                       <  1:2094933/150‑1 (MQ=255)
   cAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAggg                                                                                                                                                    <  2:229583/150‑1 (MQ=255)
     gAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCa                                                                                                                                                  <  2:4729326/150‑1 (MQ=255)
      aaCTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAGGGCAg                                                                                                                                                 >  1:2186784/1‑150 (MQ=255)
        cTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGcc                                                                                                                                               >  1:3259467/1‑150 (MQ=255)
        cTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGcc                                                                                                                                               >  1:3357101/1‑150 (MQ=255)
           ttCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCAtt                                                                                                                                            <  2:1787199/150‑1 (MQ=255)
              cGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGt                                                                                                                                         <  2:3877361/150‑1 (MQ=255)
               ggCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTg                                                                                                                                        <  1:1986842/150‑1 (MQ=255)
                  tATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATc                                                                                                                                     <  1:3127919/150‑1 (MQ=255)
                  tATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATc                                                                                                                                     >  2:4319271/1‑150 (MQ=255)
                    tCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGa                                                                                                                                   <  1:2806807/150‑1 (MQ=255)
                     ccAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAg                                                                                                                                  >  2:3750618/1‑150 (MQ=255)
                      cAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAgg                                                                                                                                 <  1:4100205/150‑1 (MQ=255)
                      cAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAgg                                                                                                                                 <  1:5073185/150‑1 (MQ=255)
                       aaCGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCAGTGGTGATCGAGGt                                                                                                                                >  1:4501467/1‑150 (MQ=255)
                        aCGCAGAGAGTAAACAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTc                                                                                                                               <  1:1299939/150‑1 (MQ=255)
                        aCGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTCCACACATCGGGCAGCCAGTGGTGAGCGAGGTc                                                                                                                               >  1:676569/1‑150 (MQ=255)
                         cGCAGAGAGTAAACAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCg                                                                                                                              <  2:4868270/150‑1 (MQ=255)
                          gCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGa                                                                                                                             <  1:5097628/150‑1 (MQ=255)
                           cAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGAt                                                                                                                            >  1:1704009/1‑150 (MQ=255)
                             gagagTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCTCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAATCATCGGGCAGCCATTGGTGATCGAGGTCGATgt                                                                                                                          <  2:1608467/150‑2 (MQ=255)
                              agagTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGac                                                                                                                         >  2:2361615/1‑150 (MQ=255)
                               gagTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGaca                                                                                                                        <  1:1021270/150‑1 (MQ=255)
                               gagTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGaca                                                                                                                        <  2:4577723/150‑1 (MQ=255)
                                 gTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACt                                                                                                                      >  1:4222980/1‑150 (MQ=255)
                                 gTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACt                                                                                                                      <  1:2591126/150‑1 (MQ=255)
                                  tAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTc                                                                                                                     >  2:5112696/1‑150 (MQ=255)
                                    aaaaGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTcgc                                                                                                                   >  1:3536776/1‑150 (MQ=255)
                                       aGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGc                                                                                                                >  2:1208168/1‑150 (MQ=255)
                                          gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAAt                                                                                                             >  1:4743448/1‑150 (MQ=255)
                                          gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCACTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAAt                                                                                                             >  2:3208827/1‑150 (MQ=255)
                                            tAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCa                                                                                                           <  2:756074/150‑1 (MQ=255)
                                            tAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCa                                                                                                           <  1:1791423/150‑1 (MQ=255)
                                              aaTCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGc                                                                                                         >  2:4625845/1‑150 (MQ=255)
                                                tCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCa                                                                                                       <  1:3483328/150‑1 (MQ=255)
                                                  ccaccaGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAgg                                                                                                     >  2:2813998/1‑150 (MQ=255)
                                                  ccaccaGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAgg                                                                                                     <  1:2448735/150‑1 (MQ=255)
                                                    accaGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCg                                                                                                   <  1:2393016/150‑1 (MQ=255)
                                                     ccaGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCgg                                                                                                  >  2:2662419/1‑150 (MQ=255)
                                                      caGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGa                                                                                                 <  2:5010895/150‑1 (MQ=255)
                                                      caGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGa                                                                                                 >  1:947496/1‑150 (MQ=255)
                                                      caGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGa                                                                                                 >  1:2999727/1‑150 (MQ=255)
                                                          aCGCCGGTGACTGCATCTTGGAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGataat                                                                                             <  2:2911998/150‑1 (MQ=255)
                                                           cGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGataata                                                                                            >  1:4328065/1‑150 (MQ=255)
                                                             ccGGTGACTGCAGCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATatc                                                                                          <  2:3397494/150‑1 (MQ=255)
                                                                gTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATatcatc                                                                                       <  1:2609088/150‑1 (MQ=255)
                                                                  gACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCACTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCAc                                                                                     >  1:2922924/1‑150 (MQ=255)
                                                                   aCTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCTATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATAATCAcc                                                                                    <  2:4240298/150‑1 (MQ=255)
                                                                   aCTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCAcc                                                                                    <  1:544719/150‑1 (MQ=255)
                                                                     tGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGc                                                                                  <  1:1505077/150‑1 (MQ=255)
                                                                      gCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGcc                                                                                 >  2:4374998/1‑150 (MQ=255)
                                                                        aTCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAg                                                                               >  1:1858282/1‑150 (MQ=255)
                                                                          cTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCAATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGcc                                                                             <  1:3038316/150‑1 (MQ=255)
                                                                                 cGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATcc                                                                      <  2:359132/150‑1 (MQ=255)
                                                                                   tGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCACTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCAc                                                                    >  2:3439607/1‑150 (MQ=255)
                                                                                      gCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACttt                                                                 <  2:4417582/150‑1 (MQ=255)
                                                                                         gACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGACCACATCCACTTTCTg                                                              <  2:3746981/150‑1 (MQ=255)
                                                                                          aCTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTgg                                                             >  1:3935051/1‑150 (MQ=255)
                                                                                            tctGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGcc                                                           >  1:625717/1‑150 (MQ=255)
                                                                                             ctgctgCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGcca                                                          <  2:5257118/150‑1 (MQ=255)
                                                                                               gctgcGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGccaga                                                        >  1:1045226/1‑150 (MQ=255)
                                                                                                     cccATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCcagaccagac                                                  >  2:2112247/1‑150 (MQ=255)
                                                                                                       cATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTc                                                >  2:4740687/1‑150 (MQ=255)
                                                                                                              cAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCa                                         >  1:3865277/1‑150 (MQ=255)
                                                                                                                aTCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaa                                       >  2:5196600/1‑150 (MQ=255)
                                                                                                                   aGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaaa                                    <  2:1482257/150‑1 (MQ=255)
                                                                                                                   aGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaaa                                    >  1:643821/1‑150 (MQ=255)
                                                                                                                   aGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaaa                                    <  2:2403418/150‑1 (MQ=255)
                                                                                                                    gTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAg                                   <  1:3603798/150‑1 (MQ=255)
                                                                                                                      ccTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCa                                 >  2:4686990/1‑150 (MQ=255)
                                                                                                                       cTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCaa                                >  1:3694766/1‑150 (MQ=255)
                                                                                                                        tGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCTCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAg                               <  1:2682915/150‑1 (MQ=255)
                                                                                                                         gTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                              <  1:74670/150‑1 (MQ=255)
                                                                                                                         gTCCACAATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                              <  1:2947805/150‑1 (MQ=255)
                                                                                                                             acaATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTg                          <  2:1960945/150‑1 (MQ=255)
                                                                                                                             acaATACGCCCTGAGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTg                          <  2:3281390/150‑1 (MQ=255)
                                                                                                                              caATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGt                         >  1:2290234/1‑150 (MQ=255)
                                                                                                                               aaTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTc                        <  2:1745091/150‑1 (MQ=255)
                                                                                                                               aaTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTc                        >  1:5162601/1‑150 (MQ=255)
                                                                                                                               aaTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTc                        >  2:1565666/1‑150 (MQ=255)
                                                                                                                                aTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCa                       >  1:3427165/1‑150 (MQ=255)
                                                                                                                                aTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCa                       >  1:1632197/1‑150 (MQ=255)
                                                                                                                                aTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCa                       >  2:2533362/1‑150 (MQ=255)
                                                                                                                                aTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCa                       >  2:3523051/1‑150 (MQ=255)
                                                                                                                                 tACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCaa                      >  2:1396482/1‑150 (MQ=255)
                                                                                                                                   cGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCACGTGCTGTCAATa                    >  1:455325/1‑150 (MQ=255)
                                                                                                                                     cccTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAAc                  <  2:3365030/150‑1 (MQ=255)
                                                                                                                                      ccTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAAcc                 <  1:827144/150‑2 (MQ=255)
                                                                                                                                      ccTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACt                 <  1:5075331/150‑1 (MQ=255)
                                                                                                                                          aGTAAAAACATTGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccac             <  1:317573/150‑1 (MQ=255)
                                                                                                                                          aGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccac             >  1:1400206/1‑150 (MQ=255)
                                                                                                                                          aGTAAAAACATAGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccac             >  1:4629360/1‑150 (MQ=255)
                                                                                                                                           gTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccacc            >  1:1225041/1‑150 (MQ=255)
                                                                                                                                               aaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCg        >  2:1098505/1‑150 (MQ=255)
                                                                                                                                                aaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCgg       >  2:40403/1‑150 (MQ=255)
                                                                                                                                                aaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCgg       >  1:3458119/1‑150 (MQ=255)
                                                                                                                                                  cATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGAt     <  1:1970585/150‑1 (MQ=255)
                                                                                                                                                     cGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAAc  <  1:2361615/150‑1 (MQ=255)
                                                                                                                                                     cGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAAc  <  2:1851526/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAAC  >  CP009273/3107112‑3107410

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: