Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 245,306 C→A 38.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,3060CA38.0% 293.0 / 207.4 220pseudogene (1252/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (71/65);  new base A (42/42);  total (113/107)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

AGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGC  >  CP009273/245158‑245452
                                                                                                                                                    |                                                                                                                                                  
aGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGACACGCCTTc                                                                                                                                                    >  2:2818074/1‑149 (MQ=255)
    ccTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTc                                                                                                                                                >  1:1570556/1‑149 (MQ=255)
       gCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcagc                                                                                                                                             >  2:807025/1‑149 (MQ=255)
        cTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcagca                                                                                                                                            <  1:3049742/149‑1 (MQ=255)
         ttGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTcagcag                                                                                                                                           >  2:419857/1‑149 (MQ=255)
         ttGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTcagcag                                                                                                                                           >  1:1867470/1‑149 (MQ=255)
             gaacgaacGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGcc                                                                                                                                       <  2:1716316/149‑1 (MQ=255)
                  aacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTa                                                                                                                                  >  2:2424185/1‑149 (MQ=255)
                  aacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTa                                                                                                                                  <  1:3594798/149‑1 (MQ=255)
                  aacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTa                                                                                                                                  >  2:1255594/1‑149 (MQ=255)
                     ggATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTacac                                                                                                                               <  1:3839391/149‑1 (MQ=255)
                       aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACtt                                                                                                                             <  1:3529073/149‑1 (MQ=255)
                       aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACtt                                                                                                                             <  1:159634/149‑1 (MQ=255)
                       aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACtt                                                                                                                             <  1:1477522/149‑1 (MQ=255)
                        tGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACttt                                                                                                                            <  1:3937944/149‑1 (MQ=255)
                            tAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATCTCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGc                                                                                                                        >  2:2098837/1‑149 (MQ=255)
                            tAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGc                                                                                                                        >  1:3718208/1‑149 (MQ=255)
                              aTGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCaa                                                                                                                      <  2:2758599/149‑1 (MQ=255)
                                 ctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTg                                                                                                                   <  2:927353/149‑1 (MQ=255)
                                 ctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTg                                                                                                                   <  2:3756084/149‑1 (MQ=255)
                                   gcgcCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGc                                                                                                                 >  1:3779423/1‑149 (MQ=255)
                                      ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGCTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGGTGGTGACGATATCGCGCAGGGACACGCCTTCGGTCAGCAGCGCCCGGTAGACGTTCAGCAACTGGCTGt                                                                                                              >  2:1517923/1‑149 (MQ=255)
                                      ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGt                                                                                                              <  2:2452801/149‑1 (MQ=255)
                                      ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTGCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGt                                                                                                              >  2:1517912/1‑149 (MQ=255)
                                       cgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTa                                                                                                             <  2:3130622/149‑1 (MQ=255)
                                        gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                            >  1:778066/1‑149 (MQ=255)
                                          aaTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAtt                                                                                                          >  2:3537634/1‑149 (MQ=255)
                                          aaTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAtt                                                                                                          >  2:1036563/1‑149 (MQ=255)
                                            tGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGa                                                                                                        <  2:675387/149‑1 (MQ=255)
                                             gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAg                                                                                                       <  2:1320076/149‑1 (MQ=255)
                                                      acaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                              >  2:2579561/1‑149 (MQ=255)
                                                      acaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                              >  2:1215006/1‑149 (MQ=255)
                                                       caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTc                                                                                             >  2:871162/1‑149 (MQ=255)
                                                       caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCACCTGGCTGTAATTGAGCGCCGCGCTc                                                                                             >  1:232585/1‑149 (MQ=255)
                                                         tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                           >  1:876613/1‑149 (MQ=255)
                                                             ccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCt                                                                                       >  1:2721189/1‑149 (MQ=255)
                                                                   gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                                                                                 >  1:4032350/1‑149 (MQ=255)
                                                                   gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                                                                                 >  1:2609728/1‑149 (MQ=255)
                                                                   gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                                                                                 >  1:195145/1‑149 (MQ=255)
                                                                      gAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                              >  2:1213410/1‑149 (MQ=255)
                                                                       aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAc                                                                             >  1:3107564/1‑149 (MQ=255)
                                                                           tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                         >  1:466662/1‑149 (MQ=255)
                                                                           tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                         >  1:3115233/1‑149 (MQ=255)
                                                                            gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                        <  2:3959917/149‑1 (MQ=255)
                                                                             aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                                                                       <  2:3463/149‑1 (MQ=255)
                                                                             aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                                                                       >  2:386342/1‑149 (MQ=255)
                                                                             aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                                                                       >  1:3407373/1‑149 (MQ=255)
                                                                               cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGcc                                                                     >  1:2302039/1‑149 (MQ=255)
                                                                                   cGTCCCCGTGCTAGGGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAACCGCCTTCGGTCAGCAGCGCCCGGTACACTTACAGCAAGTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                                                                 <  1:1042382/149‑1 (MQ=255)
                                                                                   cGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                                                                 >  1:3698742/1‑149 (MQ=255)
                                                                                    gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcga                                                                <  2:2609728/149‑1 (MQ=255)
                                                                                     tCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  2:1867489/149‑1 (MQ=255)
                                                                                     tCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  1:1162037/149‑1 (MQ=255)
                                                                                      cACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacg                                                              >  2:3260117/1‑149 (MQ=255)
                                                                                       aCCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga                                                             <  2:2714897/149‑1 (MQ=255)
                                                                                       aCCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga                                                             >  2:406415/1‑149 (MQ=255)
                                                                                        ccGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGAt                                                            <  1:3199387/149‑1 (MQ=255)
                                                                                           tGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGTTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATaaa                                                         <  2:1224513/149‑1 (MQ=255)
                                                                                            gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc                                                        >  1:2124241/1‑149 (MQ=255)
                                                                                            gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc                                                        >  1:1645626/1‑149 (MQ=255)
                                                                                            gCTAATGGCGACCAGCTCGGTGGCGTTGGTGACGATATCGCGGAGGGAAACGCCTTAGGTCATCATCGCCCGGTACACTTTCAGCAACTGGGTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc                                                        <  1:3537719/149‑1 (MQ=255)
                                                                                                cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                    >  2:530732/1‑149 (MQ=255)
                                                                                                cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                    <  1:2423989/149‑1 (MQ=255)
                                                                                                cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                    <  2:1375432/149‑1 (MQ=255)
                                                                                                 tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                   <  1:1282888/149‑1 (MQ=255)
                                                                                                 tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                   <  1:201317/149‑1 (MQ=255)
                                                                                                 tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                   <  2:2496008/149‑1 (MQ=255)
                                                                                                   gCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATg                                                 <  1:1160761/149‑1 (MQ=255)
                                                                                                      aCCCGCCCGGGGGCGATGGAGACGATATCGCGCAGGGAAACGCCGTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                              <  2:3831827/149‑1 (MQ=255)
                                                                                                      aCCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                              >  1:733085/1‑149 (MQ=255)
                                                                                                       ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAAc                                             >  1:3906185/1‑149 (MQ=255)
                                                                                                       ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAAc                                             <  2:1955730/149‑1 (MQ=255)
                                                                                                        cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACt                                            >  1:1882247/1‑149 (MQ=255)
                                                                                                              ggTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTaa                                      >  2:568057/1‑149 (MQ=255)
                                                                                                                 ggCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAatat                                   >  1:579413/1‑149 (MQ=255)
                                                                                                                     aTGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACAATAAACGGTTATGCAACTGCGTAATatcatc                               <  2:2300421/149‑1 (MQ=255)
                                                                                                                      tGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATAtcatca                              >  1:3110291/1‑149 (MQ=255)
                                                                                                                      tGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATAtcatca                              >  2:1748162/1‑149 (MQ=255)
                                                                                                                       ggTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACGGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATcatcat                             <  2:3125391/149‑1 (MQ=255)
                                                                                                                        gTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                            <  1:4006389/149‑1 (MQ=255)
                                                                                                                         tGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAg                           <  1:2705327/149‑1 (MQ=255)
                                                                                                                          gACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGt                          <  2:1072517/149‑1 (MQ=255)
                                                                                                                             gATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaa                       <  1:991065/149‑1 (MQ=255)
                                                                                                                               tatCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAAc                     >  2:2927993/1‑149 (MQ=255)
                                                                                                                                atCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCACCTGCGTAATATCATCATAGTTAAACa                    >  2:1568553/1‑149 (MQ=255)
                                                                                                                                   gcgcAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAAt                 <  1:1982330/149‑1 (MQ=255)
                                                                                                                                      cAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAg              <  1:737760/149‑1 (MQ=255)
                                                                                                                                       aGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAgg             >  1:1789591/1‑149 (MQ=255)
                                                                                                                                        gggAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGa            >  2:3785549/1‑149 (MQ=255)
                                                                                                                                          gAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat          <  1:1215006/149‑1 (MQ=255)
                                                                                                                                           aaaCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAata         <  2:184749/149‑1 (MQ=255)
                                                                                                                                            aaCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAatat        >  1:2435098/1‑149 (MQ=255)
                                                                                                                                             aCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAAtata       >  1:3767467/1‑149 (MQ=255)
                                                                                                                                              cGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg      <  2:3436233/149‑1 (MQ=255)
                                                                                                                                               gCCTTCGGTCAGCAGCGCCCGGTAGACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGGGTAATATCATCTTAGTTATACAAATCAGGAATATAGc     <  2:2033755/149‑1 (MQ=255)
                                                                                                                                               gCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCTCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGc     <  2:1536128/149‑1 (MQ=255)
                                                                                                                                                ccTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCt    <  2:876613/149‑1 (MQ=255)
                                                                                                                                                ccTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCt    >  2:3595884/1‑149 (MQ=255)
                                                                                                                                                 cTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTg   <  2:4012101/149‑1 (MQ=255)
                                                                                                                                                  ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc  >  2:2234590/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGC  >  CP009273/245158‑245452

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: