Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,699,478 C→A 11.8% A45D (GCC→GAC) ‡ hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,699,4780CA11.8% 635.3 / 55.9 211A45D (GCC→GAC) ‡hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (78/108);  new base A (11/14);  total (89/122)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 212 total aligned reads displayed.

AAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGAC  >  NC_002947/5699329‑5699626
                                                                                                                                                     |                                                                                                                                                    
   gACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGa                                                                                                                                                   >  2:3886055/1‑150 (MQ=255)
    aCCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                                  <  2:1808938/150‑1 (MQ=255)
    aCCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                                  <  1:3844419/150‑1 (MQ=255)
      cGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGt                                                                                                                                                <  1:717106/150‑1 (MQ=255)
        tATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCGTGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCa                                                                                                                                              >  2:3621236/1‑150 (MQ=255)
         aTTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAc                                                                                                                                             <  1:1176300/150‑1 (MQ=255)
         aTTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAc                                                                                                                                             >  1:399820/1‑150 (MQ=255)
          ttGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAcc                                                                                                                                            >  1:4022941/1‑150 (MQ=255)
           tGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAccc                                                                                                                                           >  1:3356501/1‑150 (MQ=255)
              tCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCa                                                                                                                                        >  2:558619/1‑150 (MQ=255)
               cATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaa                                                                                                                                       <  2:922779/150‑1 (MQ=255)
                aTGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaa                                                                                                                                      >  2:3928991/1‑150 (MQ=255)
                 tGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaa                                                                                                                                     >  2:2809861/1‑150 (MQ=255)
                 tGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaa                                                                                                                                     <  2:618263/150‑1 (MQ=255)
                      aTGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATc                                                                                                                                <  1:422263/150‑1 (MQ=255)
                           tccccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCtgat                                                                                                                           <  1:2139521/149‑1 (MQ=255)
                           tccccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCtgat                                                                                                                           >  1:1676413/2‑150 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                           >  2:1070144/1‑150 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                           <  2:3545901/150‑1 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                           >  1:599573/1‑150 (MQ=255)
                            ccccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCtgatg                                                                                                                          <  1:3190163/150‑1 (MQ=255)
                             cccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgatga                                                                                                                         >  2:4975226/1‑150 (MQ=255)
                              ccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTgatgat                                                                                                                        >  1:1374766/1‑150 (MQ=255)
                               cgtgaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATc                                                                                                                       >  2:844578/7‑150 (MQ=255)
                                 cgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCgg                                                                                                                     >  1:879640/1‑150 (MQ=255)
                                  gaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCgga                                                                                                                    <  1:768841/147‑2 (MQ=255)
                                   agaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGcc                                                                                                                   <  2:3104097/150‑1 (MQ=255)
                                       cGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTAc                                                                                                               >  2:510102/1‑150 (MQ=255)
                                        gTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTAcc                                                                                                              >  1:4705529/1‑150 (MQ=255)
                                         tCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCg                                                                                                             >  1:4471102/1‑150 (MQ=255)
                                           cACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGc                                                                                                           <  1:1780595/150‑1 (MQ=255)
                                            aCGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCg                                                                                                          >  1:4945757/1‑150 (MQ=255)
                                              gAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTc                                                                                                        <  1:2020838/150‑1 (MQ=255)
                                                 cttgaCCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGc                                                                                                     <  1:3787437/146‑1 (MQ=255)
                                                 cTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGc                                                                                                     <  2:1964810/150‑1 (MQ=255)
                                                   tgaCCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCaa                                                                                                   <  1:3683077/148‑1 (MQ=255)
                                                   cAACCGCCACCTCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCaa                                                                                                   <  2:1434691/150‑1 (MQ=255)
                                                         cttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGa                                                                                             >  1:1251938/4‑150 (MQ=255)
                                                          ttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGaa                                                                                            >  2:165518/3‑150 (MQ=255)
                                                           acaTCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGaaa                                                                                           <  2:372428/150‑1 (MQ=255)
                                                           acaTCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGaaa                                                                                           <  1:3697609/150‑1 (MQ=255)
                                                             atcatcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATc                                                                                         <  2:2272823/150‑1 (MQ=255)
                                                               catcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgc                                                                                       >  2:3037752/1‑150 (MQ=255)
                                                               catcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATctc                                                                                       <  2:4725314/150‑3 (MQ=255)
                                                                atcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgcc                                                                                      >  1:1534110/1‑150 (MQ=255)
                                                                 tcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATctccc                                                                                     >  2:3205474/1‑150 (MQ=255)
                                                                    gCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgcccgcc                                                                                  >  2:279919/1‑150 (MQ=255)
                                                                     ccAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCccgtcg                                                                                 >  2:2426375/1‑150 (MQ=255)
                                                                      cAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCccgtcgt                                                                                >  2:1000105/1‑149 (MQ=255)
                                                                       aGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGTGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGc                                                                               >  1:1981949/1‑150 (MQ=255)
                                                                         gacgacGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTg                                                                             <  2:2060827/150‑1 (MQ=255)
                                                                         aacgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCTg                                                                             <  1:1726045/149‑1 (MQ=255)
                                                                         aacgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCTg                                                                             <  2:3273528/149‑1 (MQ=255)
                                                                             acgCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCa                                                                         >  1:662725/1‑150 (MQ=255)
                                                                                  aaGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTg                                                                    <  2:510997/150‑1 (MQ=255)
                                                                                       cgcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAGTTTGGCCag                                                               <  2:1152698/150‑2 (MQ=255)
                                                                                        gcCGTGGCCATTGCCCTGCGCACCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCACAAAGATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGATATCGCCCGCCGGCTGGCCAGTTTGGCCagc                                                              <  2:3021396/150‑3 (MQ=255)
                                                                                        gcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAAc                                                              <  1:4673398/150‑1 (MQ=255)
                                                                                               ccATTGCCCCGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAACCATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGt                                                       <  1:2580056/150‑1 (MQ=255)
                                                                                                 aTTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTc                                                     <  2:4184716/150‑1 (MQ=255)
                                                                                                     cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                                 >  2:1960725/1‑150 (MQ=255)
                                                                                                      ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                                <  2:2822609/150‑1 (MQ=255)
                                                                                                          cgcAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTc                                            <  1:1982161/150‑1 (MQ=255)
                                                                                                             aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                         <  1:3303302/150‑1 (MQ=255)
                                                                                                             aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                         <  1:5046456/150‑1 (MQ=255)
                                                                                                              aCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAgg                                        <  2:1526933/150‑1 (MQ=255)
                                                                                                               ccGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGc                                       <  1:4104516/150‑1 (MQ=255)
                                                                                                                  cTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGccac                                    <  2:650046/150‑1 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 >  1:4788651/1‑150 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 <  2:4284291/150‑1 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 <  2:442268/150‑1 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 >  1:2837601/1‑150 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 >  2:4117326/1‑150 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                 <  1:4376851/150‑1 (MQ=255)
                                                                                                                       ggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGt                               <  1:3484772/150‑1 (MQ=255)
                                                                                                                         cgcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTc                             <  2:2548975/150‑1 (MQ=255)
                                                                                                                          gcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGtt                              >  1:3500551/1‑148 (MQ=255)
                                                                                                                          gcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCa                            >  1:1016921/1‑150 (MQ=255)
                                                                                                                          gcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCa                            >  1:2360949/1‑150 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                           <  1:2974415/150‑1 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                           <  1:2478605/150‑1 (MQ=255)
                                                                                                                            aTGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAcc                          <  1:3058246/150‑1 (MQ=255)
                                                                                                                            aTGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCAAAAAACATCCTGATGATCGGCCCTACCTGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAcc                          <  1:3134195/150‑1 (MQ=255)
                                                                                                                               cAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGaa                       <  2:4238622/150‑1 (MQ=255)
                                                                                                                               cAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGaa                       <  2:1224852/150‑1 (MQ=255)
                                                                                                                                aGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAg                      >  2:2932553/1‑150 (MQ=255)
                                                                                                                                aGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAg                      >  2:1383906/1‑150 (MQ=255)
                                                                                                                                 gCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGt                     <  1:3017695/150‑1 (MQ=255)
                                                                                                                                  cTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCAGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTg                    <  2:4452331/150‑1 (MQ=255)
                                                                                                                                    cccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTggg                  <  2:4581301/150‑1 (MQ=255)
                                                                                                                                     ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAACACATCCTGATGATCGGCCCTACCGGCGTCGACAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc                 <  2:3905455/150‑1 (MQ=255)
                                                                                                                                     ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc                 >  2:2589675/1‑150 (MQ=255)
                                                                                                                                     ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc                 >  1:539070/1‑150 (MQ=255)
                                                                                                                                     ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc                 >  1:243732/1‑150 (MQ=255)
                                                                                                                                          cGAGCTGCGTGCCGTAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTTGGCTACGt            <  1:1128048/150‑1 (MQ=255)
                                                                                                                                          cGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGt            <  2:1407264/150‑1 (MQ=255)
                                                                                                                                          cGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGt            <  2:399820/150‑1 (MQ=255)
                                                                                                                                               tGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGcc       >  1:4299362/1‑150 (MQ=255)
                                                                                                                                                 cGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCcgc     <  2:3418678/150‑1 (MQ=255)
                                                                                                                                                    gCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGAc  <  1:4627910/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGAC  >  NC_002947/5699329‑5699626

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: