Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,654,389 T→C 5.2% K794E (AAG→GAG)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,654,3890TC 5.2% 230.9 / 2.3 77K794E (AAG→GAG) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base C (2/2);  total (41/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 = 9.65e-01

CATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGCGGGTCGC  >  NC_002947/1654242‑1654527
                                                                                                                                                   |                                                                                                                                          
catgGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCtt                                                                                                                                            <  2:88704/148‑1 (MQ=255)
 atgGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTg                                                                                                                                           <  7:309454/148‑1 (MQ=255)
 atgGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTg                                                                                                                                           <  1:9722/148‑1 (MQ=255)
 atgGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTg                                                                                                                                           <  4:257105/148‑1 (MQ=255)
 atgGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTg                                                                                                                                           <  8:297437/148‑1 (MQ=255)
           tGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATc                                                                                                                                >  1:50482/1‑149 (MQ=255)
             aGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATc                                                                                                                                <  4:358108/147‑1 (MQ=255)
                  cccGGTTGTCTATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGTCGTGCGCGGCATCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGc                                                                                                                         <  5:363453/147‑1 (MQ=255)
                   aCGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGc                                                                                                                             >  6:402377/1‑144 (MQ=255)
                   aCGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGc                                                                                                                             <  5:402377/144‑1 (MQ=255)
                       ttGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTc                                                                                                                    >  8:353945/1‑149 (MQ=255)
                         gTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTcg                                                                                                                   <  7:71407/148‑1 (MQ=255)
                            gATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCgg                                                                                                                <  1:366539/148‑1 (MQ=255)
                             aTGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGtt                                                                                                              >  7:107168/1‑149 (MQ=255)
                               gCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTtgt                                                                                                            >  7:92034/1‑149 (MQ=255)
                                 gTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGa                                                                                                          >  5:285772/1‑149 (MQ=255)
                                  tttGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTtgtg                                                                                                           >  8:48989/1‑147 (MQ=255)
                                  tttGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTtgtg                                                                                                           <  7:48989/147‑1 (MQ=255)
                                        ggTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCtt                                                                                                   >  1:99041/1‑149 (MQ=255)
                                            aCCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCtt                                                                                                                           >  6:17508/1‑121 (MQ=255)
                                            aCCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCtt                                                                                                                           <  5:17508/121‑1 (MQ=255)
                                              cTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAAc                                                                                              >  4:370863/1‑148 (MQ=255)
                                               tGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCgg                                                                                                                >  3:374023/1‑129 (MQ=255)
                                               tGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCgg                                                                                                                <  4:374023/129‑1 (MQ=255)
                                               tGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAAc                                                                                              >  8:236802/1‑147 (MQ=255)
                                                    ggTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGc                                                                                       >  6:237885/1‑149 (MQ=255)
                                                      tgtTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGcc                                                                                      >  5:169685/1‑148 (MQ=255)
                                                          gCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGc                                                                                 >  6:413004/1‑149 (MQ=255)
                                                           ccGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGc                                                                                 <  3:370863/148‑1 (MQ=255)
                                                           ccGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCg                                                                                <  3:187147/149‑1 (MQ=255)
                                                           ccGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCg                                                                                >  6:84184/1‑149 (MQ=255)
                                                             ggCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGt                                                                              >  6:351337/1‑149 (MQ=255)
                                                             ggCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGt                                                                              <  6:285772/149‑1 (MQ=255)
                                                               ccTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCg                                                                            <  7:353945/149‑1 (MQ=255)
                                                                 tCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGAt                                                                          >  6:168086/1‑149 (MQ=255)
                                                                   gACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGc                                                                        >  6:271798/1‑149 (MQ=255)
                                                                         cgcgGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCg                                                                   <  2:371607/148‑1 (MQ=255)
                                                                          gcgGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGt                                                                 <  5:237885/149‑1 (MQ=255)
                                                                           cgGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCgg                                                                  <  4:408409/147‑1 (MQ=255)
                                                                           cgGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCgg                                                                  >  3:408409/1‑147 (MQ=255)
                                                                                tCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGt                                                           >  4:39394/1‑149 (MQ=255)
                                                                                   gTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCaa                                                        <  7:236802/149‑1 (MQ=255)
                                                                                      aGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAAccc                                                     <  2:99041/149‑1 (MQ=255)
                                                                                       gCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCg                                                    >  8:319137/1‑149 (MQ=255)
                                                                                             cATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGccacca                                              >  7:212861/1‑149 (MQ=255)
                                                                                             atgcACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGccacca                                              >  3:61627/4‑149 (MQ=255)
                                                                                              aTCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCcaccac                                             >  4:36182/1‑149 (MQ=255)
                                                                                                   tGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTc                                         >  1:144363/1‑148 (MQ=255)
                                                                                                     cgcAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACaa                                                                                            <  7:21053/95‑1 (MQ=255)
                                                                                                     cgcAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACaa                                                                                            >  8:21053/1‑95 (MQ=255)
                                                                                                     cgcAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCggt                                      <  5:351337/149‑1 (MQ=255)
                                                                                                      gcAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGAAAGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCggtg                                     >  1:318959/1‑149 (MQ=255)
                                                                                                        aaTTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATgg                                                                                                       >  6:222023/1‑81 (MQ=255)
                                                                                                        aaTTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATgg                                                                                                       <  5:222023/81‑1 (MQ=255)
                                                                                                             cAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACttt                                          >  8:327271/1‑137 (MQ=255)
                                                                                                             cAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACttt                                          <  7:327271/137‑1 (MQ=255)
                                                                                                              aGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCtggt                             >  3:39437/1‑149 (MQ=255)
                                                                                                              aGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCtggt                             <  5:84184/149‑1 (MQ=255)
                                                                                                              aGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCtggt                             <  5:413004/149‑1 (MQ=255)
                                                                                                              aGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCtggt                             <  8:212861/149‑1 (MQ=255)
                                                                                                               gATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCtggtg                            <  5:168086/149‑1 (MQ=255)
                                                                                                                     cgcTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGcc                                                  >  6:12168/1‑121 (MQ=255)
                                                                                                                     cgcTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGcc                                                  <  5:12168/121‑1 (MQ=255)
                                                                                                                                aaCTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGtt           >  4:293790/1‑149 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  >  3:114104/1‑124 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  <  4:114104/124‑1 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGc         >  4:302355/1‑149 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGc         >  6:24266/1‑149 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTCGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  <  8:203096/124‑1 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTCGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  <  8:209011/124‑1 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTCGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  >  7:209011/1‑124 (MQ=255)
                                                                                                                                  cTCGCTGCGCTGCACCTCGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGgtggtg                                  >  7:203096/1‑124 (MQ=255)
                                                                                                                                   tCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAg                                                                                           >  4:349877/1‑66 (MQ=255)
                                                                                                                                   tCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAg                                                                                           <  3:349877/66‑1 (MQ=255)
                                                                                                                                     gctgcgctgcACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGCgg       <  4:39437/148‑1 (MQ=255)
                                                                                                                                      ctgcgctgcACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACgg                >  8:265103/1‑138 (MQ=255)
                                                                                                                                      ctgcgctgcACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACgg                <  7:265103/138‑1 (MQ=255)
                                                                                                                                         cgctgcACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGCGGGTCGc  >  7:52945/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCACCTGCAGGGTGTTGCCGGCCTCGACGTGCGCGGCCTCGGTGAGCAGTGCATCACTGCGCAATTTCAGATCGCCGCTGGCCTGCAACTCGCTGCGCTGCACCTTGCTGTCGCGATCGGCTTGCACCTCGCGGTTGTGATGGCCTTCGAACAAGCTGCCCAGGCGGTCGATGCCACCGGTTACTGTCAACCCGCCACCACTTTCGGTGGTGCTGGTGGCTACGTCACGGTCGTTGCGGGTCGC  >  NC_002947/1654242‑1654527

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