Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 245,306 C→A 93.4% 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,3060CA93.4% 765.0 / 38.9 272pseudogene (1252/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base A (141/113);  total (151/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 272 total aligned reads displayed.

GCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGC  >  CP009273/245159‑245452
                                                                                                                                                   |                                                                                                                                                  
gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAgg                                                                                                                                                  >  2:1139181/1‑150 (MQ=255)
gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAgg                                                                                                                                                  >  2:3720081/1‑150 (MQ=255)
  tCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCGTAGGTc                                                                                                                                                >  2:59860/1‑150 (MQ=255)
    cTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcag                                                                                                                                              <  2:470088/150‑1 (MQ=255)
      gCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcagca                                                                                                                                            <  2:1418266/150‑1 (MQ=255)
      gCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcagca                                                                                                                                            >  2:1557200/1‑150 (MQ=255)
       cTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTcagcag                                                                                                                                           <  2:5076880/150‑1 (MQ=255)
         tGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAgcg                                                                                                                                         <  1:3253196/150‑1 (MQ=255)
            gaacgaacGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGccc                                                                                                                                      <  1:4161372/150‑1 (MQ=255)
              acgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCgg                                                                                                                                    >  2:4736178/1‑150 (MQ=255)
                      aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACttt                                                                                                                            >  2:491554/1‑150 (MQ=255)
                      aTGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACttt                                                                                                                            >  2:3060058/1‑150 (MQ=255)
                          gTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGc                                                                                                                        >  2:4625944/1‑150 (MQ=255)
                          gTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGc                                                                                                                        >  2:465017/1‑150 (MQ=255)
                           tAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGACACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCa                                                                                                                       >  1:3464701/1‑150 (MQ=255)
                             aTGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGACCGGTACACTTTCAGCAAc                                                                                                                     >  1:1921627/1‑150 (MQ=255)
                              tgctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACt                                                                                                                    <  2:485025/150‑1 (MQ=255)
                              tgctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACt                                                                                                                    >  2:539007/1‑150 (MQ=255)
                              tgctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACt                                                                                                                    >  1:21768/1‑150 (MQ=255)
                                ctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTgg                                                                                                                  <  2:4314056/150‑1 (MQ=255)
                                   cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                                                                               >  2:3661934/1‑150 (MQ=255)
                                   cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                                                                               >  1:1562050/1‑150 (MQ=255)
                                     ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCACCTGGCTGTa                                                                                                             >  1:1406620/1‑150 (MQ=255)
                                      cgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                            <  2:2808863/150‑1 (MQ=255)
                                      cgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTaa                                                                                                            >  1:44239/1‑150 (MQ=255)
                                        cAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAtt                                                                                                          >  2:4542060/1‑150 (MQ=255)
                                         aaTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTg                                                                                                         <  1:367649/150‑1 (MQ=255)
                                         aaTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTg                                                                                                         <  1:697685/150‑1 (MQ=255)
                                          aTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGa                                                                                                        >  1:3052515/1‑150 (MQ=255)
                                           tGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAg                                                                                                       >  1:4749106/1‑150 (MQ=255)
                                            gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgc                                                                                                      <  2:4117037/150‑1 (MQ=255)
                                                    cacaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                              <  1:1483035/150‑1 (MQ=255)
                                                       aTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTGAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                           >  2:1777064/1‑150 (MQ=255)
                                                         cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCATGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAAt                                                                                         >  2:4770553/1‑150 (MQ=255)
                                                            ccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCtt                                                                                      >  1:656171/1‑150 (MQ=255)
                                                              gccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTcc                                                                                    >  2:396503/1‑150 (MQ=255)
                                                              gccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTcc                                                                                    >  2:4392671/1‑150 (MQ=255)
                                                               ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                                                                   >  1:52664/1‑150 (MQ=255)
                                                               ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                                                                   >  1:2595971/1‑150 (MQ=255)
                                                                cagcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGc                                                                                  >  1:1112756/1‑150 (MQ=255)
                                                                 agcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                                                                                 >  1:1911801/1‑150 (MQ=255)
                                                                    agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                                                              <  1:54626/150‑1 (MQ=255)
                                                                       atatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                                                                           >  2:2847179/1‑150 (MQ=255)
                                                                        tatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgcg                                                                          <  2:3911132/150‑1 (MQ=255)
                                                                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                         >  1:4825323/1‑150 (MQ=255)
                                                                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                         >  1:2215767/1‑150 (MQ=255)
                                                                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                         <  2:637509/150‑1 (MQ=255)
                                                                           gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGta                                                                       <  2:1328462/150‑2 (MQ=255)
                                                                             tCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGcc                                                                     >  2:2575407/1‑150 (MQ=255)
                                                                              cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                    >  2:1631049/1‑150 (MQ=255)
                                                                              cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                    >  1:446852/1‑150 (MQ=255)
                                                                                 tCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                                                                 >  1:4262043/1‑150 (MQ=255)
                                                                                 tCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                                                                 >  1:4397531/1‑150 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  1:3725910/150‑1 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  1:4983049/150‑1 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  2:5103628/150‑1 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  1:1317292/150‑1 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  2:2014056/150‑1 (MQ=255)
                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                               <  1:515475/150‑1 (MQ=255)
                                                                                           gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg                                                       <  2:4724516/150‑1 (MQ=255)
                                                                                             tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGt                                                     <  2:3320592/150‑1 (MQ=255)
                                                                                               cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                   >  2:4099932/1‑150 (MQ=255)
                                                                                                tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTAt                                                  <  1:2769372/150‑1 (MQ=255)
                                                                                                 ggCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATg                                                 >  2:3836436/1‑150 (MQ=255)
                                                                                                   cGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCa                                               <  1:2642021/150‑1 (MQ=255)
                                                                                                   cGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCa                                               >  2:3665002/1‑150 (MQ=255)
                                                                                                      ccAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACt                                            <  2:1848921/150‑1 (MQ=255)
                                                                                                       cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTg                                           >  1:1139525/1‑150 (MQ=255)
                                                                                                       cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTg                                           <  2:4931393/150‑1 (MQ=255)
                                                                                                       cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTg                                           <  1:3833238/150‑1 (MQ=255)
                                                                                                        aGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGc                                          <  1:3627285/150‑1 (MQ=255)
                                                                                                          cACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGt                                        <  2:3508428/150‑1 (MQ=255)
                                                                                                             ggTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAat                                     >  1:4948839/1‑150 (MQ=255)
                                                                                                              gTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAata                                    <  1:2896316/150‑1 (MQ=255)
                                                                                                                ggCGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatc                                  <  1:3740843/150‑1 (MQ=255)
                                                                                                                  cGATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcat                                >  1:4802178/1‑150 (MQ=255)
                                                                                                                   gATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                               >  1:3192977/1‑150 (MQ=255)
                                                                                                                   gATGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                               >  2:4442739/1‑150 (MQ=255)
                                                                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATAtcatca                              >  1:5071668/1‑150 (MQ=255)
                                                                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATAtcatca                              <  1:4504988/150‑1 (MQ=255)
                                                                                                                      tgTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                            >  1:2465618/2‑150 (MQ=255)
                                                                                                                       gTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAg                           <  2:4221410/150‑1 (MQ=255)
                                                                                                                       gTGACGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCAAAg                           >  2:2485715/1‑150 (MQ=255)
                                                                                                                          aCGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTa                        <  2:2993466/150‑1 (MQ=255)
                                                                                                                           cGATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaa                       >  1:5122184/1‑150 (MQ=255)
                                                                                                                            gATATCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaaa                      <  1:3985241/150‑1 (MQ=255)
                                                                                                                             atatCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAAc                     <  2:2155764/150‑1 (MQ=255)
                                                                                                                              tatCGCGCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACa                    >  2:1009225/1‑150 (MQ=255)
                                                                                                                              tatCGCGCAGGGAAACGCCTTAGGTCAGCAGAGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACa                    >  2:2756787/1‑150 (MQ=255)
                                                                                                                                 cgcgCAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAAt                 <  2:188126/150‑1 (MQ=255)
                                                                                                                                     cAGGGAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAgg             <  1:3109354/150‑1 (MQ=255)
                                                                                                                                       gggAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGaa           <  1:4676854/150‑1 (MQ=255)
                                                                                                                                       gggAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGaa           >  1:2136731/1‑150 (MQ=255)
                                                                                                                                        ggAAACGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat          >  1:1885533/1‑150 (MQ=255)
                                                                                                                                          aaaCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAatat        >  2:3543770/1‑150 (MQ=255)
                                                                                                                                          aaaCGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAatat        >  2:407664/1‑150 (MQ=255)
                                                                                                                                          aaaCGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAatat        >  2:3199288/1‑150 (MQ=255)
                                                                                                                                            aCGCCTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAg      <  2:2083891/150‑1 (MQ=255)
                                                                                                                                                cTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc  >  2:722001/1‑150 (MQ=255)
                                                                                                                                                cTTAGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc  >  1:3063223/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGC  >  CP009273/245159‑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: