Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,039,899:1 +G 6.0% coding (712/864 nt) PP_3560 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,039,8991.G 6.0% 336.1 / 14.4 83coding (712/864 nt)PP_3560LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (34/44);  new base G (2/3);  total (36/47)
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.19e-01

TTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGAAGC  >  NC_002947/4039758‑4040047
                                                                                                                                               |                                                                                                                                                    
ttCGACAAGCGTGCCTGGGGCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCg                                                                                                                                              <  1:722670/150‑1 (MQ=255)
   gACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCCGCCATggggat                                                                                                                                            <  2:880471/150‑5 (MQ=255)
   gACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGAt                                                                                                                                           >  1:1533489/1‑150 (MQ=255)
   gACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGAt                                                                                                                                           <  1:1939857/150‑1 (MQ=255)
      aaGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATgac                                                                                                                                        >  1:713183/1‑150 (MQ=255)
      aaGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATgac                                                                                                                                        >  2:2024638/1‑150 (MQ=255)
          gTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGAcc                                                                                                                                    <  2:621236/150‑1 (MQ=255)
           tGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGCAGAATCGTCA‑TTATGGGTCCCTGGCGCGCAGGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCt                                                                                                                                   <  2:82012/150‑1 (MQ=255)
              cTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCCGCCATGG‑GGATgactgcca                                                                                                                                   <  1:381398/147‑6 (MQ=255)
              cTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCCGCCATGG‑GGATgactgcca                                                                                                                                   >  2:381398/1‑142 (MQ=255)
               tGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCaa                                                                                                                               >  1:1186779/1‑150 (MQ=255)
                    cggccggcAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCtt                                                                                                                          <  1:2175787/150‑1 (MQ=255)
                        cggcAGATCGAGCAGGAACGACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCa                                                                                                                      <  2:332098/150‑1 (MQ=255)
                         tgcAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCAc                                                                                                                     >  2:582206/2‑150 (MQ=255)
                              aTCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGcc                                                                                                                >  1:1064221/1‑150 (MQ=255)
                              aTCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGcc                                                                                                                >  1:852136/1‑150 (MQ=255)
                                 gAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGCTGACGGCCGACAACGCTACCCACCCCCCAcc                                                                                                             >  2:1721064/1‑150 (MQ=255)
                                       gAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGTCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCAcctgcctg                                                                                                       >  1:206728/1‑150 (MQ=255)
                                           tACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGt                                                                                                   >  1:1811955/1‑150 (MQ=255)
                                           tACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGt                                                                                                   >  1:619885/1‑150 (MQ=255)
                                           tACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGt                                                                                                   >  1:765221/1‑150 (MQ=255)
                                             cgccgcCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGttt                                                                                                 <  2:2041257/150‑1 (MQ=255)
                                             cgccgcCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGttg                                                                                                 <  2:648817/150‑2 (MQ=255)
                                              gccgccAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTg                                                                                                <  1:699522/150‑1 (MQ=255)
                                              gccgacAGCAGGCGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTg                                                                                                >  2:1171476/1‑150 (MQ=255)
                                                 gccAGCAGGTGCAGGTGAGCGAACTGATGGG‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCg                                                                                             <  1:2404700/150‑1 (MQ=255)
                                                  ccAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGa                                                                                            <  1:1372299/150‑1 (MQ=255)
                                                       aggtgcaggtgAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGc                                                                                       >  2:1009670/1‑150 (MQ=255)
                                                         gtgcaggtgAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCATGGCGCGCGCCACCTGCCTGCCGTTTGGCGAGCCGCAg                                                                                     >  2:563205/1‑150 (MQ=255)
                                                                    cGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCACAGCATCGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCcggcg                                                                          <  1:145986/150‑1 (MQ=255)
                                                                       aCTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggt                                                                       <  1:1029878/150‑1 (MQ=255)
                                                                        cTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtg                                                                      <  2:1811955/150‑1 (MQ=255)
                                                                        cTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtg                                                                      <  2:774754/150‑1 (MQ=255)
                                                                        cTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtg                                                                      <  2:774780/150‑1 (MQ=255)
                                                                         tgctGGA‑ACTGGACAGCCCGGCTGGTCTGGCGCGCTTTGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtgg                                                                     <  2:1839857/147‑1 (MQ=255)
                                                                         tGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtgg                                                                     <  1:1093836/150‑1 (MQ=255)
                                                                          gATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtggt                                                                    <  1:1526932/150‑1 (MQ=255)
                                                                          gATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtggt                                                                    <  2:770209/150‑1 (MQ=255)
                                                                           aTGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTTGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGgtggtg                                                                   <  2:1638728/150‑1 (MQ=255)
                                                                           aTGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGgtggtg                                                                   >  1:498154/1‑150 (MQ=255)
                                                                             ggA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGtt                                                                 >  1:1564283/1‑150 (MQ=255)
                                                                               a‑aCTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGc                                                               >  2:2156001/1‑150 (MQ=255)
                                                                                 aCTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCt                                                              <  1:186017/150‑1 (MQ=255)
                                                                                     gACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAg                                                          >  1:2177670/1‑150 (MQ=255)
                                                                                        aGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCAt                                                       <  1:2024638/150‑1 (MQ=255)
                                                                                          ccAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGa                                                           <  2:1897622/144‑1 (MQ=255)
                                                                                          ccAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGa                                                           >  1:1897622/1‑144 (MQ=255)
                                                                                          ccAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGa                                                     >  1:35484/1‑150 (MQ=255)
                                                                                          ccAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGa                                                     >  2:1033229/1‑150 (MQ=255)
                                                                                               aaGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGaaaaaa                                                <  1:1033229/150‑1 (MQ=255)
                                                                                               aaGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGATGTTGCTGGAGCATGaaaaaa                                                <  1:1171476/150‑1 (MQ=255)
                                                                                                 gCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCa                                                                                                                                                   >  2:1324464/1‑49 (MQ=255)
                                                                                                 gCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCa                                                                                                                                                   <  1:1324464/49‑1 (MQ=255)
                                                                                                 gCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGc                                              <  2:1962156/150‑1 (MQ=255)
                                                                                                   ccTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCtt                                                                                                                          >  1:1738398/1‑72 (MQ=255)
                                                                                                   ccTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCtt                                                                                                                          <  2:1738398/72‑1 (MQ=255)
                                                                                                   ccTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGGGGTGTTGCTGGAGCATGAAAAAAGCCa                                            >  2:1342741/1‑150 (MQ=255)
                                                                                                    cTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAAGGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCCGGAAAAAAGACAg                                           >  1:1101227/1‑150 (MQ=255)
                                                                                                     tGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAgc                                          >  1:2224643/1‑150 (MQ=255)
                                                                                                        cgcgcATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAgcggc                                       <  1:1041137/150‑1 (MQ=255)
                                                                                                         gcgcATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAgcggcg                                      <  1:1622946/150‑1 (MQ=255)
                                                                                                            cATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtggt                                                                    >  1:31598/1‑117 (MQ=255)
                                                                                                            cATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGggtggt                                                                    <  2:31598/117‑1 (MQ=255)
                                                                                                              gggCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCAtt                                 <  1:1721064/149‑1 (MQ=255)
                                                                                                                 cGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGa                              <  1:2005965/150‑1 (MQ=255)
                                                                                                                   gTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGt                            >  1:1586289/1‑150 (MQ=255)
                                                                                                                      aTGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCg                         <  2:1393625/150‑1 (MQ=255)
                                                                                                                      aTGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCg                         >  1:2172529/1‑150 (MQ=255)
                                                                                                                      aTGGGCCTGGGCGTGGCAGTCATCC‑CCATCGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCg                         <  2:455415/150‑1 (MQ=255)
                                                                                                                          gCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGgtgt                     >  2:493703/1‑150 (MQ=255)
                                                                                                                          gCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGCTGGTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGgtgt                     <  2:1533489/150‑1 (MQ=255)
                                                                                                                           ccTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCCTGCCACCCGCCACCTGCCTTCCGTCTGGCGCGAGCCAGCTCTCCGGGCTCATATGACTCGTCCAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGtt                    <  1:646268/150‑1 (MQ=255)
                                                                                                                            cTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAgcggc                                       <  2:187935/130‑1 (MQ=255)
                                                                                                                            cTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAgcggc                                       >  1:187935/1‑130 (MQ=255)
                                                                                                                              gggCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGaa                 <  1:1897016/150‑1 (MQ=255)
                                                                                                                                 cGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGAc              <  2:765221/150‑1 (MQ=255)
                                                                                                                                  gTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGAcc             <  2:1564283/150‑1 (MQ=255)
                                                                                                                                    ggtcgTCATCCGCCATGG‑GGATGACTACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTg                                                             <  1:1028557/96‑1 (MQ=38)
                                                                                                                                    ggtcgTCATCCGCCATGG‑GGATGACTACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTg                                                             >  2:1028557/5‑100 (MQ=38)
                                                                                                                                       aGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCt        <  1:2156001/150‑1 (MQ=255)
                                                                                                                                           aTCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGaa    <  2:35484/150‑1 (MQ=255)
                                                                                                                                           aTCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGaa    >  1:671753/1‑150 (MQ=255)
                                                                                                                                            tCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGAAg   >  2:1488826/1‑150 (MQ=255)
                                                                                                                                             cc‑ccATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGAAGc  >  1:1421663/1‑150 (MQ=255)
                                                                                                                                             cc‑ccATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGAAGc  >  1:1421697/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGA‑ACTGGACAGCCAGGAAGCCCTGGCGCGCATGGCGGTGATGGGCCTGGGCGTGGCAGTCATCC‑CCATGGCGGATGACGACCTGCAACGCTTGCCACCCGCCACCTGCCTGCCGTTTGGCGAGCCGCAGCTGACCCGGCGGGTGGTGTTGCTGGAGCATGAAAAAAGCCAGCGGCGGCATTTGAGTGCGGTGTTGAAGACCGCCCTGGAAGC  >  NC_002947/4039758‑4040047

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