Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,097,943 T→A 74.8% N19Y (AAT→TAT)  csgE ← curlin secretion specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,097,9430TA74.8% 414.8 / 154.5 218N19Y (AAT→TAT) csgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base A (71/92);  total (95/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 218 total aligned reads displayed.

CGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCG  >  CP009273/1097795‑1098087
                                                                                                                                                    |                                                                                                                                                
  tGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAccc                                                                                                                                               <  1:3578808/150‑1 (MQ=255)
    cACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTg                                                                                                                                             <  2:2544082/150‑1 (MQ=255)
    cACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTg                                                                                                                                             <  1:1077924/150‑1 (MQ=255)
     aCTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgc                                                                                                                                            <  2:2897566/150‑1 (MQ=255)
      cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                           <  1:4921493/150‑1 (MQ=255)
      cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                           >  1:4932695/1‑150 (MQ=255)
      cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                           >  1:182624/1‑150 (MQ=255)
      cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                           >  2:602818/1‑150 (MQ=255)
       tGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGTCTCGGTAAAAATCAAGGACAATAGATGAAACAGTATGGTCAGTTAGCAAAACCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcgg                                                                                                                                          >  1:1980850/1‑150 (MQ=255)
       tGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcgg                                                                                                                                          <  1:2259796/150‑1 (MQ=255)
         ggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcggcg                                                                                                                                        >  2:779673/1‑150 (MQ=255)
          gCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGa                                                                                                                                       <  1:4219502/150‑1 (MQ=255)
            cTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAAc                                                                                                                                     <  1:4087779/150‑1 (MQ=255)
                    aaTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGGAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTc                                                                                                                             >  2:4539699/1‑150 (MQ=255)
                     aTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCt                                                                                                                            <  2:1709699/150‑1 (MQ=255)
                       cGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTgc                                                                                                                          >  2:3649714/1‑150 (MQ=255)
                        gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgcc                                                                                                                         <  1:4332967/150‑1 (MQ=255)
                        gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgcc                                                                                                                         >  2:1120950/1‑150 (MQ=255)
                        gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTgcc                                                                                                                         <  1:2413885/150‑1 (MQ=255)
                              ttACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCcaca                                                                                                                   <  1:189697/150‑1 (MQ=255)
                              ttACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCcaca                                                                                                                   <  2:2857929/150‑1 (MQ=255)
                               tACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACaa                                                                                                                  >  1:526199/1‑150 (MQ=255)
                                aCCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAAGAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAAt                                                                                                                 <  2:2046411/150‑1 (MQ=255)
                                    gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                             <  2:1368159/150‑1 (MQ=255)
                                    gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                             <  1:3337993/150‑1 (MQ=255)
                                    gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                             >  1:1002402/1‑150 (MQ=255)
                                      ataGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCg                                                                                                           <  2:1615443/150‑1 (MQ=255)
                                      ataGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCg                                                                                                           <  2:638312/150‑1 (MQ=255)
                                               ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACg                                                                                                  >  2:4860865/1‑150 (MQ=255)
                                                tCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGt                                                                                                 >  1:1624791/1‑150 (MQ=255)
                                                    aTTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCa                                                                                             <  2:4058630/150‑1 (MQ=255)
                                                      ttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATg                                                                                           >  2:1532655/1‑150 (MQ=255)
                                                      ttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATg                                                                                           <  2:1030673/150‑1 (MQ=255)
                                                       tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATgg                                                                                          <  2:461474/150‑1 (MQ=255)
                                                         tCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCt                                                                                        <  2:72262/150‑1 (MQ=255)
                                                          cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                       <  2:562780/150‑1 (MQ=255)
                                                           aCTAAAGGCTCGGTAAGAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCttt                                                                                      >  2:4174742/1‑150 (MQ=255)
                                                            cTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTa                                                                                     >  2:1874024/1‑150 (MQ=255)
                                                             tAAAGGCTCGTTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTAt                                                                                    <  2:876926/150‑1 (MQ=255)
                                                             tAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTAt                                                                                    <  1:2208125/150‑1 (MQ=255)
                                                              aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATc                                                                                   <  2:1795382/150‑1 (MQ=255)
                                                              aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATc                                                                                   >  1:4039206/1‑150 (MQ=255)
                                                              aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATc                                                                                   <  2:3059360/150‑1 (MQ=255)
                                                               aaGGCTCGGTAAAAATCATGGACAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCg                                                                                  >  2:4908715/1‑150 (MQ=255)
                                                                    tCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGa                                                                             <  2:2392917/150‑1 (MQ=255)
                                                                      ggTAAAAATCATGGCCAATAGCTGAAACAGTATGGTCAGTTAGCCATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAAGTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                           <  2:4707114/150‑1 (MQ=255)
                                                                      ggTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                           <  2:4347799/150‑1 (MQ=255)
                                                                      ggTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAAGAACGTTTCATGGCTTTATCGCCTGAgg                                                                           <  1:1494754/150‑1 (MQ=255)
                                                                       gTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGt                                                                          <  1:1424149/150‑1 (MQ=255)
                                                                       gTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGt                                                                          <  2:4982549/150‑1 (MQ=255)
                                                                        tAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGtt                                                                         <  2:3563938/150‑1 (MQ=255)
                                                                           aaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATc                                                                      <  2:2137671/150‑1 (MQ=255)
                                                                           aaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATc                                                                      >  1:3847051/1‑150 (MQ=255)
                                                                            aaTCATGGCCAATAGATGAATCAGTATGGTCAGTTAGCTATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCg                                                                     <  2:2654029/150‑1 (MQ=255)
                                                                              tCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGtt                                                                   >  1:2553213/1‑150 (MQ=255)
                                                                                aTGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTg                                                                 >  2:1224798/1‑150 (MQ=255)
                                                                                aTGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTg                                                                 >  2:3159898/1‑150 (MQ=255)
                                                                                    ccAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCa                                                             <  1:1310602/150‑1 (MQ=255)
                                                                                     cAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAg                                                            >  1:1045001/1‑150 (MQ=255)
                                                                                      aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAgg                                                           <  1:1767992/150‑1 (MQ=255)
                                                                                      aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAgg                                                           >  1:1875077/1‑150 (MQ=255)
                                                                                       aTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGa                                                          <  1:4436482/150‑1 (MQ=255)
                                                                                        tAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaa                                                         >  1:3939847/1‑150 (MQ=255)
                                                                                           aTGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAcc                                                      <  1:1616052/150‑1 (MQ=255)
                                                                                               aaCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCtt                                                  >  1:920518/1‑150 (MQ=255)
                                                                                                aCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGAGACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtg                                                 <  1:4797503/150‑1 (MQ=255)
                                                                                                  aGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                               >  2:4292781/1‑150 (MQ=255)
                                                                                                  aGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                               <  1:3370593/150‑1 (MQ=255)
                                                                                                  aGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                               <  2:3674008/150‑1 (MQ=255)
                                                                                                     aTGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTcc                                            >  2:1435965/1‑150 (MQ=255)
                                                                                                      tGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCg                                           <  2:535990/150‑1 (MQ=255)
                                                                                                       ggTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCgg                                          >  1:5009841/1‑150 (MQ=255)
                                                                                                       ggTCAGTTAGCAAGCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCgg                                          <  2:1268161/150‑1 (MQ=255)
                                                                                                         tCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtt                                        <  1:258426/150‑1 (MQ=255)
                                                                                                          cAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGttt                                       >  2:3049069/1‑150 (MQ=255)
                                                                                                                gCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACg                                 >  2:1526786/1‑150 (MQ=255)
                                                                                                                 cAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACgg                                >  1:3954544/1‑150 (MQ=255)
                                                                                                                   aTCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCt                              <  1:1889618/150‑1 (MQ=255)
                                                                                                                    tCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                             <  1:4832342/150‑1 (MQ=255)
                                                                                                                     cccGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAt                            >  1:1055524/1‑150 (MQ=255)
                                                                                                                      ccGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtc                           >  1:2096110/1‑150 (MQ=255)
                                                                                                                      ccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtc                           <  1:4761850/150‑1 (MQ=255)
                                                                                                                       cGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtct                          <  1:4243067/150‑1 (MQ=255)
                                                                                                                        gggACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtctt                         >  1:3609223/1‑150 (MQ=255)
                                                                                                                          gACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttct                       <  1:4600186/150‑1 (MQ=255)
                                                                                                                          gACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttct                       <  2:3352914/150‑1 (MQ=255)
                                                                                                                           aCTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATcttctt                      >  2:3443130/1‑150 (MQ=255)
                                                                                                                           aCTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATcttctt                      <  2:4932695/150‑1 (MQ=255)
                                                                                                                            cTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                     >  2:4302347/1‑150 (MQ=255)
                                                                                                                               cTACCTCTCCGGCGTGAAGAGACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaa                  <  1:4265323/150‑1 (MQ=255)
                                                                                                                               cTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaa                  <  1:1687628/150‑1 (MQ=255)
                                                                                                                                    tCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAtt             >  2:4768258/1‑150 (MQ=255)
                                                                                                                                     cAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATta            <  1:1273124/150‑1 (MQ=255)
                                                                                                                                       aCGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtata          >  2:1367038/1‑150 (MQ=255)
                                                                                                                                       aCGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACCTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtata          >  2:1367195/1‑150 (MQ=255)
                                                                                                                                           cGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGa      >  1:4912319/1‑150 (MQ=255)
                                                                                                                                            gTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGAt     >  2:4056871/1‑150 (MQ=255)
                                                                                                                                            gTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGAt     >  2:3791471/1‑150 (MQ=255)
                                                                                                                                             tGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATg    <  1:81459/150‑1 (MQ=255)
                                                                                                                                               aaGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCCGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCg  <  2:4562438/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCG  >  CP009273/1097795‑1098087

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: