Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,165,671 C→T 7.3% G159R (GGG→AGG)  PP_3666 ← MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,165,6710CT 7.3% 194.4 / 4.9 68G159R (GGG→AGG) PP_3666MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGT  >  NC_002947/4165525‑4165818
                                                                                                                                                  |                                                                                                                                                   
cgccAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAgg                                                                                                                                                  >  2:692238/1‑150 (MQ=255)
 gccAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAggc                                                                                                                                                 >  2:879510/1‑150 (MQ=255)
     gcagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGgcggcg                                                                                                                                             <  2:1462024/150‑1 (MQ=255)
      cagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCTGGgcggcgt                                                                                                                                            <  2:520230/150‑2 (MQ=255)
      cagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGgcgc                                                                                                                                            <  1:573795/150‑1 (MQ=255)
       agGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGcgcg                                                                                                                                           >  2:817665/1‑150 (MQ=255)
         gAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACGGCTGAGCAGCTAATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTa                                                                                                                                         <  2:504239/150‑1 (MQ=255)
            cGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTc                                                                                                                                      <  1:1637598/150‑1 (MQ=255)
              gCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCgg                                                                                                                                    >  1:1059736/1‑150 (MQ=255)
               cAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCACGCGCCTCGGGCTCGGc                                                                                                                                   >  2:41659/1‑150 (MQ=255)
                aGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGcc                                                                                                                                  >  1:1083311/1‑150 (MQ=255)
                  gcgcCAGCTCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAg                                                                                                                                <  1:251894/150‑1 (MQ=255)
                       aGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACg                                                                                                                           >  1:1506198/1‑150 (MQ=255)
                         ccccACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGCCGATCGCGGTATGCATCACCGTGCCCTCCTCCGTGCTGGCCGGCACGATGATGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCa                                                                                                                         <  1:1192284/150‑1 (MQ=255)
                           ccACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTAc                                                                                                                                        >  2:157350/1‑133 (MQ=255)
                           ccACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTAc                                                                                                                                        <  1:157350/133‑1 (MQ=255)
                               tCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGcc                                                                                                                   <  1:33485/150‑1 (MQ=255)
                                  tGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCg                                                                                                                <  2:339380/150‑1 (MQ=255)
                                     aGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCTGGgcggcgtggc                                                                                                                                         <  1:274853/122‑5 (MQ=255)
                                     aGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCTGGgcggcgtggc                                                                                                                                         >  2:274853/1‑118 (MQ=255)
                                                   gAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGc                                                                                               <  1:1622303/150‑1 (MQ=255)
                                                   catcagcTCATGCATCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCTCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCTATGCCGCCCAGTCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGc                                                                                               <  1:801798/147‑1 (MQ=255)
                                                    agcagcTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCg                                                                                              <  1:879510/150‑1 (MQ=255)
                                                        gcTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGaaa                                                                                          <  1:971257/150‑1 (MQ=255)
                                                            aTGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAc                                                                                                      <  2:914943/134‑1 (MQ=255)
                                                            aTGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAc                                                                                                      >  1:914943/1‑134 (MQ=255)
                                                            aTGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAccc                                                                                      >  1:1200570/1‑150 (MQ=255)
                                                              gcagcaCCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTTCTAACCGATGATCGGAGGTGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCt                                                                                    <  1:692238/150‑1 (MQ=255)
                                                                agcaCCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTa                                                                                                                                         >  1:1099684/1‑95 (MQ=255)
                                                                agcaCCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTa                                                                                                                                         <  2:1099684/95‑1 (MQ=255)
                                                                  caCCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCa                                                                                                                         >  1:545672/1‑109 (MQ=255)
                                                                  caCCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCa                                                                                                                         <  2:545672/109‑1 (MQ=255)
                                                                      gcgACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAgg                                                                            >  1:447225/1‑150 (MQ=255)
                                                                         acgaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCgg                                                                         >  2:161451/1‑150 (MQ=255)
                                                                           gaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtggc                                                                                                                                         <  1:1499409/84‑5 (MQ=38)
                                                                           gaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtggc                                                                                                                                         >  2:1499409/1‑80 (MQ=38)
                                                                           gaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGAt                                                                       <  1:471783/150‑1 (MQ=255)
                                                                             aCGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATca                                                                     <  2:1506198/150‑1 (MQ=255)
                                                                              cgccgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcag                                                                    <  2:531681/150‑1 (MQ=255)
                                                                                 cgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaa                                                                 <  2:447225/150‑1 (MQ=255)
                                                                                    aCCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAcc                                                                                                     <  2:1636695/111‑1 (MQ=255)
                                                                                    aCCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAcc                                                                                                     >  1:1636695/1‑111 (MQ=255)
                                                                                      cGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacag                                                              <  1:354318/148‑1 (MQ=255)
                                                                                      cGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacag                                                              >  2:354318/1‑148 (MQ=255)
                                                                                      cGAACAACAACCCCGCCACCGTGCTTGCCGGCACGATGCTGGTCTACAAGCCACGCCGCCCAGGCGGCGCGTCCTCGGCCAGGAAAGCAGCAGCCTGAGCGCATGCACAGTACGCGGAAAACCCCTGTACCAGGCGGATcagcaacagca                                                            >  2:358047/1‑150 (MQ=255)
                                                                                          caacaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGc                                                        >  1:412184/1‑150 (MQ=255)
                                                                                            acaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCg                                                      <  1:41659/150‑1 (MQ=255)
                                                                                                 cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGcc                                                 <  2:1200570/150‑1 (MQ=255)
                                                                                                 cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGcc                                                 >  1:134011/1‑150 (MQ=255)
                                                                                                  ccgccgccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCATCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGccc                                                <  1:817665/150‑1 (MQ=255)
                                                                                                   cgccgccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCca                                               >  2:1147807/1‑150 (MQ=255)
                                                                                                        ccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGcc                                          <  2:134011/150‑1 (MQ=255)
                                                                                                              tGGCCGGCACGATGATGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGCACGCAGCGTCCCCGGCGTACTCACCCACGGCGGAAAGCCCCTGTACCAGGCGGATCAGCAACAGCAAGCCCGATGCCCACAGCCCAATCt                                    >  2:688659/1‑150 (MQ=255)
                                                                                                                gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCTTCGTACTCACCCGTAGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc                                  >  1:1349248/1‑150 (MQ=255)
                                                                                                                    gCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCAt                              >  2:690460/1‑150 (MQ=255)
                                                                                                                     cACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATa                             >  1:587044/1‑150 (MQ=255)
                                                                                                                      aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAg                            <  2:1059736/150‑1 (MQ=255)
                                                                                                                        gATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGcc                          >  1:1475177/1‑150 (MQ=255)
                                                                                                                        gATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGcc                          >  1:1255504/1‑150 (MQ=255)
                                                                                                                           gCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCgg                       <  2:1475177/150‑1 (MQ=255)
                                                                                                                             tGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAgg                                                                            >  2:417953/1‑95 (MQ=255)
                                                                                                                             tGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAgg                                                                            <  1:417953/95‑1 (MQ=255)
                                                                                                                               gtgtACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAgg                                                                            >  2:345843/1‑93 (MQ=255)
                                                                                                                               gtgtACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAgg                                                                            <  1:345843/93‑1 (MQ=255)
                                                                                                                                 gtACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCAGGcaacaa                 <  2:587044/150‑1 (MQ=255)
                                                                                                                                   aCAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGc               >  2:927388/1‑150 (MQ=255)
                                                                                                                                     aaGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCGTCATCCCGATCTGCGCATAGCCCGGCAACAAGCCg             >  2:543573/1‑150 (MQ=255)
                                                                                                                                       gCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGAt           >  2:1317333/1‑150 (MQ=255)
                                                                                                                                                cccAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGt  >  2:144284/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGT  >  NC_002947/4165525‑4165818

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: