Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,876,733 G→T 5.9% R214R (CGC→CGA)  PP_2531 ← 5‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,876,7330GT 5.9% 217.4 / 2.7 68R214R (CGC→CGA) PP_25315‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base T (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

AAAAAAGCACTCATCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCAATCACGCCA  >  NC_002947/2876586‑2876880
                                                                                                                                                   |                                                                                                                                                   
aaaaaaGCACTCATCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCCGCATGc                                                                                                                                                    >  6:131395/1‑149 (MQ=255)
 aaaaaGCACTCATCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGc                                                                                                                                                    <  6:298408/148‑1 (MQ=255)
  aaaaGCACTCATCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCACGGCCAGCTTGCGGTCCTCCAGCAAGCg                                                                                                                                                   >  8:180348/1‑148 (MQ=255)
            aTCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGa                                                                                                                                          >  4:159743/1‑147 (MQ=255)
                   gATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCg                                                                                                                                 <  7:350260/149‑1 (MQ=255)
                    aTACGGCCCAATCCTTCAAGGTCGTCACCACCGGCTCATGGGAAGGCGCGTTGCGCGCGGTGACGTTGGCAATGCGGAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCTGTGCTCCAGCATGCGTATCCGCTATAGTTTGCgg                                                                                                                                <  7:205762/149‑1 (MQ=255)
                         ccccAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAg                                                                                                                           >  8:297645/1‑149 (MQ=255)
                          cccAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTCACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAg                                                                                                                           >  7:208312/1‑148 (MQ=255)
                             aaTCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCt                                                                                                                       <  1:378032/149‑1 (MQ=255)
                                 cTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGc                                                                                                                     >  7:348295/1‑147 (MQ=255)
                                 cTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCa                                                                                                                   <  3:159743/149‑1 (MQ=255)
                                  ttCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGAGAGTTAGCGGTAAAGGTCTGCCAg                                                                                                                  >  5:12362/1‑149 (MQ=255)
                                   tCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGt                                                                                                                                      <  2:100793/128‑1 (MQ=255)
                                   tCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGACAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGt                                                                                                                                      >  1:100793/1‑128 (MQ=255)
                                    cAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTg                                                                                                                <  7:50667/149‑1 (MQ=255)
                                       ggTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCAGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGcc                                                                                                             >  4:82853/1‑149 (MQ=255)
                                            tcACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGc                                                                                                        >  6:330684/1‑149 (MQ=255)
                                            tcACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGc                                                                                                        <  5:187644/149‑1 (MQ=255)
                                             caccacCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGc                                                                                                        >  7:154793/1‑148 (MQ=255)
                                             caccacCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGc                                                                                                        >  7:172017/1‑148 (MQ=255)
                                             caccacCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGc                                                                                                        >  8:351885/1‑148 (MQ=255)
                                              accaccCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCg                                                                                                       >  4:17477/1‑148 (MQ=255)
                                                caccCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAAt                                                                                                    <  4:110745/149‑1 (MQ=255)
                                                    cGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCgg                                                                                                <  1:366017/149‑1 (MQ=255)
                                                    cGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCgg                                                                                                <  1:177822/149‑1 (MQ=255)
                                                    cGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCgg                                                                                                <  6:174418/149‑1 (MQ=255)
                                                     gCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCggg                                                                                               >  7:103795/1‑149 (MQ=255)
                                                         aTGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTGGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGt                                                                                             <  3:9920/147‑1 (MQ=255)
                                                                 gcgcGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCg                                                                                                 <  7:157736/135‑1 (MQ=255)
                                                                 gcgcGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCg                                                                                                 >  8:157736/1‑135 (MQ=255)
                                                                         cgcgcCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTc                                                                           <  4:208363/149‑1 (MQ=255)
                                                                          gcgcCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCg                                                                          >  3:196108/1‑149 (MQ=255)
                                                                           cgcCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCgt                                                                         <  7:320886/149‑1 (MQ=255)
                                                                           cgcCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCGCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCgt                                                                         <  7:326140/149‑1 (MQ=255)
                                                                                tGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTcc                                                                    <  1:155043/149‑1 (MQ=255)
                                                                                tGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCACGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTcc                                                                    <  5:330684/149‑1 (MQ=255)
                                                                                    gATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGa                                                                 >  6:88062/1‑148 (MQ=255)
                                                                                     aTGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGaaa                                                               <  8:208312/149‑1 (MQ=255)
                                                                                            tGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTcc                                                        <  7:297645/149‑1 (MQ=255)
                                                                                                 aGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGtt                                                   <  5:357672/149‑1 (MQ=255)
                                                                                                 aGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGtt                                                   <  4:196108/149‑1 (MQ=255)
                                                                                                     tCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTtgt                                                                                     >  4:66639/1‑111 (MQ=255)
                                                                                                     tCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTtgt                                                                                     <  3:66639/111‑1 (MQ=255)
                                                                                                       cGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGtt                                             <  3:82853/149‑1 (MQ=255)
                                                                                                       cGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGtt                                             <  8:103795/149‑1 (MQ=255)
                                                                                                        gCTGGTCCGTCGGGTCGTCGGGCCGCGCGCGGTCCGGGAGCTTGCGTATCCGCGATAGTTTGCGGTATAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGGGCTCGTGTTCCTGAAAGTCGTCCTGGTTGTCGTTg                                            <  8:329542/149‑1 (MQ=255)
                                                                                                        gCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTg                                            >  5:228545/1‑149 (MQ=255)
                                                                                                            gTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTgcgc                                         <  3:17477/148‑1 (MQ=255)
                                                                                                             tACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATTCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCtt                                       >  8:91606/1‑149 (MQ=255)
                                                                                                             tACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATTCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCtt                                       >  8:327745/1‑149 (MQ=255)
                                                                                                               cGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCa                                                                                  >  6:103253/1‑104 (MQ=255)
                                                                                                               cGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCa                                                                                  <  5:103253/104‑1 (MQ=255)
                                                                                                                gTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGt                                     >  2:21479/1‑148 (MQ=255)
                                                                                                                     gTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAacac                               <  7:351885/149‑1 (MQ=255)
                                                                                                                     gTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAacac                               <  1:21479/149‑1 (MQ=255)
                                                                                                                         tCGGCCAGCTTGCGGTCGTCCTGGATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCg                                               <  4:220179/129‑1 (MQ=255)
                                                                                                                         tCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCg                                               >  3:220179/1‑129 (MQ=255)
                                                                                                                                  ttGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTaaa                                  <  7:345441/133‑1 (MQ=255)
                                                                                                                                  ttGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTaaa                                  >  8:345441/1‑133 (MQ=255)
                                                                                                                                  ttGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACActctct                          <  2:187304/141‑1 (MQ=255)
                                                                                                                                  ttGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACActctct                          >  1:187304/1‑141 (MQ=255)
                                                                                                                                        tcctccAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCgg             >  1:93453/1‑148 (MQ=255)
                                                                                                                                        tcctccAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCgg             <  2:93453/148‑1 (MQ=255)
                                                                                                                                         cctccAGCATTCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCa           <  7:91606/149‑1 (MQ=255)
                                                                                                                                         cctccAGCATTCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCa           <  7:327745/149‑1 (MQ=255)
                                                                                                                                         cctccAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGc            >  7:214888/1‑148 (MQ=255)
                                                                                                                                             cAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCAATCa       >  7:89695/1‑149 (MQ=255)
                                                                                                                                                  tGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCAATCACGCCa  <  8:214888/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAAAAGCACTCATCCGGGGATACGCCCCAATCCTTCAAGGTCGTCACCACCCGCTCATGCGAAGGCGCGTTGCGCGCCGTGACGATGGCAATGCGCAGAATCCGCTGGTACGTCGGGTCCTCGGCCAGCTTGCGGTCCTCCAGCATGCGTATCCGCGATAGTTTGCGGTAAAGGTCTGCCAGTGGCCCCGGCGAATGCGGGGTGCTCTTGTGCATGCGCTCGTGTTCCTGAAAGTCGTCCAGGTTGTCGTTGCGCTTGTAAACACTCTCTGCCTCGTCGTCGGCAATCACGCCA  >  NC_002947/2876586‑2876880

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