Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,293,261 T→G 5.9% A286A (GCT→GCG PP_1128 → OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,293,2610TG 5.9% 232.4 / 4.7 84A286A (GCT→GCGPP_1128OmpA family protein
Reads supporting (aligned to +/- strand):  ref base T (33/46);  new base G (3/2);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

CGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGC  >  NC_002947/1293117‑1293410
                                                                                                                                                 |                                                                                                                                                     
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:1884445/150‑1 (MQ=255)
                                                                                    ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                              <  2:83104/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                                                                                                <  2:47778/108‑1 (MQ=255)
                                                                                            cGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAg                                                                                                >  1:47778/1‑108 (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)
                                                                                                                                                 |                                                                                                                                                     
CGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGC  >  NC_002947/1293117‑1293410

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: