Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,699,478 C→A 12.7% 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,4780CA12.7% 739.5 / 77.8 251A45D (GCC→GAC) ‡hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (129/90);  new base A (22/10);  total (151/100)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 255 total aligned reads displayed.

AAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGACG  >  NC_002947/5699329‑5699627
                                                                                                                                                     |                                                                                                                                                     
aaGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGc                                                                                                                                                       >  1:3558046/1‑150 (MQ=255)
    aCCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                                   <  1:1871045/150‑1 (MQ=255)
     ccGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAg                                                                                                                                                  >  1:2722728/1‑150 (MQ=255)
       gTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTc                                                                                                                                                <  2:5203150/150‑1 (MQ=255)
         aTTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAc                                                                                                                                              <  2:2694002/150‑1 (MQ=255)
          ttGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAcc                                                                                                                                             >  2:4220143/1‑150 (MQ=255)
          ttGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAcc                                                                                                                                             >  1:2607187/1‑150 (MQ=255)
                  gTCCATGACCCCCCGCGTGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaaa                                                                                                                                     <  2:877858/150‑1 (MQ=255)
                     cATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACAt                                                                                                                                  >  1:1799761/1‑150 (MQ=255)
                         aCCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtg                                                                                                                              >  2:4583039/1‑150 (MQ=255)
                           tccccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCtgat                                                                                                                            >  1:2580381/2‑150 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                            <  1:911561/150‑1 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                            >  2:2064974/1‑150 (MQ=255)
                           cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                                            >  1:4573243/1‑150 (MQ=255)
                              ccGTGAAATCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTgatgat                                                                                                                         >  2:4488940/1‑150 (MQ=255)
                                gcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCg                                                                                                                       <  1:3681068/150‑1 (MQ=255)
                                  gagaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGc                                                                                                                     <  1:3168401/150‑1 (MQ=255)
                                  gaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCgga                                                                                                                     <  2:4237855/147‑2 (MQ=255)
                                         tCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCg                                                                                                              >  2:74336/1‑150 (MQ=255)
                                           cACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGc                                                                                                            >  2:1624995/1‑150 (MQ=255)
                                            aCGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCg                                                                                                           >  2:1752729/1‑150 (MQ=255)
                                            aCGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCg                                                                                                           >  1:3485422/1‑150 (MQ=255)
                                             cGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGt                                                                                                          >  2:2424392/1‑150 (MQ=255)
                                                  tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                                                     >  1:1688256/1‑150 (MQ=255)
                                                  tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                                                     >  1:306892/1‑150 (MQ=255)
                                                  tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                                                     >  2:3465842/1‑150 (MQ=255)
                                                   cAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCaa                                                                                                    <  1:3007590/150‑1 (MQ=255)
                                                     aCCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGa                                                                                                  >  2:393759/1‑150 (MQ=255)
                                                     aCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCg                                                                                                              >  1:4879317/1‑138 (MQ=255)
                                                      ccGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATcc                                                                                                                                <  1:3767285/119‑1 (MQ=255)
                                                      ccGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATcc                                                                                                                                >  2:3767285/1‑119 (MQ=255)
                                                       cgcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcg                                                                                                >  1:1488113/6‑149 (MQ=255)
                                                       cgcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcg                                                                                                >  2:2551371/6‑149 (MQ=255)
                                                        gcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcgg                                                                                               >  2:1902224/5‑148 (MQ=255)
                                                           acaTCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGaaa                                                                                            <  1:222403/150‑1 (MQ=255)
                                                            catcatCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAAt                                                                                           <  2:546572/150‑1 (MQ=255)
                                                             atcatcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATc                                                                                          >  2:4868491/1‑150 (MQ=255)
                                                              tcatcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcg                                                                                         >  2:1237822/1‑150 (MQ=255)
                                                              tcatcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATct                                                                                         >  1:1403884/1‑149 (MQ=255)
                                                                atcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgcc                                                                                       >  1:3981306/1‑150 (MQ=255)
                                                                atcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATctcc                                                                                       >  1:637722/1‑150 (MQ=255)
                                                                     ccAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCccgccg                                                                                  >  1:5106924/1‑150 (MQ=255)
                                                                     ccAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCccgccg                                                                                  <  2:3027085/150‑1 (MQ=255)
                                                                     ccAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCccgccg                                                                                  >  1:681162/1‑150 (MQ=255)
                                                                     ccAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCccgccg                                                                                  >  1:106029/1‑150 (MQ=255)
                                                                     ccAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCccgtcg                                                                                  >  1:4429517/1‑150 (MQ=255)
                                                                        gaacgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCt                                                                               >  1:88404/3‑150 (MQ=255)
                                                                         gacgacGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTg                                                                              <  1:2082276/150‑1 (MQ=255)
                                                                         gaagacgCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                                   >  2:778026/4‑81 (MQ=255)
                                                                         aacgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCTAGACGGAAATCTCCCGTCGTCTg                                                                              >  2:4803334/2‑150 (MQ=255)
                                                                          gcgccgCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTgg                                                                             >  1:2037067/5‑150 (MQ=255)
                                                                          acgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCTgg                                                                             <  2:2700720/150‑1 (MQ=255)
                                                                               gCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAGAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCaaa                                                                        >  2:3949958/1‑150 (MQ=255)
                                                                                   aGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTgg                                                                    <  1:865751/150‑1 (MQ=255)
                                                                                   aGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTgg                                                                    <  1:3405251/150‑1 (MQ=255)
                                                                                         ccGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAACACATCCGGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAGTTTGGCCagcc                                                              <  2:5014271/150‑4 (MQ=255)
                                                                                         ccGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAAcg                                                              >  2:290674/1‑150 (MQ=255)
                                                                                          cGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAAcgc                                                             <  1:5426653/150‑1 (MQ=255)
                                                                                                 aTTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTc                                                      >  1:749201/1‑150 (MQ=255)
                                                                                                     cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                                  <  1:2759137/150‑1 (MQ=255)
                                                                                                     cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                                  <  2:525701/150‑1 (MQ=255)
                                                                                                     cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                                  >  2:3686336/1‑150 (MQ=255)
                                                                                                      ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                                 <  2:4450240/150‑1 (MQ=255)
                                                                                                      ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                                 >  1:3245835/1‑150 (MQ=255)
                                                                                                      ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                                 >  2:3286872/1‑150 (MQ=255)
                                                                                                        tGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAgg                                               >  1:2477131/1‑150 (MQ=255)
                                                                                                         gcgcAACCGCTGGCGGCGCATGCAGCTCCCCGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGt                                              <  1:1451479/150‑1 (MQ=255)
                                                                                                            cAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGa                                           >  2:2170427/1‑150 (MQ=255)
                                                                                                            cAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGa                                           <  2:530151/150‑1 (MQ=255)
                                                                                                             aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                          >  2:4339029/1‑150 (MQ=255)
                                                                                                             aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                          <  1:2766762/150‑1 (MQ=255)
                                                                                                             aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                          <  1:39405/150‑1 (MQ=255)
                                                                                                                 gCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGcca                                      <  1:5410047/150‑1 (MQ=255)
                                                                                                                 gCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGcca                                      >  2:2826219/1‑150 (MQ=255)
                                                                                                                     gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAATCATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                                  >  2:3867374/1‑150 (MQ=255)
                                                                                                                        gcgcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGtt                               >  2:1249185/1‑150 (MQ=255)
                                                                                                                        gcgcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGtt                               >  2:1549873/1‑150 (MQ=255)
                                                                                                                          gcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCa                             >  1:4150179/1‑150 (MQ=255)
                                                                                                                          gcACGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                                          >  1:5443095/1‑137 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                            <  1:1295488/150‑1 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                            <  2:349088/150‑1 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                            <  1:1947549/150‑1 (MQ=255)
                                                                                                                           cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc                            <  2:1843312/150‑1 (MQ=255)
                                                                                                                            aTGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAcc                           >  2:4356949/1‑150 (MQ=255)
                                                                                                                              gCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGa                         >  1:4605861/1‑150 (MQ=255)
                                                                                                                               cAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGaa                        >  2:3670176/1‑150 (MQ=255)
                                                                                                                                 gCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGt                      <  2:4072899/150‑1 (MQ=255)
                                                                                                                                  cTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTgg                                                                    >  1:284598/1‑103 (MQ=255)
                                                                                                                                     ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc                  <  1:3465860/150‑1 (MQ=255)
                                                                                                                                       tGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTa                >  1:5142758/1‑150 (MQ=255)
                                                                                                                                       tGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTa                <  1:4220143/150‑1 (MQ=255)
                                                                                                                                           gAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGc                                                                   >  2:464387/1‑95 (MQ=255)
                                                                                                                                           gAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTg            >  1:1106170/1‑150 (MQ=255)
                                                                                                                                           gAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTg            >  2:4007679/1‑150 (MQ=255)
                                                                                                                                            aGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTgg           <  2:3834836/150‑1 (MQ=255)
                                                                                                                                               tGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGcc        <  2:5440157/150‑1 (MQ=255)
                                                                                                                                               tGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGcc        >  1:4693489/1‑150 (MQ=255)
                                                                                                                                                  gTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCcgcg     >  1:1612184/1‑150 (MQ=255)
                                                                                                                                                     ccGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGCCg  >  2:303795/1‑150 (MQ=255)
                                                                                                                                                     ccGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGACg  >  2:4958860/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTACGTGGGCCGCGACG  >  NC_002947/5699329‑5699627

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: