Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 245,306 C→A 46.0% pseudogene (1252/1713 nt) lfhA ← pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273245,3060CA46.0% 133.0 / 217.7 180pseudogene (1252/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (52/45);  new base A (32/51);  total (84/96)
Fisher's exact test for biased strand distribution p-value = 5.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

CTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTG  >  CP009273/245160‑245451
                                                                                                                                                  |                                                                                                                                                 
cTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCg                                                                                                                                                  >  1:177056/1‑148 (MQ=255)
  ccTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGt                                                                                                                                                >  2:1357395/1‑148 (MQ=255)
  ccTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGt                                                                                                                                                <  1:1482750/148‑1 (MQ=255)
   cTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTc                                                                                                                                               <  1:2844579/148‑1 (MQ=255)
    tGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTca                                                                                                                                              >  1:459091/1‑148 (MQ=255)
    tGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTca                                                                                                                                              >  1:1923723/1‑148 (MQ=255)
            aacgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGcc                                                                                                                                      >  1:3493794/1‑148 (MQ=255)
             acgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGccc                                                                                                                                     >  1:971196/1‑148 (MQ=255)
             acgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGccc                                                                                                                                     >  1:1421900/1‑148 (MQ=255)
             acgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGccc                                                                                                                                     >  1:3325058/1‑148 (MQ=255)
              cgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCg                                                                                                                                    >  2:1099490/1‑148 (MQ=255)
                    gATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTacac                                                                                                                              >  1:3361278/1‑148 (MQ=255)
                     aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACt                                                                                                                             <  1:109915/148‑1 (MQ=255)
                     aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACt                                                                                                                             <  2:3307986/148‑1 (MQ=255)
                     aTGGGTAATGCTGCGCCGCAATGCCAGACGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACt                                                                                                                             >  2:1096570/1‑148 (MQ=255)
                      tGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACtt                                                                                                                            <  2:2615471/148‑1 (MQ=255)
                       gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACttt                                                                                                                           >  1:797974/1‑148 (MQ=255)
                             tgctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCaa                                                                                                                     <  1:2593994/148‑1 (MQ=255)
                                  cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGc                                                                                                                <  1:123649/148‑1 (MQ=255)
                                  cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGc                                                                                                                <  2:2343466/148‑1 (MQ=255)
                                   gccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCt                                                                                                               <  2:3659224/148‑1 (MQ=255)
                                      gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTa                                                                                                            >  1:2264596/1‑148 (MQ=255)
                                       cAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                           >  2:2722927/1‑148 (MQ=255)
                                       cAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                           >  1:1594365/1‑148 (MQ=255)
                                       cAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                           >  2:1291531/1‑148 (MQ=255)
                                           gCCAGGCGCACCTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGa                                                                                                       <  2:1594343/148‑1 (MQ=255)
                                           gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGACCCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGa                                                                                                       <  2:1593660/148‑1 (MQ=255)
                                               ggCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgcgc                                                                                                   >  1:3600237/1‑148 (MQ=255)
                                                   cacaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCcgcg                                                                                               <  2:3330447/148‑1 (MQ=255)
                                                    acaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCgcgc                                                                                              >  2:3217469/1‑148 (MQ=255)
                                                     caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                             >  1:838698/1‑148 (MQ=255)
                                                     caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                             >  1:1162995/1‑148 (MQ=255)
                                                     caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                             <  1:435343/148‑1 (MQ=255)
                                                       tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCa                                                                                           >  2:3556249/1‑148 (MQ=255)
                                                       tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCa                                                                                           >  2:257206/1‑148 (MQ=255)
                                                       tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCa                                                                                           >  1:1037567/1‑148 (MQ=255)
                                                       tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCa                                                                                           >  2:2010223/1‑148 (MQ=255)
                                                        cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                          <  2:128630/148‑1 (MQ=255)
                                                        cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                          >  1:3604865/1‑148 (MQ=255)
                                                         ggCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaaa                                                                                         <  1:3820109/148‑1 (MQ=255)
                                                           ccgccAGCAGAAGATGATCTTGCGTCACCGTGCGACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATc                                                                                       <  2:2390880/148‑1 (MQ=255)
                                                             gccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCtt                                                                                     <  2:340654/148‑1 (MQ=255)
                                                                 gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGc                                                                                 >  1:857110/1‑148 (MQ=255)
                                                                    gAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAg                                                                              >  2:715008/1‑148 (MQ=255)
                                                                     aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                             >  1:131204/1‑148 (MQ=255)
                                                                     aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                             <  2:3270639/148‑1 (MQ=255)
                                                                     aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                             >  1:34865/1‑148 (MQ=255)
                                                                     aaTATGATCTTTCGTCACCGTGCTAATGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                             >  1:3400638/1‑148 (MQ=255)
                                                                       tatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcg                                                                           >  2:1987254/1‑148 (MQ=255)
                                                                        atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                                                                          >  1:1186076/1‑148 (MQ=255)
                                                                         tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgcg                                                                         >  2:2797854/1‑148 (MQ=255)
                                                                           aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                       >  2:2774357/1‑148 (MQ=255)
                                                                            tCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                                                                      >  1:663287/1‑148 (MQ=255)
                                                                                tCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGt                                                                  <  1:3433522/148‑1 (MQ=255)
                                                                                  gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                                                                <  1:2336114/148‑1 (MQ=255)
                                                                                      ccGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga                                                            >  2:1145078/1‑148 (MQ=255)
                                                                                       cGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGAt                                                           <  2:1711365/148‑1 (MQ=255)
                                                                                         tGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATaa                                                         >  2:1215431/1‑148 (MQ=255)
                                                                                            tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg                                                      >  2:1874679/1‑148 (MQ=255)
                                                                                            tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg                                                      >  2:1874855/1‑148 (MQ=255)
                                                                                              cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGt                                                    <  2:3325058/148‑1 (MQ=255)
                                                                                              cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGt                                                    <  2:797974/148‑1 (MQ=255)
                                                                                               tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                   <  2:59100/148‑1 (MQ=255)
                                                                                                ggCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                  <  2:1037567/148‑1 (MQ=255)
                                                                                                 gCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTAt                                                 >  2:803653/1‑148 (MQ=255)
                                                                                                     ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                             <  1:749289/148‑1 (MQ=255)
                                                                                                     ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                             <  1:1663063/148‑1 (MQ=255)
                                                                                                     ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                             <  2:2592138/148‑1 (MQ=255)
                                                                                                        gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTg                                          <  2:3079566/148‑1 (MQ=255)
                                                                                                            ggTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTa                                      <  2:2796948/148‑1 (MQ=255)
                                                                                                            ggTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTa                                      <  1:3688152/148‑1 (MQ=255)
                                                                                                             gTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTaa                                     >  2:102597/1‑148 (MQ=255)
                                                                                                             gTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTaa                                     <  2:1834492/148‑1 (MQ=255)
                                                                                                             gTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTaa                                     <  2:1834020/148‑1 (MQ=255)
                                                                                                                gCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAatat                                  <  2:2306436/148‑1 (MQ=255)
                                                                                                                       tGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                           <  2:3775120/148‑1 (MQ=255)
                                                                                                                       tGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                           <  1:2552285/148‑1 (MQ=255)
                                                                                                                        gACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAg                          <  2:3586858/148‑1 (MQ=255)
                                                                                                                          cGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGtt                        <  2:1260857/148‑1 (MQ=255)
                                                                                                                          cGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGtt                        <  2:34899/148‑1 (MQ=255)
                                                                                                                            atatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaa                      <  1:2958086/148‑1 (MQ=255)
                                                                                                                                cgcgCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                  <  1:3217469/148‑1 (MQ=255)
                                                                                                                                cgcgCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                  <  1:3174468/148‑1 (MQ=255)
                                                                                                                                cgcgCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                  <  2:1923723/148‑1 (MQ=255)
                                                                                                                                cgcgCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                  <  2:1594365/148‑1 (MQ=255)
                                                                                                                                 gcgcAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaaa                 <  1:2695892/148‑1 (MQ=255)
                                                                                                                                  cgcAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAAt                <  1:3025618/148‑1 (MQ=255)
                                                                                                                                   gcAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATc               <  2:3455480/148‑1 (MQ=255)
                                                                                                                                   gcAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATc               <  2:1566025/148‑1 (MQ=255)
                                                                                                                                    cAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCa              <  2:3235898/148‑1 (MQ=255)
                                                                                                                                     aGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAg             <  2:493746/148‑1 (MQ=255)
                                                                                                                                         aaaCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat         <  2:971196/148‑1 (MQ=255)
                                                                                                                                         aaaCGCCTTAGGTCAACAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat         <  1:1273144/148‑1 (MQ=255)
                                                                                                                                           aCGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAatat       >  2:1043169/1‑148 (MQ=255)
                                                                                                                                             gCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg     >  2:71960/1‑148 (MQ=255)
                                                                                                                                             gCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg     >  2:2281272/1‑148 (MQ=255)
                                                                                                                                             gCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg     <  1:2799989/148‑1 (MQ=255)
                                                                                                                                             gCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACGGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg     <  2:2894193/148‑1 (MQ=255)
                                                                                                                                              ccTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGc    <  2:3694777/148‑1 (MQ=255)
                                                                                                                                                ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTg  >  2:2638624/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTG  >  CP009273/245160‑245451

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: