Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,497,820 T→G 5.0% F134C (TTC→TGC)  PP_1310 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,497,8200TG 5.0% 238.3 / 3.4 81F134C (TTC→TGC) PP_1310hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/34);  new base G (2/2);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CCGATGCCGAGGCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGC  >  NC_002947/1497676‑1497969
                                                                                                                                                |                                                                                                                                                      
ccGATGCCGAGGCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAACCTGCAGGCGTCGAGTTGCAGCACCTGTTCATGa                                                                                                                                                   >  2:2362567/1‑150 (MQ=255)
  gATGCCGAGGCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATc                                                                                                                                                 <  1:1996515/150‑1 (MQ=255)
          ggCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCcggc                                                                                                                                         <  2:347458/150‑1 (MQ=255)
           gCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCcggcg                                                                                                                                        >  2:1064939/1‑150 (MQ=255)
           gCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATAGCCCcggcg                                                                                                                                        >  2:173969/1‑150 (MQ=255)
                tGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAt                                                                                                                                   >  2:1690079/1‑150 (MQ=255)
                    gcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATgcgc                                                                                                                               <  1:1534205/150‑1 (MQ=255)
                    ccaTGGCTCCGTCCCAGGCACCCCGGGCTACCGCTTATGCTGGTGGGCTTATCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGCTGCAGCACCTGTCCATAATCGCCCCGGCGCTCATgccc                                                                                                                               >  2:418778/2‑148 (MQ=255)
                      aTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGgcgcat                                                                                                                             >  2:1763856/1‑148 (MQ=255)
                        ggCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTGGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGa                                                                                                                           <  2:324615/150‑1 (MQ=255)
                           tGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccg                                                                                                                        >  1:107110/1‑150 (MQ=255)
                              gTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccg                                                                                                                     <  2:1368450/150‑1 (MQ=255)
                              gTGCAAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccg                                                                                                                     <  2:1266696/150‑1 (MQ=255)
                               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccgg                                                                                                                    <  1:1763856/150‑2 (MQ=255)
                                        ccccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATc                                                                                                                                                 >  2:16480/1‑112 (MQ=255)
                                        ccccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATc                                                                                                                                                 <  1:16480/112‑1 (MQ=255)
                                        ccccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGAt                                                                                                          >  2:1113608/1‑150 (MQ=255)
                                         cccGGGCTACCGCTGGTGCTGATTGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                              <  2:1349288/130‑5 (MQ=255)
                                         cccGGGCTACCGCTGGTGCTGATTGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                              >  1:1349288/1‑126 (MQ=255)
                                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATg                                                                                                         >  2:2164653/1‑150 (MQ=255)
                                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                              <  2:1348296/130‑5 (MQ=255)
                                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                              >  1:1348296/1‑126 (MQ=255)
                                          ccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGc                                                                                                        <  1:1690585/150‑1 (MQ=255)
                                               cTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAg                                                                                                   <  1:1690079/150‑1 (MQ=255)
                                                tACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAgc                                                                                                  >  2:907590/1‑150 (MQ=255)
                                                    gCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCTTGACCGCCG‑AATTCCCGATGCCGCAGCGCTg                                                                                              >  1:2187566/1‑150 (MQ=255)
                                                    gCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTg                                                                                              >  2:175116/1‑150 (MQ=255)
                                                      tGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGACGCGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGcc                                                                                            >  1:174322/1‑150 (MQ=255)
                                                         tGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGc                                                                                         >  1:1920833/1‑150 (MQ=255)
                                                          gCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCt                                                                                        <  1:1113608/150‑1 (MQ=255)
                                                            tGATGGGCTTTTTTTTCGGGGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTcc                                                                                      >  2:704719/1‑149 (MQ=255)
                                                              aTGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCAcc                                                                                    >  1:672466/1‑150 (MQ=255)
                                                                 ggCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtg                                                                                 <  2:248017/150‑1 (MQ=255)
                                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACGGCCG‑AATTCCCAAGAGGGCAGCGCTGCCCGCTCACCgtgg                                                                                >  1:611828/1‑150 (MQ=255)
                                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtgg                                                                                <  2:1065589/150‑1 (MQ=255)
                                                                   cTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtggt                                                                               >  2:2213243/1‑150 (MQ=255)
                                                                         gttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGcc                                                                         >  1:1539260/1‑150 (MQ=255)
                                                                          ttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGGCg                                                                        >  1:223321/1‑150 (MQ=255)
                                                                          ttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCg                                                                        <  1:632354/150‑1 (MQ=255)
                                                                               tggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGc                                                                   >  2:1307670/1‑150 (MQ=255)
                                                                                  tttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCTTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGa                                                                <  2:2187566/150‑1 (MQ=255)
                                                                                  tttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGa                                                                >  2:1630737/1‑150 (MQ=255)
                                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                               <  2:1920833/150‑1 (MQ=255)
                                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                               <  2:567560/150‑1 (MQ=255)
                                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTCGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                               <  1:1443835/150‑1 (MQ=255)
                                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                            >  1:987253/1‑150 (MQ=255)
                                                                                       tttccaccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCTCGCTCACCGTGATGCAGCCGGACGCCCATGCAg                                                           >  2:1046579/4‑150 (MQ=255)
                                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGCACTTGGCG‑TACCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                           >  2:1269791/1‑150 (MQ=255)
                                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                           >  2:208649/1‑150 (MQ=255)
                                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                           <  1:173969/150‑1 (MQ=255)
                                                                                       ttGCCAACAACCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                           >  2:833909/1‑150 (MQ=255)
                                                                                           caccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgt                                                       >  1:753097/1‑150 (MQ=255)
                                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                      >  2:2381330/1‑150 (MQ=255)
                                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAATTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                      >  1:1066541/1‑150 (MQ=255)
                                                                                                 ccTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGcgcc                                                 <  1:1064939/150‑1 (MQ=255)
                                                                                                        gggCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTg                                          <  2:987253/150‑1 (MQ=255)
                                                                                                            gcCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGttt                                      <  2:223321/150‑1 (MQ=255)
                                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTaa                                    >  2:1081619/1‑149 (MQ=255)
                                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                    >  2:167772/1‑150 (MQ=255)
                                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                    >  1:868342/1‑150 (MQ=255)
                                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                    >  2:2256235/1‑150 (MQ=255)
                                                                                                                  aaaCTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                >  1:1477522/1‑150 (MQ=255)
                                                                                                                  aaaCTGCAGGCGGCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                <  2:1018259/150‑1 (MQ=255)
                                                                                                                      tGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTc                            <  1:175116/150‑1 (MQ=255)
                                                                                                                         aGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCcgg                                                                                                                                          >  1:807896/1‑38 (MQ=255)
                                                                                                                         aGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCcgg                                                                                                                                          <  2:807896/38‑1 (MQ=255)
                                                                                                                         aGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACACCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGa                         >  1:1139738/1‑150 (MQ=255)
                                                                                                                           gCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACa                       <  1:208649/150‑1 (MQ=255)
                                                                                                                              tCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                >  1:1361389/1‑138 (MQ=255)
                                                                                                                              tCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                <  2:1361389/138‑1 (MQ=255)
                                                                                                                                gAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTg                  <  1:833909/150‑1 (MQ=255)
                                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                <  1:2164653/150‑1 (MQ=255)
                                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                <  1:1307670/150‑1 (MQ=255)
                                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                <  1:1129770/150‑1 (MQ=255)
                                                                                                                                    tGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcgc              <  1:167772/150‑1 (MQ=255)
                                                                                                                                     gcagcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcgcg             >  2:460024/1‑150 (MQ=255)
                                                                                                                                           ccTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCa       >  1:2138642/1‑150 (MQ=255)
                                                                                                                                            cTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAg                      <  2:679838/134‑1 (MQ=255)
                                                                                                                                            cTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAg                      >  1:679838/1‑134 (MQ=255)
                                                                                                                                              gTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGt                 <  2:1408626/137‑1 (MQ=255)
                                                                                                                                              gTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGt                 >  1:1408626/1‑137 (MQ=255)
                                                                                                                                               ttCATGATCGCCCCGGCGGTCATGCGCATGACAGCCG‑AATTCCCGATGCCGCAGGGCTGCCCACTCACCGTGGTGCAGCCGGACGCCGATGTAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAg   >  2:276483/1‑150 (MQ=255)
                                                                                                                                                tCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgc  >  1:522455/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CCGATGCCGAGGCCGCTGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCG‑AATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGC  >  NC_002947/1497676‑1497969

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: ‑