Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,497,841 T→G 6.6% M141R (ATG→AGG)  PP_1310 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,497,8410TG 6.6% 261.6 / 9.7 91M141R (ATG→AGG) PP_1310hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (48/37);  new base G (3/3);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1497692‑1497990
                                                                                                                                                     |                                                                                                                                                       
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)
                        ccccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGAt                                                                                                                                >  2:1113608/1‑150 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATTGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                    <  2:1349288/130‑5 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATTGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                    >  1:1349288/1‑126 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATg                                                                                                                               >  2:2164653/1‑150 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                    <  2:1348296/130‑5 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                    >  1:1348296/1‑126 (MQ=255)
                          ccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGc                                                                                                                              <  1:1690585/150‑1 (MQ=255)
                               cTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAg                                                                                                                         <  1:1690079/150‑1 (MQ=255)
                                tACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAgc                                                                                                                        >  2:907590/1‑150 (MQ=255)
                                    gCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCTTGACCGCCGAATTCCCGATGCCGCAGCGCTg                                                                                                                    >  1:2187566/1‑150 (MQ=255)
                                    gCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTg                                                                                                                    >  2:175116/1‑150 (MQ=255)
                                      tGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGACGCGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGcc                                                                                                                  >  1:174322/1‑150 (MQ=255)
                                         tGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGc                                                                                                               >  1:1920833/1‑150 (MQ=255)
                                          gCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCt                                                                                                              <  1:1113608/150‑1 (MQ=255)
                                            tGATGGGCTTTTTTTTCGGGGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTcc                                                                                                            >  2:704719/1‑149 (MQ=255)
                                              aTGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCAcc                                                                                                          >  1:672466/1‑150 (MQ=255)
                                                 ggCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtg                                                                                                       <  2:248017/150‑1 (MQ=255)
                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACGGCCGAATTCCCAAGAGGGCAGCGCTGCCCGCTCACCgtgg                                                                                                      >  1:611828/1‑150 (MQ=255)
                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtgg                                                                                                      <  2:1065589/150‑1 (MQ=255)
                                                   cTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtggt                                                                                                     >  2:2213243/1‑150 (MQ=255)
                                                         gttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGcc                                                                                               >  1:1539260/1‑150 (MQ=255)
                                                          ttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGGCg                                                                                              >  1:223321/1‑150 (MQ=255)
                                                          ttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCg                                                                                              <  1:632354/150‑1 (MQ=255)
                                                               tggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGc                                                                                         >  2:1307670/1‑150 (MQ=255)
                                                                  tttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCTTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGa                                                                                      <  2:2187566/150‑1 (MQ=255)
                                                                  tttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGa                                                                                      >  2:1630737/1‑150 (MQ=255)
                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                                                     <  2:567560/150‑1 (MQ=255)
                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                                                     <  2:1920833/150‑1 (MQ=255)
                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTCGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                                                     <  1:1443835/150‑1 (MQ=255)
                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                  >  1:987253/1‑150 (MQ=255)
                                                                       tttccaccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCTCGCTCACCGTGATGCAGCCGGACGCCCATGCAg                                                                                 >  2:1046579/4‑150 (MQ=255)
                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGCACTTGGCGTACCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                                                 >  2:1269791/1‑150 (MQ=255)
                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                                                 <  1:173969/150‑1 (MQ=255)
                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                                                 >  2:208649/1‑150 (MQ=255)
                                                                       ttGCCAACAACCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                                                 >  2:833909/1‑150 (MQ=255)
                                                                           caccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgt                                                                             >  1:753097/1‑150 (MQ=255)
                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                                            >  2:2381330/1‑150 (MQ=255)
                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAATTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                                            >  1:1066541/1‑150 (MQ=255)
                                                                                 ccTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGcgcc                                                                       <  1:1064939/150‑1 (MQ=255)
                                                                                        gggCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTg                                                                <  2:987253/150‑1 (MQ=255)
                                                                                            gcCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGttt                                                            <  2:223321/150‑1 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTaa                                                          >  2:1081619/1‑149 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                          >  1:868342/1‑150 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                          >  2:2256235/1‑150 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                          >  2:167772/1‑150 (MQ=255)
                                                                                                  aaaCTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                      >  1:1477522/1‑150 (MQ=255)
                                                                                                  aaaCTGCAGGCGGCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                      <  2:1018259/150‑1 (MQ=255)
                                                                                                      tGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTc                                                  <  1:175116/150‑1 (MQ=255)
                                                                                                         aGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACACCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGa                                               >  1:1139738/1‑150 (MQ=255)
                                                                                                           gCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACa                                             <  1:208649/150‑1 (MQ=255)
                                                                                                              tCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                      >  1:1361389/1‑138 (MQ=255)
                                                                                                              tCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                      <  2:1361389/138‑1 (MQ=255)
                                                                                                                gAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTg                                        <  1:833909/150‑1 (MQ=255)
                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                                      <  1:1129770/150‑1 (MQ=255)
                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                                      <  1:1307670/150‑1 (MQ=255)
                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                                      <  1:2164653/150‑1 (MQ=255)
                                                                                                                    tGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcgc                                    <  1:167772/150‑1 (MQ=255)
                                                                                                                     gcagcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcgcg                                   >  2:460024/1‑150 (MQ=255)
                                                                                                                           ccTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCa                            >  1:2138642/1‑150 (MQ=255)
                                                                                                                            cTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAg                                            >  1:679838/1‑134 (MQ=255)
                                                                                                                            cTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAg                                            <  2:679838/134‑1 (MQ=255)
                                                                                                                              gTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGt                                       >  1:1408626/1‑137 (MQ=255)
                                                                                                                              gTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGt                                       <  2:1408626/137‑1 (MQ=255)
                                                                                                                               ttCATGATCGCCCCGGCGGTCATGCGCATGACAGCCGAATTCCCGATGCCGCAGGGCTGCCCACTCACCGTGGTGCAGCCGGACGCCGATGTAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAg                        >  2:276483/1‑150 (MQ=255)
                                                                                                                                tCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgc                       >  1:522455/1‑150 (MQ=255)
                                                                                                                                 cATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgct                      >  1:1534392/1‑150 (MQ=255)
                                                                                                                                  aTGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgctg                     >  1:197652/1‑150 (MQ=255)
                                                                                                                                   tGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                  <  1:1158883/89‑1 (MQ=255)
                                                                                                                                   tGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                  >  2:1158883/1‑89 (MQ=255)
                                                                                                                                   tGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgctgc                    >  2:1077060/1‑150 (MQ=255)
                                                                                                                                        gCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGaaa               <  1:1630737/150‑1 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGGCAGCCTGTCGCGTCCCCATGAGCTGCTGAAAg              >  1:1795169/1‑150 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAg              <  2:1534392/150‑1 (MQ=255)
                                                                                                                                          cccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGc                                           <  1:20768/121‑1 (MQ=255)
                                                                                                                                          cccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGc                                           >  2:20768/1‑121 (MQ=255)
                                                                                                                                           ccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGGCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTg            <  1:39209/150‑1 (MQ=255)
                                                                                                                                            cggcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTgg           <  1:2213243/150‑1 (MQ=255)
                                                                                                                                            cggcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCATCCCCATGAGCTGCGGGAAGTgg           >  1:1707284/1‑150 (MQ=255)
                                                                                                                                              gctgTCATGCGTCTGACCGCCGAATTCCCGATGACGCAGGGCTGGCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         <  2:1066541/147‑1 (MQ=255)
                                                                                                                                              gcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTACGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         >  1:1139514/1‑150 (MQ=255)
                                                                                                                                              gcggTCAGGCGCATGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         >  1:2248793/1‑150 (MQ=255)
                                                                                                                                                 gTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCCGTACCGCATGAGCTGCTGAAAGTGGCAGaa      <  1:1269791/150‑1 (MQ=255)
                                                                                                                                                  tCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAAt     >  1:376745/1‑150 (MQ=255)
                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  <  2:672466/150‑1 (MQ=255)
                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  <  2:2138642/150‑1 (MQ=255)
                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  >  2:725205/1‑150 (MQ=255)
                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  <  1:2381330/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1497692‑1497990

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: