Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,632,620 C→G 6.9% P31R (CCC→CGC)  PP_3201 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,632,6200CG 6.9% 218.6 / 4.4 69P31R (CCC→CGC) PP_3201BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/3);  total (34/35)
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.97e-01

CAGGCGCGCTGCGCAGGTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGCCACTGCCGT  >  NC_002947/3632476‑3632761
                                                                                                                                                  |                                                                                                                                             
cAGGCGCGCTGCGCAGGTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCcgcg                                                                                                                                           >  7:36185/1‑149 (MQ=255)
   gcgcgcTGCGCAGGTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTg                                                                                                                                        >  2:80900/1‑149 (MQ=255)
     gcgcTGCGCAGGTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCa                                                                                                                                      >  7:87724/1‑149 (MQ=255)
        cTGCGCAGGTTCTGTACCGTCTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGTGCAGCGCGGGCTGcgcga                                                                                                                                   <  3:60709/149‑2 (MQ=255)
            gcAGGTGCTGTACAGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCTCAGCCAGCGCTGCACGGcccc                                                                                                                               >  8:88517/1‑149 (MQ=255)
                gTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAgg                                                                                                                           >  4:108820/1‑149 (MQ=255)
                   cTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCCCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGctc                                                                                                                        <  5:80805/149‑1 (MQ=255)
                         cGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTACCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgcgtg                                                                                                                  >  5:53576/1‑149 (MQ=255)
                          gACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgcgtgc                                                                                                                 <  8:87724/149‑1 (MQ=255)
                           aCTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgcgtgcg                                                                                                                <  4:243510/149‑1 (MQ=255)
                            ctctAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGGCTCGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgcgtgc                                                                                                                 <  1:80900/149‑1 (MQ=255)
                             tctAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGgt                                                                                                              >  6:163918/1‑149 (MQ=255)
                               tAAGTACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCTGTGACAGCGTGGATGTGCTCGCGAAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGgtgt                                                                                                            <  7:88517/149‑1 (MQ=255)
                               tAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGccccc                                                                                                                              <  8:189066/131‑1 (MQ=255)
                               tAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGgtgt                                                                                                            >  4:1234/1‑149 (MQ=255)
                                aaGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgc                                                                                                                     >  1:68183/1‑139 (MQ=255)
                                aaGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCtgc                                                                                                                     <  2:68183/139‑1 (MQ=255)
                                 aGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGtt                                                                                                           <  3:108820/148‑1 (MQ=255)
                                  gAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGc                                                                                                         <  8:36185/149‑1 (MQ=255)
                                      aaaaGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGa                                                                                                     >  2:259625/1‑149 (MQ=255)
                                         aGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGAtgt                                                                                                  >  7:166162/1‑149 (MQ=255)
                                                gtgcgATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAg                                                                                                       >  1:156244/1‑137 (MQ=255)
                                                gtgcgATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAg                                                                                                       <  2:156244/137‑1 (MQ=255)
                                                      tGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGAt                                                                                     <  3:1234/149‑1 (MQ=255)
                                                       gAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATg                                                                                    >  4:188601/1‑149 (MQ=255)
                                                          cAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCAc                                                                                 >  3:105506/1‑149 (MQ=255)
                                                            aGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGc                                                                               <  5:142827/149‑1 (MQ=255)
                                                             gATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGcc                                                                              >  7:216629/1‑149 (MQ=255)
                                                                  tctACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTc                                                                         >  3:147991/1‑149 (MQ=255)
                                                                  tctACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTc                                                                         >  2:214150/1‑149 (MQ=255)
                                                                    tACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGc                                                                                                         <  4:77265/115‑1 (MQ=255)
                                                                     aCGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTcgcg                                                                      >  8:75595/1‑149 (MQ=255)
                                                                      caacaggctgCCAGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCgg                                                                     >  4:67303/6‑149 (MQ=255)
                                                                      cGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGCTGACAGCGTGGATGTGCTAGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCgg                                                                     >  8:61648/1‑149 (MQ=255)
                                                                         ctggctgCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCggtgg                                                                  <  8:216629/149‑1 (MQ=255)
                                                                         ctggctgCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCggtgg                                                                  <  8:166162/149‑1 (MQ=255)
                                                                                 tgGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAAcgt                                                          >  4:84680/1‑149 (MQ=255)
                                                                                   gCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCCGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGCTGCACGCCTGGTCGCggt                                                                    >  8:11734/1‑137 (MQ=255)
                                                                                   gCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCggt                                                                    <  7:11734/137‑1 (MQ=255)
                                                                                   gCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAAcgtgg                                                        <  4:147991/149‑1 (MQ=255)
                                                                                   gCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAAc                                                            <  7:98794/145‑1 (MQ=255)
                                                                                   gCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAAc                                                            >  8:98794/1‑145 (MQ=255)
                                                                                     aGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAAcgtggcg                                                      <  5:153781/149‑1 (MQ=255)
                                                                                           ttGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTa                                                 <  8:172918/148‑1 (MQ=255)
                                                                                              cGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTActgc                                             >  4:24265/1‑149 (MQ=255)
                                                                                                  ctaccgccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTc                                         <  6:53576/146‑1 (MQ=255)
                                                                                                    gccgccGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTcc                                        >  3:24144/1‑148 (MQ=255)
                                                                                                     ccgccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTcc                                        <  7:75595/147‑1 (MQ=255)
                                                                                                     ccgccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCg                                       <  1:214150/148‑1 (MQ=255)
                                                                                                       gccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCggtg                                                                   <  1:167892/118‑1 (MQ=255)
                                                                                                       gccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCggtg                                                                   >  2:167892/1‑118 (MQ=255)
                                                                                                       gccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGAt                                     >  2:33756/1‑148 (MQ=255)
                                                                                                       gccgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATg                                    <  1:259625/149‑1 (MQ=255)
                                                                                                         cgG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGa                                   >  3:76690/1‑148 (MQ=255)
                                                                                                                       gTGACAGCGTGGATGTGCTCGCGCAGCCCGAGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTTGTACTGCTGTCCGATGACAATGGTGCCAcc                      <  1:161717/149‑1 (MQ=255)
                                                                                                                            aGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGgcg                 >  7:52091/1‑149 (MQ=255)
                                                                                                                              cGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGgcgc                >  5:240222/1‑148 (MQ=255)
                                                                                                                                tGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAg             <  3:67303/149‑1 (MQ=255)
                                                                                                                                 ggATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAgg            <  8:52091/149‑1 (MQ=255)
                                                                                                                                   aTGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGcc          <  4:76690/149‑1 (MQ=255)
                                                                                                                                    tgtgCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATgg                             >  2:47624/1‑129 (MQ=255)
                                                                                                                                    tgtgCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATgg                             <  1:47624/129‑1 (MQ=255)
                                                                                                                                       gCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAg             <  3:198383/142‑1 (MQ=255)
                                                                                                                                       gCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAg             >  4:198383/1‑142 (MQ=255)
                                                                                                                                        gccgtgcAGCGCGGGCTGCGCGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGCCACtgc     >  2:255901/6‑149 (MQ=255)
                                                                                                                                         ccgtgcAGCGCGGGCTGCGCGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGCCACtgcc    <  1:255901/145‑1 (MQ=255)
                                                                                                                                          cgtgcAGCGCGGGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACCGCTGTCCGATGAc                                  >  3:170054/4‑118 (MQ=255)
                                                                                                                                          cgtgcAGCGCGGGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACCGCTGTCCGATGAc                                  <  4:170054/115‑1 (MQ=255)
                                                                                                                                           gcgcAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGCCACtgccgt  <  7:61648/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CAGGCGCGCTGCGCAGGTGCTGTACCGACTCTAAGAACAAAAGGTTGCGTGCGATGAACAAGATGCTCTACGGCTGGCTGCTGGCAGCCTGTTGCGGGCTGCCGCCGG‑‑CGTTCGCCGGTGACAGCGTGGATGTGCTCGCGCAGCCCGCGCTGCACGGCCCCCAGGCTCTGCGTGCGGTGTTGCAGGATGTGGCCCGGGCCGATGCACGCCTGGTCGCGGTGGGTGAACGTGGCGTGGTACTGCTGTCCGATGACAATGGTGCCACCTGGCGCCAGGCCACTGCCGT  >  NC_002947/3632476‑3632761

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: