Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 384,783 C→G 6.9% R421G (CGC→GGC)  PP_0320 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947384,7830CG 6.9% 230.3 / 5.8 72R421G (CGC→GGC) PP_0320methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base G (3/2);  total (38/34)
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.67e-01

GAGTACCTCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCC  >  NC_002947/384635‑384929
                                                                                                                                                    |                                                                                                                                                  
gAGTACCTCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGccc                                                                                                                                                    >  2:82601/1‑149 (MQ=255)
  gTACCTCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCg                                                                                                                                                   >  4:244617/1‑148 (MQ=255)
      cTCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATg                                                                                                                                              >  5:172172/1‑149 (MQ=255)
         cGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGcc                                                                                                                                            <  1:133822/148‑1 (MQ=255)
         cGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCg                                                                                                                                           >  2:52370/1‑149 (MQ=255)
             cGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGa                                                                                                                                       <  3:244617/149‑1 (MQ=255)
                 tGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGAAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGACATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCg                                                                                                                                   >  6:169576/1‑149 (MQ=255)
                    gACGGGCCGTGGAGCAGGTAACCCAGGTTGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGcc                                                                                                                                <  1:308151/149‑1 (MQ=255)
                       gggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAg                                                                                                                             >  5:108765/1‑149 (MQ=255)
                        ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGccc                                                                                                                               <  8:6910/146‑1 (MQ=255)
                        ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGccc                                                                                                                               >  7:6910/1‑146 (MQ=255)
                         gCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCgg                                                                                                                                          <  6:120688/134‑1 (MQ=255)
                         gCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCgg                                                                                                                                          >  5:120688/1‑134 (MQ=255)
                          ccGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCg                                                                                                                          >  8:97067/1‑149 (MQ=255)
                           cGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGt                                                                                                                         <  4:285225/149‑1 (MQ=255)
                              ggAGCAGGTAACCCAGGTGGTGGACAACACCGCCGCGCGCTCCCCGCGCAAGCAGGAGATGACCGAAATGGTCGCAAATGCTGAACACGAGATGGGCCTGACCGTGCAGGCCATCTCCCGCAATGCAGGCGATGCGGCCCAGGCGTCGc                                                                                                                      >  6:147417/1‑149 (MQ=255)
                              ggAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCg                                                                                                                       >  6:210848/1‑148 (MQ=255)
                                 gCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGt                                                                                                                   <  1:82601/149‑1 (MQ=255)
                                    ggTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCgg                                                                                                                <  1:52370/149‑1 (MQ=255)
                                    ggTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCgg                                                                                                                >  5:44719/1‑149 (MQ=255)
                                     gTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTc                                                                                                                        <  7:188928/140‑1 (MQ=255)
                                     gTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTc                                                                                                                        >  8:188928/1‑140 (MQ=255)
                                        aCCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCac                                                                                                             <  6:108765/148‑1 (MQ=255)
                                         cccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacga                                                                                                           <  4:217439/149‑1 (MQ=255)
                                          ccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCGGCATTGCGGGCGAt                                                                                                                                      <  3:128507/121‑1 (MQ=255)
                                          ccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCGGCATTGCGGGCGAt                                                                                                                                      >  4:128507/1‑121 (MQ=255)
                                          ccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCac                                                                                                             >  7:112174/1‑146 (MQ=255)
                                           cAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgaga                                                                                                         >  5:160566/1‑149 (MQ=255)
                                                 gtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTc                                                                                                                  >  1:96284/1‑134 (MQ=255)
                                                 gtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTc                                                                                                                  <  2:96284/134‑1 (MQ=255)
                                                 gtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAgg                                                                                                    <  6:130028/148‑1 (MQ=255)
                                                    gACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGTGGCGtt                                                                                                <  8:178589/149‑1 (MQ=255)
                                                    gACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGtt                                                                                                <  5:232897/149‑1 (MQ=255)
                                                        acacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAg                                                                                            <  8:112174/149‑1 (MQ=255)
                                                         cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgagacg                                                                                                       <  7:146782/137‑1 (MQ=255)
                                                         cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgagacg                                                                                                       >  8:146782/1‑137 (MQ=255)
                                                          acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                          >  5:135454/1‑149 (MQ=255)
                                                             gccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCcg                                                                                       >  1:46883/1‑149 (MQ=255)
                                                                gAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGa                                                                                    <  6:172172/149‑1 (MQ=255)
                                                                        cGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGTAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGgtggtg                                                                             >  6:306944/1‑148 (MQ=255)
                                                                        cGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                              <  7:97067/147‑1 (MQ=255)
                                                                        cGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGgtggtg                                                                             >  8:43111/1‑148 (MQ=255)
                                                                           gcgcCAGCAGGAGATGACCGAAATGGTCGCCACCGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAAc                                                                         >  8:186380/1‑149 (MQ=255)
                                                                               cagcagGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATc                                                                     <  7:43111/149‑1 (MQ=255)
                                                                                agcagGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATc                                                                     >  4:309070/1‑148 (MQ=255)
                                                                                   agGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCAtt                                                                 >  1:6539/1‑149 (MQ=255)
                                                                                    tgttttGACCGTAATGGTTGGCACTGCTGTACACGAGTTGGGCCTTACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGTCGTCGCAGTCGGCACTAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTc                                                                <  2:171/144‑1 (MQ=255)
                                                                                       gaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTg                                                             <  6:160566/149‑1 (MQ=255)
                                                                                       gaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTg                                                             <  6:44719/149‑1 (MQ=255)
                                                                                                  tGGTCGCCACTGGTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGGTGCAGGCCCGCGAGGTGGTGCAACGAGCCATTCGTGGCATCGAGGGc                                                  <  5:147417/149‑1 (MQ=255)
                                                                                                  tGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGc                                                  <  5:210848/149‑1 (MQ=255)
                                                                                                               tgTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGc                                         >  1:48822/1‑145 (MQ=255)
                                                                                                               tgTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGc                                         <  2:48822/145‑1 (MQ=255)
                                                                                                                 tACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATc                                   >  6:247463/1‑149 (MQ=255)
                                                                                                                   cacGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCgg                                 >  5:306243/1‑149 (MQ=255)
                                                                                                                    acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACGTCgg                                 >  6:163180/1‑148 (MQ=255)
                                                                                                                    acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCgg                                 >  8:61629/1‑148 (MQ=255)
                                                                                                                      gagaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGc                                                                            <  8:209070/103‑1 (MQ=255)
                                                                                                                      gagaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGc                                                                            >  7:209070/1‑103 (MQ=255)
                                                                                                                           gggCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAggcgg                         >  6:198274/1‑149 (MQ=255)
                                                                                                                            ggCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAggcggc                        <  6:135454/149‑1 (MQ=255)
                                                                                                                              ccTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGcc                       >  8:172136/1‑148 (MQ=255)
                                                                                                                                  aCCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGc                  >  6:59912/1‑149 (MQ=255)
                                                                                                                                    cGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATATCGGGCGACATCGGCAAGGCGGCCGATGCgg                >  3:55120/1‑149 (MQ=255)
                                                                                                                                     gTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGt               <  6:306243/149‑1 (MQ=255)
                                                                                                                                        cAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcag            >  2:217876/1‑149 (MQ=255)
                                                                                                                                        cAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTca             <  2:6539/148‑1 (MQ=255)
                                                                                                                                        cAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTca             <  7:61629/148‑1 (MQ=255)
                                                                                                                                            gcATCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCgg                <  8:100362/140‑1 (MQ=255)
                                                                                                                                            gcATCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCgg                >  7:100362/2‑141 (MQ=255)
                                                                                                                                            gcATCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagccag        >  2:257247/2‑149 (MQ=255)
                                                                                                                                              aTCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGt                                              >  6:23797/1‑109 (MQ=255)
                                                                                                                                              aTCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGt                                              <  5:23797/109‑1 (MQ=255)
                                                                                                                                                   ccGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGcc  >  3:64734/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGTACCTCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCC  >  NC_002947/384635‑384929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑