Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,885,730 C→A 6.9% R50L (CGC→CTC)  PP_2541 ← transcriptional factor‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,885,7300CA 6.9% 231.3 / 6.0 72R50L (CGC→CTC) PP_2541transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base A (2/3);  total (36/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.42e-01

GACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGGG  >  NC_002947/2885588‑2885875
                                                                                                                                              |                                                                                                                                                 
gacgTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGa                                                                                                                                             >  5:51131/1‑149 (MQ=255)
   gggtgtCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCGCGACCTGGTCGAGCTCGCCgc                                                                                                                                          <  2:191441/147‑2 (MQ=255)
    tgtgtCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGt                                                                                                                                         <  5:325469/149‑1 (MQ=255)
      tgtCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGc                                                                                                                                       >  6:173239/1‑149 (MQ=255)
         cGGTGGTGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGtt                                                                                                                                    <  8:83381/149‑1 (MQ=255)
           gTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTg                                                                                                                                  >  3:78363/1‑149 (MQ=255)
            tGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGc                                                                                                                                 >  1:120291/1‑149 (MQ=255)
              gggTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtg                                                                                                                               <  5:173239/149‑1 (MQ=255)
                gtcaGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAACCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTCAAGCGGGAACAGTGCGGCGAGCTCGACCAGGGCGCGCTCGACGGTGCGTTTGCtggt                                                                                                                             <  1:396655/146‑1 (MQ=255)
                gTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggt                                                                                                                             <  7:316603/149‑1 (MQ=255)
                      gCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATg                                                                                                                       >  4:76170/1‑149 (MQ=255)
                       cGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATgg                                                                                                                      >  4:356365/1‑149 (MQ=255)
                       cGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATg                                                                                                                       <  2:120291/148‑1 (MQ=255)
                            ggTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCAc                                                                                                                 >  4:310011/1‑149 (MQ=255)
                             gTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGCAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACg                                                                                                                >  5:385150/1‑149 (MQ=255)
                                    ctgctgCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGt                                                                                                           >  2:200581/1‑147 (MQ=255)
                                     tgctgcGCTTCGCCCAGGCACAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtg                                                                                                                               <  8:252191/126‑1 (MQ=255)
                                     tgctgcGCTTCGCCCAGGCACAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtg                                                                                                                               >  7:252191/1‑126 (MQ=255)
                                      gctgcGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAgg                                                                                                       >  8:357615/1‑149 (MQ=255)
                                           gcTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGc                                                                                                  <  7:110863/149‑1 (MQ=255)
                                           gcTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGc                                                                                                  >  4:103973/1‑149 (MQ=255)
                                                 cccAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                            >  5:55394/1‑149 (MQ=255)
                                                 cccAGGCGCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTGTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                            >  6:21193/1‑149 (MQ=255)
                                                   cAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACt                                                                                          >  2:21049/1‑149 (MQ=255)
                                                     ggCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTgg                                                                                        <  1:200581/149‑1 (MQ=255)
                                                      gCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTggg                                                                                       <  6:51131/149‑1 (MQ=255)
                                                       cTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGt                                                                                      <  3:103973/149‑1 (MQ=255)
                                                        tCAAGCCCGGTTGCCAGTACCCGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGCCGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                     >  6:239122/1‑149 (MQ=255)
                                                        tCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                     >  3:128847/1‑149 (MQ=255)
                                                          aaGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                     <  4:128847/147‑1 (MQ=255)
                                                          aaGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                     <  6:55394/147‑1 (MQ=255)
                                                          aaGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                     <  4:78363/147‑1 (MQ=255)
                                                             cccGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCt                                                                                <  3:310011/149‑1 (MQ=255)
                                                                  ttGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTc                                                                           >  8:66389/1‑149 (MQ=255)
                                                                     ccAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCg                                                                        >  4:139049/1‑149 (MQ=255)
                                                                      cAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCgg                                                                       >  4:115986/1‑149 (MQ=255)
                                                                       aGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCggg                                                                      <  7:357615/149‑1 (MQ=255)
                                                                          aCCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTgg                                                                   >  8:321496/1‑149 (MQ=255)
                                                                              gCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACt                                                                                          >  6:41378/1‑122 (MQ=255)
                                                                              gCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACt                                                                                          <  5:41378/122‑1 (MQ=255)
                                                                                cATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCg                                                             >  3:348449/1‑149 (MQ=255)
                                                                                 aTAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCACGGGTGGCGGCCgg                                                            >  5:326645/1‑149 (MQ=255)
                                                                                     ggCATGCCTTTACTGTTGCACTGAAGCGGGAACTGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGgcag                                                        >  4:356522/1‑149 (MQ=255)
                                                                                      gCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTg                                                                    <  4:24324/136‑1 (MQ=255)
                                                                                      gCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTg                                                                    >  3:24324/1‑136 (MQ=255)
                                                                                             tttACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACtt                                                  >  4:190220/1‑147 (MQ=255)
                                                                                               tACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGa                                                <  3:190220/147‑1 (MQ=255)
                                                                                               tACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAg                                               <  3:115986/148‑1 (MQ=255)
                                                                                                aCTGTTGCACGGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGGGCGACGGGGCGTTGGCTGGTGGAATGGCCCACGGTGGTAAGGGCGGCGTGCAACTGGGTGGAGCGCATGCCCGGGTGGCGGCCGGGCCGCAACtggctca                                             >  4:215545/1‑143 (MQ=255)
                                                                                                     tGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGcc                                                                                                 <  6:321248/92‑1 (MQ=255)
                                                                                                     tGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGcc                                                                                                 >  5:321248/1‑92 (MQ=255)
                                                                                                             aGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTggcg                                <  4:348449/149‑1 (MQ=255)
                                                                                                              gCGGGAACCGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCGCTggcgg                               <  5:370890/149‑1 (MQ=255)
                                                                                                                   aaCAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGAt                          <  3:356365/149‑1 (MQ=255)
                                                                                                                    aCAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATg                         <  3:76170/149‑1 (MQ=255)
                                                                                                                       gTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgtg                      >  5:78055/1‑149 (MQ=255)
                                                                                                                       gTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgtg                      <  7:66389/149‑1 (MQ=255)
                                                                                                                         ccgtcgAGCGCGACCTGGTCGAGCTCGCCGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATg                         >  8:85593/5‑144 (MQ=255)
                                                                                                                         ccgtcgAGCGCGACCTGGTCGAGCTCGCCGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATg                         <  7:85593/140‑1 (MQ=255)
                                                                                                                            gcgAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgtgt                     <  2:124692/145‑1 (MQ=255)
                                                                                                                            gcgAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgtgt                     >  1:124692/1‑145 (MQ=255)
                                                                                                                              gAGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGtt                                          >  3:113635/1‑122 (MQ=255)
                                                                                                                              gAGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGtt                                          <  4:113635/122‑1 (MQ=255)
                                                                                                                                 cTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcg            >  4:348585/1‑149 (MQ=255)
                                                                                                                                 cTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCGCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGGTGCCACTGGCGGGCGATGGTGTGGcggggcg            <  3:215545/149‑1 (MQ=255)
                                                                                                                                    gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt         <  7:321496/149‑1 (MQ=255)
                                                                                                                                    gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgg           <  3:348585/147‑1 (MQ=255)
                                                                                                                                      ccAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt         >  2:316518/1‑147 (MQ=255)
                                                                                                                                      ccAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt         <  1:316518/147‑1 (MQ=255)
                                                                                                                                         ggTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTg    <  3:139049/149‑1 (MQ=255)
                                                                                                                                          gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCCGCCTGCAACTGGGTGGAGCTCATTCGCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCCCTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTgg   >  2:339127/1‑149 (MQ=255)
                                                                                                                                           tCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTggg  <  1:1175/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGGG  >  NC_002947/2885588‑2885875

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: