Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,434,260 C→A 94.1% P111H (CCT→CAT)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,434,2600CA94.1% 493.7 / 21.5 170P111H (CCT→CAT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base A (76/84);  total (80/90)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 182 total aligned reads displayed.

GTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTC  >  CP009273/4434115‑4434404
                                                                                                                                                 |                                                                                                                                                
 ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                               <  1:2525662/148‑1 (MQ=255)
 ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                               <  1:677045/148‑1 (MQ=255)
 ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                               <  1:969109/148‑1 (MQ=255)
 ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                               <  2:1558299/148‑1 (MQ=255)
  tcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcct                                                                                                                                              >  1:2483899/1‑148 (MQ=255)
  tcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcct                                                                                                                                              >  2:2785884/1‑148 (MQ=255)
   cgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCtt                                                                                                                                             <  2:735033/148‑1 (MQ=255)
        tGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGaa                                                                                                                                        <  2:2279125/148‑1 (MQ=255)
         gAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGAATCCTTCTGAAc                                                                                                                                       >  1:2989647/1‑148 (MQ=255)
          agcagcCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACa                                                                                                                                      <  2:3170223/148‑1 (MQ=255)
             agcCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGt                                                                                                                                   <  1:2759766/148‑1 (MQ=255)
                   gCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTgaaga                                                                                                                             >  1:2486770/1‑148 (MQ=255)
                    cGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGATTGAGTGCCTGTGGAACAGCATTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCATCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaagaa                                                                                                                            <  2:3529409/148‑1 (MQ=255)
                     gTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGaagaag                                                                                                                           <  2:1982649/148‑1 (MQ=255)
                     gTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGaagaag                                                                                                                           <  1:2823495/148‑1 (MQ=255)
                     gTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGaagaag                                                                                                                           <  2:1696049/148‑1 (MQ=255)
                     gTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGaagaag                                                                                                                           <  1:1048482/148‑1 (MQ=255)
                      tGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAagaaga                                                                                                                          <  1:3159629/148‑1 (MQ=255)
                        gCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAgg                                                                                                                        >  1:174023/1‑148 (MQ=255)
                           gTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAg                                                                                                                     <  2:1407299/148‑1 (MQ=255)
                            tAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGa                                                                                                                    <  2:1494642/148‑1 (MQ=255)
                            tAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGa                                                                                                                    <  2:1018092/148‑1 (MQ=255)
                            tAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGa                                                                                                                    >  2:858722/1‑148 (MQ=255)
                             aaaaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGAt                                                                                                                   <  2:2804001/148‑1 (MQ=255)
                               aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAg                                                                                                                 <  1:1848162/148‑1 (MQ=255)
                                aGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGc                                                                                                                >  2:653579/1‑148 (MQ=255)
                                 gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                               <  2:2681628/148‑1 (MQ=255)
                                 gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                               >  2:663934/1‑148 (MQ=255)
                                    ggCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTc                                                                                                            >  2:3373981/1‑148 (MQ=255)
                                        aTCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCt                                                                                                        >  1:782975/1‑148 (MQ=255)
                                        aTCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCt                                                                                                        >  1:2564133/1‑148 (MQ=255)
                                            gCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAt                                                                                                    <  1:1421895/148‑1 (MQ=255)
                                              aTCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtct                                                                                                  <  2:3475836/148‑1 (MQ=255)
                                              aTCTGGCTGTAGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAAGGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtct                                                                                                  <  1:2517368/148‑1 (MQ=255)
                                               tCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtctc                                                                                                 >  1:3067034/1‑148 (MQ=255)
                                                cTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtctct                                                                                                >  1:860844/1‑148 (MQ=255)
                                                 tGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtctctc                                                                                               <  1:2414762/148‑1 (MQ=255)
                                                   gCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCa                                                                                             <  1:1348303/148‑1 (MQ=255)
                                                   gCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCa                                                                                             >  2:432162/1‑148 (MQ=255)
                                                       tCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGcc                                                                                         >  1:344308/1‑148 (MQ=255)
                                                       tCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGcc                                                                                         >  1:2161186/1‑148 (MQ=255)
                                                          ggCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTa                                                                                      >  2:1123175/1‑148 (MQ=255)
                                                            cTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATc                                                                                    <  1:2620852/148‑1 (MQ=255)
                                                             ttGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATcc                                                                                   <  2:1093231/148‑1 (MQ=255)
                                                             ttGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATcc                                                                                   >  1:757283/1‑148 (MQ=255)
                                                               gAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCa                                                                                 >  2:1975301/1‑148 (MQ=255)
                                                                  tGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCa                                                                              >  2:1707466/1‑148 (MQ=255)
                                                                    ccTGTGGAACAGCAGTGTTTGTATTAATGCTCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAcc                                                                            <  1:2197545/148‑1 (MQ=255)
                                                                    ccTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAcc                                                                            <  1:2228783/148‑1 (MQ=255)
                                                                          gAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGaca                                                                      >  2:261213/1‑148 (MQ=255)
                                                                          gAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGaca                                                                      >  1:438981/1‑148 (MQ=255)
                                                                          gAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGaca                                                                      <  2:2962442/148‑1 (MQ=255)
                                                                           aaCAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGacac                                                                     <  1:1816874/148‑1 (MQ=255)
                                                                            aCAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACg                                                                    <  2:3067034/148‑1 (MQ=255)
                                                                                cagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTg                                                                >  2:220501/1‑148 (MQ=255)
                                                                                  gtgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGc                                                              >  2:2660889/1‑148 (MQ=255)
                                                                                   tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCa                                                             >  2:1498861/1‑148 (MQ=255)
                                                                                       tGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTgg                                                         >  1:2276058/1‑148 (MQ=255)
                                                                                       tGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTgg                                                         <  2:879744/148‑1 (MQ=255)
                                                                                          aTTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGTTCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGTCACTTGACa                                                      <  2:814734/148‑1 (MQ=255)
                                                                                          aTTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACa                                                      >  2:1337495/1‑148 (MQ=255)
                                                                                           ttAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACaa                                                     <  1:1937664/148‑1 (MQ=255)
                                                                                               tGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTc                                                 >  1:6405/1‑148 (MQ=255)
                                                                                                 aTCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCaa                                               >  2:2064700/1‑148 (MQ=255)
                                                                                                  tCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAAc                                              <  1:3540486/148‑1 (MQ=255)
                                                                                                  tCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAAc                                              <  2:918652/148‑1 (MQ=255)
                                                                                                    tGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACAt                                            >  1:1592217/1‑148 (MQ=255)
                                                                                                      gcgcTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTa                                          <  2:1623926/148‑1 (MQ=255)
                                                                                                       cgcTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTaa                                         <  2:3420238/148‑1 (MQ=255)
                                                                                                          tGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGAt                                      >  1:970208/1‑148 (MQ=255)
                                                                                                             aaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtga                                   <  1:2390453/148‑1 (MQ=255)
                                                                                                             aaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtga                                   >  2:3323200/1‑148 (MQ=255)
                                                                                                              aGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACTACGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgat                                  >  2:1000518/1‑148 (MQ=255)
                                                                                                              aGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgat                                  >  1:2916885/1‑148 (MQ=255)
                                                                                                              aGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgat                                  >  1:1100995/1‑148 (MQ=255)
                                                                                                               gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatg                                 <  1:2692784/148‑1 (MQ=255)
                                                                                                               gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatg                                 >  1:3301861/1‑148 (MQ=255)
                                                                                                                  aTTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGTACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGaca                              >  2:2429440/1‑148 (MQ=255)
                                                                                                                   ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGacac                             >  1:3481067/1‑148 (MQ=255)
                                                                                                                   ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGacac                             >  1:1034899/1‑148 (MQ=255)
                                                                                                                    tCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacc                            <  1:3128765/148‑1 (MQ=255)
                                                                                                                       ggTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACg                         <  2:1402872/148‑1 (MQ=255)
                                                                                                                       ggTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACg                         <  2:1403307/148‑1 (MQ=255)
                                                                                                                       ggTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACg                         <  1:1404726/148‑1 (MQ=255)
                                                                                                                        gTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                        <  2:971835/148‑1 (MQ=255)
                                                                                                                         tCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGt                       <  1:321878/148‑1 (MQ=255)
                                                                                                                           aaCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTAt                     <  2:2601826/148‑1 (MQ=255)
                                                                                                                              aTCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCgg                  <  1:3391693/148‑1 (MQ=255)
                                                                                                                                cGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg                >  1:795874/1‑148 (MQ=255)
                                                                                                                                cGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg                >  2:1250081/1‑148 (MQ=255)
                                                                                                                                cGACAGTAAAAAAATGCATCCTTCAGAACAGGTTGAAGAAGAGGAAGATATCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCAGACACGGGTGGCACTGGACAACGTCAACATTAACATAGATGACACCACGGTAACGCTg                >  1:2501264/1‑148 (MQ=255)
                                                                                                                                 gACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGa               <  1:1686873/148‑1 (MQ=255)
                                                                                                                                   cAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATg             <  1:2984287/148‑1 (MQ=255)
                                                                                                                                   cAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATg             <  2:2365498/148‑1 (MQ=255)
                                                                                                                                     gTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGa           >  1:2981374/1‑148 (MQ=255)
                                                                                                                                     gTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGa           >  1:1176218/1‑148 (MQ=255)
                                                                                                                                       aaaaaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACt         <  1:2460037/148‑1 (MQ=255)
                                                                                                                                          aaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCa      >  1:348745/1‑148 (MQ=255)
                                                                                                                                           aaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCAc     >  1:2082832/1‑148 (MQ=255)
                                                                                                                                              tGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTc  >  1:2595623/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTC  >  CP009273/4434115‑4434404

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: