Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,714 C→T 5.2% A21A (GCC→GCT PP_1087 → OmpA family outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,245,7140CT 5.2% 687.9 / 17.5 211A21A (GCC→GCTPP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (96/104);  new base T (4/7);  total (100/111)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 230 total aligned reads displayed.

AATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGCG  >  NC_002947/1245566‑1245863
                                                                                                                                                    |                                                                                                                                                     
aaTCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccggccg                                                                                                                                                      <  1:3534540/150‑1 (MQ=255)
  tCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCcggccggc                                                                                                                                                    >  1:3754829/1‑150 (MQ=255)
       ccGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAca                                                                                                                                               <  2:4014154/150‑2 (MQ=255)
        cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTg                                                                                                                                              >  2:1867233/1‑150 (MQ=255)
         ggAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGa                                                                                                                                             >  2:1589466/1‑150 (MQ=255)
         ggAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCAGTGCGCCGGCCGGCAGGCTGGCGAAcagg                                                                                                                                             <  2:4271203/150‑4 (MQ=255)
               ccccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTcc                                                                                                                                       <  1:2548425/150‑1 (MQ=255)
               ccccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTcc                                                                                                                                       >  2:3772792/1‑150 (MQ=255)
               ccccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTcc                                                                                                                                       <  1:4163127/150‑1 (MQ=255)
                cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCa                                                                                                                                      >  2:980481/1‑150 (MQ=255)
                cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCa                                                                                                                                      >  1:3409673/1‑150 (MQ=255)
                 ccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAg                                                                                                                                     >  2:523055/1‑150 (MQ=255)
                         aGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCAt                                                                                                                             <  1:2700269/150‑1 (MQ=255)
                         aGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCAt                                                                                                                             <  1:2843743/150‑1 (MQ=255)
                          gCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATg                                                                                                                            <  2:1438455/150‑1 (MQ=255)
                          gCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATg                                                                                                                            <  2:4239529/150‑1 (MQ=255)
                             cGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGgag                                                                                                                         <  2:1414426/150‑1 (MQ=255)
                                cTATCAAGTTTTCCGGACAGCGACCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAAc                                                                                                                      <  2:3590007/150‑1 (MQ=255)
                                 tATCAAGTGTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACa                                                                                                                     <  2:1230297/150‑1 (MQ=255)
                                  aTCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACAt                                                                                                                    >  1:47609/1‑150 (MQ=255)
                                    cAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTg                                                                                                                  <  1:1138495/150‑1 (MQ=255)
                                     aaGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGc                                                                                                                 >  1:90652/1‑150 (MQ=255)
                                       gTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCt                                                                                                               <  1:946560/150‑1 (MQ=255)
                                           tCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGaa                                                                                                           <  1:3109975/150‑1 (MQ=255)
                                             cGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAgg                                                                                                         >  2:1648930/1‑150 (MQ=255)
                                              ggACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGt                                                                                                        >  2:2193034/1‑150 (MQ=255)
                                              ggACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGt                                                                                                        >  1:2440496/1‑150 (MQ=255)
                                               gACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTc                                                                                                       <  1:322677/150‑1 (MQ=255)
                                               gACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTc                                                                                                       >  1:2033584/1‑150 (MQ=255)
                                                  aGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAg                                                                                                    >  2:486345/1‑150 (MQ=255)
                                                   gCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                                   <  2:2117525/150‑1 (MQ=255)
                                                   gCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                                   >  2:1845829/1‑150 (MQ=255)
                                                     gTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGt                                                                                                 >  1:1063356/1‑150 (MQ=255)
                                                      tCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTg                                                                                                >  2:2330576/1‑150 (MQ=255)
                                                         gATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAcc                                                                                             <  1:1665913/150‑1 (MQ=255)
                                                            aGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGt                                                                                          <  2:1401138/150‑1 (MQ=255)
                                                              ccTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTc                                                                                        >  1:3043570/1‑150 (MQ=255)
                                                               ctgctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCg                                                                                       >  1:2287608/1‑150 (MQ=255)
                                                               ctgctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCg                                                                                       >  1:2233807/1‑150 (MQ=255)
                                                                tgctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGt                                                                                                 <  2:852027/139‑1 (MQ=255)
                                                                tgctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGa                                                                                      <  1:4190752/150‑1 (MQ=255)
                                                                  ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgt                                                                                    <  2:1706874/150‑1 (MQ=255)
                                                                  ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgt                                                                                    <  2:3497577/150‑1 (MQ=255)
                                                                  ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGGTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgt                                                                                    <  2:3696892/150‑1 (MQ=255)
                                                                   tgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgtg                                                                                   >  1:492957/1‑150 (MQ=255)
                                                                     gTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtc                                                                                 <  2:2145632/150‑1 (MQ=255)
                                                                     gTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtc                                                                                 <  2:3839283/150‑1 (MQ=255)
                                                                      tAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtcg                                                                                >  2:3267445/1‑150 (MQ=255)
                                                                          ccGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTg                                                                            <  1:2844615/150‑1 (MQ=255)
                                                                          ccGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTg                                                                            >  1:2765485/1‑150 (MQ=255)
                                                                           cGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGa                                                                           <  1:3938777/150‑1 (MQ=255)
                                                                            gTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGag                                                                          >  1:2469160/1‑150 (MQ=255)
                                                                             tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagc                                                                         >  2:4105438/1‑150 (MQ=255)
                                                                             tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagc                                                                         >  2:1917035/1‑150 (MQ=255)
                                                                               aGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagcca                                                                       >  2:2438875/1‑150 (MQ=255)
                                                                                gAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagccag                                                                      >  1:4007488/1‑150 (MQ=255)
                                                                                         gAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                             <  2:1245150/150‑1 (MQ=255)
                                                                                         gAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                             >  1:171645/1‑150 (MQ=255)
                                                                                          aGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAg                                                            >  1:1185683/1‑150 (MQ=255)
                                                                                            ccgccCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGt                                                          >  2:3067065/1‑150 (MQ=255)
                                                                                               cccGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTc                                                       >  2:1985941/1‑150 (MQ=255)
                                                                                                  gTGCGCCAGCGTTACCTAGCCCTGCTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGc                                                    >  1:1405802/1‑150 (MQ=255)
                                                                                                   tGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGACGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCa                                                   >  2:4456957/1‑150 (MQ=255)
                                                                                                    gcgcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAg                                                  <  2:720445/150‑1 (MQ=255)
                                                                                                      gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                                <  1:930277/150‑1 (MQ=255)
                                                                                                      gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                                >  1:1271615/1‑150 (MQ=255)
                                                                                                          gCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTAc                                            >  2:2688038/1‑150 (MQ=255)
                                                                                                                ccTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                      <  1:3667077/150‑1 (MQ=255)
                                                                                                                 cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTg                                     >  1:4504790/1‑150 (MQ=255)
                                                                                                                 cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTg                                     >  1:605180/1‑150 (MQ=255)
                                                                                                                 cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTg                                     >  1:1756389/1‑150 (MQ=255)
                                                                                                                  tAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGc                                    >  1:700977/1‑150 (MQ=255)
                                                                                                                  tAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGc                                    >  2:28925/1‑150 (MQ=255)
                                                                                                                     cccTGTTGACCCTGTTCGCCAGCCTTCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCTGTACGTTCGTGCGTc                                 <  1:324636/150‑1 (MQ=255)
                                                                                                                       cTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCgg                               <  1:3819840/150‑1 (MQ=255)
                                                                                                                        tGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTGTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCggg                              <  2:2941543/150‑1 (MQ=255)
                                                                                                                           tGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAg                           <  2:279501/150‑1 (MQ=255)
                                                                                                                            gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAgg                          <  2:4144759/150‑1 (MQ=255)
                                                                                                                            gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAgg                          >  1:1541447/1‑150 (MQ=255)
                                                                                                                               ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa                       <  1:1867233/150‑1 (MQ=255)
                                                                                                                                 tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGagc                     >  2:2569276/1‑150 (MQ=255)
                                                                                                                                 tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATCGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGagc                     >  1:4214224/1‑150 (MQ=255)
                                                                                                                                   ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGagcag                   >  1:4417696/1‑150 (MQ=255)
                                                                                                                                   ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGagcag                   >  1:648690/1‑150 (MQ=255)
                                                                                                                                   tgcGCCGGCCGGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGagcag                   >  2:795416/3‑150 (MQ=255)
                                                                                                                                          gCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGgtgt            <  2:4203510/150‑1 (MQ=255)
                                                                                                                                            cTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGttt          <  2:3248387/150‑1 (MQ=255)
                                                                                                                                            cTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGttt          <  1:2151673/150‑1 (MQ=255)
                                                                                                                                             tGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTc         <  1:3160429/150‑1 (MQ=255)
                                                                                                                                             tGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTc         >  1:3485355/1‑150 (MQ=255)
                                                                                                                                             tGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTc         >  1:1393120/1‑150 (MQ=255)
                                                                                                                                             tGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTc         >  1:4427679/1‑150 (MQ=255)
                                                                                                                                              gccggccgGCGCACTGACCTCCCAGACGCGCATGGAGACCATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCTATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCa        <  2:4031029/150‑1 (MQ=255)
                                                                                                                                                cggccggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGc      >  2:4514615/1‑150 (MQ=255)
                                                                                                                                                cggccggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGc      >  1:1964580/1‑150 (MQ=255)
                                                                                                                                                  gccggcGAACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGAGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGATg    >  2:1064689/1‑150 (MQ=255)
                                                                                                                                                   ccggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGc   >  1:2375793/1‑150 (MQ=255)
                                                                                                                                                   ccggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGc   <  1:466161/150‑1 (MQ=255)
                                                                                                                                                    cggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGCg  <  1:432365/150‑1 (MQ=255)
                                                                                                                                                    cggcGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGCg  <  2:1257034/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCAGCTGCG  >  NC_002947/1245566‑1245863

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: