Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,699,466 C→G 12.5% A41G (GCC→GGC)  hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,699,4660CG12.5% 762.5 / 82.2 256A41G (GCC→GGC) hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (132/92);  new base G (22/10);  total (154/102)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 256 total aligned reads displayed.

TCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTAC  >  NC_002947/5699319‑5699614
                                                                                                                                                   |                                                                                                                                                    
tCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCg                                                                                                                                                    >  2:1085855/1‑150 (MQ=255)
  aCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                                  >  1:3112342/1‑150 (MQ=255)
  aCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                                  >  2:2858413/1‑150 (MQ=255)
  aCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                                  >  1:252242/1‑150 (MQ=255)
     tGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTg                                                                                                                                               <  1:2307631/150‑1 (MQ=255)
      gCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGc                                                                                                                                              <  1:1041147/150‑1 (MQ=255)
      gCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGc                                                                                                                                              <  2:1162711/150‑1 (MQ=255)
       cacAAGGACCGTAATGCTCATGACCCTGACCCCCCGCGAGAGCGTCCACGACCTCAACCGCCACATCATCGGCCAGGTCGACGCCAAGCGCGCCGAGGCCATTGCCCTGCGCAACCGCTGGCGGCGCCTGCAGCTCCCGGCCGAGCTGCg                                                                                                                                             >  2:3957784/3‑150 (MQ=255)
         cAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTg                                                                                                                                           >  2:2832634/1‑150 (MQ=255)
          aaGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGc                                                                                                                                          >  1:3558046/1‑150 (MQ=255)
              aCCGTATTGTTCATGTCCATGACCCCCCGAGAGATCTTCGACTAACTCTACAGCGACATCATCGGCCAGGACGACTCCAAGCGCGCGGTGGCCATTGGCCTGCGCAACCGCTGGCTGCGCTTTCATCTCGCTGCCTAGATGCGTgaataa                                                                                                                                      <  2:4584528/150‑6 (MQ=255)
              aCCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                      <  2:3665689/150‑1 (MQ=255)
                cGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGt                                                                                                                                    >  1:5440157/1‑150 (MQ=255)
                cGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGt                                                                                                                                    <  2:3758801/150‑1 (MQ=255)
                  tATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCa                                                                                                                                  <  2:4500276/150‑1 (MQ=255)
                   aTTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAc                                                                                                                                 <  1:1896456/150‑1 (MQ=255)
                    ttGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAcc                                                                                                                                >  2:4220143/1‑150 (MQ=255)
                       aTCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAccccc                                                                                                                             >  1:4937676/1‑150 (MQ=255)
                       aTCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAccccc                                                                                                                             >  1:4126628/1‑150 (MQ=255)
                       aTCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAccccc                                                                                                                             >  1:3517002/1‑150 (MQ=255)
                        tCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCa                                                                                                                            >  2:5438820/1‑150 (MQ=255)
                           tGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaa                                                                                                                         <  1:879947/150‑1 (MQ=255)
                            gTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaaa                                                                                                                        >  2:4358796/1‑150 (MQ=255)
                              ccATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTTCCGAAGTCACCCCCAAAAACa                                                                                                                      >  2:2270354/1‑150 (MQ=255)
                              ccATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACa                                                                                                                      >  2:4794812/1‑150 (MQ=255)
                                aTGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                                  >  2:2978205/1‑120 (MQ=255)
                                aTGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                                  <  1:2978205/120‑1 (MQ=255)
                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                               >  2:784838/1‑150 (MQ=255)
                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                               >  2:2064974/1‑150 (MQ=255)
                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                               >  2:5004130/1‑150 (MQ=255)
                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                               >  1:4573243/1‑150 (MQ=255)
                                          gcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAAc                                                                                                                       <  2:2592333/137‑1 (MQ=255)
                                          gcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAAc                                                                                                                       >  1:2592333/1‑137 (MQ=255)
                                          gcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCg                                                                                                          <  1:3681080/150‑1 (MQ=255)
                                            gagaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGc                                                                                                        <  1:3168401/150‑1 (MQ=255)
                                            gaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCgga                                                                                                        <  2:4237855/147‑2 (MQ=255)
                                              gaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGccc                                                                                                      >  2:3244016/1‑150 (MQ=255)
                                                tCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGAGCGGCCCTa                                                                                                    >  2:1734748/1‑150 (MQ=255)
                                                  gTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTAcc                                                                                                  >  1:263867/1‑150 (MQ=255)
                                                   tCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCg                                                                                                 >  2:74336/1‑150 (MQ=255)
                                                     cACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGc                                                                                               >  2:1624995/1‑150 (MQ=255)
                                                          gcttgaCCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTgg                                                                                          <  1:1719710/145‑1 (MQ=255)
                                                           cTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGc                                                                                         >  2:3888105/1‑150 (MQ=255)
                                                            ttgaCCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCa                                                                                        >  2:3480544/4‑150 (MQ=255)
                                                            tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                                        >  2:2766762/1‑150 (MQ=255)
                                                            tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                                        >  2:3465842/1‑150 (MQ=255)
                                                              aaCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAg                                                                                      >  1:3743068/1‑150 (MQ=255)
                                                               aCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCg                                                                                                 <  2:4879317/138‑1 (MQ=255)
                                                                 cgcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcg                                                                                   >  2:2551371/6‑149 (MQ=255)
                                                                 cgcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcg                                                                                   >  1:1488113/6‑149 (MQ=255)
                                                                  gcttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGAcgg                                                                                  >  2:1898712/5‑148 (MQ=255)
                                                                   cttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGa                                                                                 <  2:5029901/147‑1 (MQ=255)
                                                                      catcatCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAAt                                                                              <  2:546572/150‑1 (MQ=255)
                                                                        tcatcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATct                                                                            >  1:1403884/1‑149 (MQ=255)
                                                                          atcGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATctcc                                                                          <  1:4381366/150‑1 (MQ=255)
                                                                           tcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgccc                                                                         <  1:1624995/150‑1 (MQ=255)
                                                                           tcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgccc                                                                         <  1:5425472/150‑1 (MQ=255)
                                                                              gCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATctcccgtc                                                                      >  1:89606/1‑148 (MQ=255)
                                                                               ccAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCccgccg                                                                     >  1:5135473/1‑150 (MQ=255)
                                                                                   gacgacGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                                      <  1:778026/81‑1 (MQ=255)
                                                                                   aacgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCTAGACGGAAATCTCCCGTCGTCTg                                                                 >  2:4803334/2‑150 (MQ=255)
                                                                                    acgacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCTgg                                                                <  2:2700708/150‑1 (MQ=255)
                                                                                         gCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCaaa                                                           >  2:5226035/1‑150 (MQ=255)
                                                                                         gCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCaaa                                                           <  1:3034785/150‑1 (MQ=255)
                                                                                            aaGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTg                                                        <  1:5252754/150‑1 (MQ=255)
                                                                                                      tGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGCAATCGCCCGCCGGCTGGCCAGTTTGGCCagccggc                                              <  1:2812712/150‑7 (MQ=255)
                                                                                                          cATTGCCCTGCGTAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGtt                                          >  2:3948200/1‑150 (MQ=255)
                                                                                                          cATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGtt                                          >  2:880576/1‑150 (MQ=255)
                                                                                                           aTTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTc                                         >  1:1843312/1‑150 (MQ=255)
                                                                                                               cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                     <  1:2759137/150‑1 (MQ=255)
                                                                                                               cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCa                                     <  1:549099/150‑1 (MQ=255)
                                                                                                                ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                    >  1:3245835/1‑150 (MQ=255)
                                                                                                                ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                                    <  2:4450240/150‑1 (MQ=255)
                                                                                                                  tGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAgg                                  >  1:2477131/1‑150 (MQ=255)
                                                                                                                   gcgcAACCGCTGGCGGCGCATGCAGCTCCCCGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGt                                 <  1:1451479/150‑1 (MQ=255)
                                                                                                                      cAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGa                              >  2:2170427/1‑150 (MQ=255)
                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAGAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                             <  1:3949958/150‑1 (MQ=255)
                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                             <  1:654006/150‑1 (MQ=255)
                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                             >  2:5283371/1‑150 (MQ=255)
                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                             >  1:3795744/1‑150 (MQ=255)
                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg                             >  2:4339029/1‑150 (MQ=255)
                                                                                                                          cGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGcc                          >  2:1526554/1‑150 (MQ=255)
                                                                                                                           gCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGcca                         <  2:3035506/150‑1 (MQ=255)
                                                                                                                               gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa                     >  2:386447/1‑150 (MQ=255)
                                                                                                                                  gcgcATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGtt                  >  2:1249185/1‑150 (MQ=255)
                                                                                                                                     cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc               <  1:1295488/150‑1 (MQ=255)
                                                                                                                                     cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc               <  2:1843312/150‑1 (MQ=255)
                                                                                                                                     cATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCAc               <  1:1947558/150‑1 (MQ=255)
                                                                                                                                        gCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGa            >  1:4605455/1‑150 (MQ=255)
                                                                                                                                         cAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGaa           <  2:4025009/150‑1 (MQ=255)
                                                                                                                                         cAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGaa           >  2:3670275/1‑150 (MQ=255)
                                                                                                                                          aGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAg          >  1:10612/1‑150 (MQ=255)
                                                                                                                                          aGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAg          >  1:719836/1‑150 (MQ=255)
                                                                                                                                              cccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTggg      >  1:4723597/1‑150 (MQ=255)
                                                                                                                                               ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc     <  1:3465793/150‑1 (MQ=255)
                                                                                                                                               ccTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGc     <  1:3465860/150‑1 (MQ=255)
                                                                                                                                                 tGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTa   <  1:4220143/150‑1 (MQ=255)
                                                                                                                                                  gCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTAc  >  1:208199/1‑150 (MQ=255)
                                                                                                                                                  gCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTAc  >  1:208434/1‑150 (MQ=255)
                                                                                                                                                  gCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTAc  <  2:4126628/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTCACCGAAGTGGGCTAC  >  NC_002947/5699319‑5699614

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: