Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,293,265 A→T 5.7% K288* (AAG→TAG) ‡ PP_1128 → OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,293,2650AT 5.7% 261.0 / 3.7 86K288* (AAG→TAG) ‡PP_1128OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (35/46);  new base T (3/2);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGC  >  NC_002947/1293117‑1293413
                                                                                                                                                     |                                                                                                                                                    
cGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCaa                                                                                                                                                     <  1:445509/150‑1 (MQ=255)
cGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCaa                                                                                                                                                     <  1:1347460/150‑1 (MQ=255)
  aGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCCGGCCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGATGCCAAgg                                                                                                                                                   <  1:398378/150‑1 (MQ=255)
   gggCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGc                                                                                                                                                  <  2:1230446/150‑1 (MQ=255)
      ccGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACAAACAGGCCGCCCCCGACTCGCCCCTGTCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGc                                                                                                                                               <  1:233065/150‑1 (MQ=255)
       cgcgCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCt                                                                                                                                              >  1:624144/1‑150 (MQ=255)
           cccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGAAGACCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCtggtg                                                                                                                                          <  1:128074/150‑1 (MQ=255)
            ccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCtggtgg                                                                                                                                         >  2:1926164/1‑150 (MQ=255)
              tGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCa                                                                                                                                       >  2:1588590/1‑150 (MQ=255)
                      cgcgTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGcc                                                                                                                               <  2:1457220/150‑1 (MQ=255)
                       gcgTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGcc                                                                                                                               >  1:480429/1‑149 (MQ=255)
                       gcgTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGcc                                                                                                                               <  2:480429/149‑1 (MQ=255)
                                   tGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgcgg                                                                                                                  <  2:624144/150‑1 (MQ=255)
                                       gACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCTCGGTTATGTTGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaa                                                                                                              <  1:1554621/150‑1 (MQ=255)
                                        aCCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaag                                                                                                             >  1:1428274/1‑150 (MQ=255)
                                         ccccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCtggtg                                                                                                                                          <  2:689370/120‑1 (MQ=255)
                                         ccccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCtggtg                                                                                                                                          >  1:689370/1‑120 (MQ=255)
                                         ccccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaaga                                                                                                            >  1:916768/1‑150 (MQ=255)
                                          cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaagaa                                                                                                           >  1:797845/1‑150 (MQ=255)
                                           ccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGaaa                                                                                                          <  2:916768/150‑1 (MQ=255)
                                           ccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGaaa                                                                                                          >  1:714086/1‑150 (MQ=255)
                                            cGCAAAACACCATGGTCGCGCTGCTGGCCAAATCCCAGCAGGTCGTGCCAACCATAAAGCCGATACGCATGGTGGAGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAAc                                                                                                         <  1:697059/150‑1 (MQ=255)
                                              cAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                       <  2:1255500/150‑1 (MQ=255)
                                              cAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                       <  2:797845/150‑1 (MQ=255)
                                               aaaaCACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTg                                                                                                      <  2:551294/150‑1 (MQ=255)
                                                 aaCACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTa                                                                                                    >  2:866311/1‑150 (MQ=255)
                                                  acacCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAg                                                                                                   <  2:955416/150‑1 (MQ=255)
                                                      cATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCaa                                                                                               >  2:356781/1‑150 (MQ=255)
                                                      cATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCaa                                                                                               >  1:532447/1‑150 (MQ=255)
                                                       aTGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCt                                                                                                                                  >  1:622611/1‑114 (MQ=255)
                                                       aTGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCt                                                                                                                                  <  2:622611/114‑1 (MQ=255)
                                                        tGGTCGCTCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCAGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGCTGGCAAACCTGCCGTCAAGGCTGCGGCGAACAAACCTGTAGCCAAgg                                                                                             <  2:977506/150‑1 (MQ=255)
                                                        tGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                             <  1:461026/150‑1 (MQ=255)
                                                           tCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCa                                                                                          >  2:4033/1‑150 (MQ=255)
                                                                 tgctgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGc                                                                                    >  1:1426743/1‑150 (MQ=255)
                                                                   ctgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAg                                                                                                                         <  2:698866/111‑1 (MQ=255)
                                                                   ctgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGATGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAg                                                                                                                         >  1:698866/1‑111 (MQ=255)
                                                                    tgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTg                                                                                 <  2:624653/150‑1 (MQ=255)
                                                                       ccAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTTGGCAGCtgg                                                                                                                                            >  1:564542/1‑88 (MQ=255)
                                                                       ccAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTTGGCAGCtgg                                                                                                                                            <  2:564542/88‑1 (MQ=255)
                                                                        cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGaa                                                                             >  1:565821/1‑150 (MQ=255)
                                                                        atttggccAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCACGGCCGCGGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCaaaa                                                                                       <  2:1205670/134‑1 (MQ=255)
                                                                        atttggccAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCaaaa                                                                                       >  1:1205670/7‑140 (MQ=255)
                                                                          aaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGa                                                                           <  1:866311/150‑1 (MQ=255)
                                                                           aTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGAc                                                                          >  1:1601859/1‑150 (MQ=255)
                                                                              ccAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGc                                                                       >  1:1844708/1‑150 (MQ=255)
                                                                              ccAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGc                                                                       >  2:1875355/1‑150 (MQ=255)
                                                                               cagcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCa                                                                      <  1:1141405/150‑1 (MQ=255)
                                                                                 gcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGCTGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGc                                                                    >  2:1161888/1‑150 (MQ=255)
                                                                                    ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                 <  2:83104/150‑1 (MQ=255)
                                                                                    ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                 <  2:1884445/150‑1 (MQ=255)
                                                                                    ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTAAGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                 <  1:800999/150‑1 (MQ=255)
                                                                                       cgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGGCGTCAAGGCTGCGGCGAAGAAAc                                                                                                         <  2:940873/107‑1 (MQ=255)
                                                                                       cgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAAc                                                                                                         >  1:940873/1‑107 (MQ=255)
                                                                                           ccGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGaaaaaa                                                          <  2:1428274/150‑1 (MQ=255)
                                                                                            cGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAg                                                                                                   >  1:47778/1‑108 (MQ=255)
                                                                                            cGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAg                                                                                                   <  2:47778/108‑1 (MQ=255)
                                                                                              gCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCt                                                       <  1:971655/150‑1 (MQ=255)
                                                                                                    gCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCg                                                 >  2:166609/1‑150 (MQ=255)
                                                                                                       ccGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGt                                                                                                     >  2:827087/1‑95 (MQ=255)
                                                                                                       ccGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGt                                                                                                     <  1:827087/95‑1 (MQ=255)
                                                                                                       ccGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCCGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGcc                                              >  2:1421843/1‑150 (MQ=255)
                                                                                                             cTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAGAGCTGCGCCGGCCAAGcc                                         >  2:1003850/1‑150 (MQ=255)
                                                                                                               cGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTg                                       <  2:155517/150‑1 (MQ=255)
                                                                                                                      ggcggcGCAGGATGCCAAGGCACCAGCTGCCAAGCCCGCGGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGTCCAAACCGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCaa                                <  1:810091/150‑1 (MQ=255)
                                                                                                                      ggCGTTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTCGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCaa                                <  1:483137/150‑1 (MQ=255)
                                                                                                                          gTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGaaa                            <  1:1173513/150‑1 (MQ=255)
                                                                                                                             cAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCt                         <  1:1503735/150‑1 (MQ=255)
                                                                                                                               ggATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGc                       <  1:795920/150‑1 (MQ=255)
                                                                                                                                gATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGTAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGcc                      >  2:1900415/1‑150 (MQ=255)
                                                                                                                                     cAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 >  2:708217/1‑150 (MQ=255)
                                                                                                                                     cAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGNTGa                 <  2:1194131/150‑1 (MQ=255)
                                                                                                                                      aaGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 <  1:656148/149‑1 (MQ=255)
                                                                                                                                      aaGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 >  2:656148/1‑149 (MQ=255)
                                                                                                                                          cACCAGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 >  2:1582753/1‑145 (MQ=255)
                                                                                                                                          cACCAGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 <  1:1582753/145‑1 (MQ=255)
                                                                                                                                           aCCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGaa           <  2:317834/150‑1 (MQ=255)
                                                                                                                                           aCCAGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGAaaaaaa           >  1:1012571/1‑150 (MQ=255)
                                                                                                                                            ccAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGaaa          <  1:166609/150‑1 (MQ=255)
                                                                                                                                            ccAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGaaa          <  1:1421843/150‑1 (MQ=255)
                                                                                                                                             cAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgc                                                                                                                    <  1:513293/43‑1 (MQ=255)
                                                                                                                                             cAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgc                                                                                                                    >  2:513293/1‑43 (MQ=255)
                                                                                                                                             cAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAg         >  1:1896037/1‑150 (MQ=255)
                                                                                                                                                cTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAg      <  2:1346828/150‑1 (MQ=255)
                                                                                                                                                 tGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAACAGCTGTGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGc     <  2:1601859/150‑1 (MQ=255)
                                                                                                                                                 tGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGc     >  2:1304306/1‑150 (MQ=255)
                                                                                                                                                   ccAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTg   >  2:686046/1‑150 (MQ=255)
                                                                                                                                                    cAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGc  <  2:1426743/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGC  >  NC_002947/1293117‑1293413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: