Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,064 A→T 7.0% intergenic (‑158/+25) PP_0537 ← / ← ppa transcriptional regulator/inorganic pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947625,0640AT 7.0% 246.7 / 5.2 82intergenic (‑158/+25)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base T (4/2);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGC  >  NC_002947/624920‑625209
                                                                                                                                                  |                                                                                                                                                 
gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCg                                                                                                                                               >  3:304014/1‑149 (MQ=255)
gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCg                                                                                                                                               >  6:146016/1‑149 (MQ=255)
 cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGTCATGAAAAACCCCGCCAAAGCgg                                                                                                                                              <  2:369200/149‑1 (MQ=255)
       gTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTc                                                                                                                                        >  2:106225/1‑149 (MQ=255)
           ccAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTGAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCt                                                                                                                                     >  1:287568/1‑149 (MQ=255)
           ccAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGTGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTg                                                                                                                                    <  1:81917/149‑1 (MQ=255)
           ccAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGTGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTg                                                                                                                                    <  1:86980/149‑1 (MQ=255)
           ccAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttcat                                                                                                                                    <  3:142224/149‑6 (MQ=255)
           ccAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCt                                                                                                                                     <  5:146016/148‑1 (MQ=255)
            cAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTg                                                                                                                                    >  1:195705/1‑148 (MQ=255)
              gTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTccc                                                                                                                                      >  7:24585/1‑144 (MQ=255)
              gTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTccc                                                                                                                                      <  8:24585/144‑1 (MQ=255)
                  ggAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCg                                                                                                                              >  3:270968/1‑148 (MQ=255)
                     aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAAt                                                                                                                          >  4:175783/1‑149 (MQ=255)
                        gAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttcatg                                                                                                                                   <  5:193539/137‑7 (MQ=255)
                        gAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttcatg                                                                                                                                   >  6:193539/1‑131 (MQ=255)
                         aaCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAg                                                                                                                       <  4:304014/148‑1 (MQ=255)
                          aCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGcc                                                                                                                     <  2:287568/149‑1 (MQ=255)
                             gggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTGCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCt                                                                                                                   <  1:106225/149‑1 (MQ=255)
                                tgtTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                               <  6:128352/149‑1 (MQ=255)
                                    gACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGc                                                                                                            <  2:188808/148‑1 (MQ=255)
                                    gACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCa                                                                                                           >  7:21555/1‑149 (MQ=255)
                                                 cGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTgg                                                                                              <  7:265076/149‑1 (MQ=255)
                                                  gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGt                                                                                             >  7:37961/1‑149 (MQ=255)
                                                   aaaaGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAg                                                                                                                       <  7:243135/122‑1 (MQ=255)
                                                   aaaaGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAg                                                                                                                       >  8:243135/1‑122 (MQ=255)
                                                      agagGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATc                                                                                         <  8:21555/149‑1 (MQ=255)
                                                         ggCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATcgcg                                                                                      >  5:321469/1‑149 (MQ=255)
                                                            tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCa                                                                                   >  5:355257/1‑149 (MQ=255)
                                                             gTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                  >  5:167974/1‑149 (MQ=255)
                                                                             gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                   >  3:339372/1‑149 (MQ=255)
                                                                               gtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCt                                                                 >  8:176930/1‑149 (MQ=255)
                                                                                     cTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTg                                                                                                                                    <  6:88918/76‑1 (MQ=255)
                                                                                     cTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTg                                                                                                                                    >  5:88918/1‑76 (MQ=255)
                                                                                        tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCt                                                        >  7:299326/1‑149 (MQ=255)
                                                                                         gtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                       >  8:354720/1‑149 (MQ=255)
                                                                                         gtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                       >  3:128627/1‑149 (MQ=255)
                                                                                          tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCAAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTc                                                      >  5:177818/1‑149 (MQ=255)
                                                                                           aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTc                                                      >  8:259972/1‑148 (MQ=255)
                                                                                            ggATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                     <  3:175783/148‑1 (MQ=255)
                                                                                             gATTTGTTCCCAGCAGTGAGCGCGTTCCAATCAGGCATGAAAACCCCCTCCAAAGCTGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCTGCGCCTTCCCAGCCTTCGAt                                                   <  5:288548/149‑1 (MQ=255)
                                                                                              aTTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                  <  6:167974/149‑1 (MQ=255)
                                                                                               tttGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCt                                                 <  8:37961/149‑1 (MQ=255)
                                                                                               tttGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCt                                                 >  4:321445/1‑149 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggcgg                                                                            <  8:7480/115‑1 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggcgg                                                                            >  7:7480/1‑115 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCa                                          >  6:235359/1‑149 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCa                                          >  2:255055/1‑149 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCa                                          >  8:159033/1‑149 (MQ=255)
                                                                                                       cagcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAccc                                           <  8:163411/147‑1 (MQ=255)
                                                                                                       cagcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAccc                                           >  7:163411/1‑147 (MQ=255)
                                                                                                        agcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACt                                        >  7:179247/1‑149 (MQ=255)
                                                                                                        agcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACt                                        >  5:226295/1‑149 (MQ=255)
                                                                                                          cagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTg                                      >  6:72207/1‑149 (MQ=255)
                                                                                                             tGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGccc                                   >  3:183804/1‑149 (MQ=255)
                                                                                                                gcgcGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCtt                                 >  8:152297/1‑148 (MQ=255)
                                                                                                                gcgcGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                >  5:199090/1‑149 (MQ=255)
                                                                                                                  gcgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                              >  2:233933/1‑149 (MQ=255)
                                                                                                                  gcgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                              <  4:270968/149‑1 (MQ=255)
                                                                                                                     ttCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                              >  2:112717/1‑146 (MQ=255)
                                                                                                                         cATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGAt                         <  4:339372/148‑1 (MQ=255)
                                                                                                                            tCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                     <  7:354720/149‑1 (MQ=255)
                                                                                                                             cAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                     <  7:152297/148‑1 (MQ=255)
                                                                                                                             cAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                     <  4:183804/148‑1 (MQ=255)
                                                                                                                              aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTg                    <  6:199090/148‑1 (MQ=255)
                                                                                                                               ggCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                  <  6:321469/149‑1 (MQ=255)
                                                                                                                                gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGCCCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAg                 >  5:103547/1‑149 (MQ=255)
                                                                                                                                gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAg                 <  6:355257/149‑1 (MQ=255)
                                                                                                                                 gcagggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGGGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGt                >  1:118488/7‑149 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                     <  1:148688/143‑1 (MQ=255)
                                                                                                                                   tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg              >  1:303950/1‑149 (MQ=255)
                                                                                                                                    gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                               <  2:337190/67‑1 (MQ=255)
                                                                                                                                    gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                               >  1:337190/1‑67 (MQ=255)
                                                                                                                                     ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg             >  6:135221/3‑148 (MQ=255)
                                                                                                                                     ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg             >  3:131454/3‑148 (MQ=255)
                                                                                                                                     aaaaaCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg             <  8:179247/148‑1 (MQ=255)
                                                                                                                                     aaaaaCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg             <  2:119582/148‑1 (MQ=255)
                                                                                                                                           cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCTACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa      <  5:235359/149‑1 (MQ=255)
                                                                                                                                              gCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTg   <  5:72207/149‑1 (MQ=255)
                                                                                                                                              gCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTg   <  4:128627/149‑1 (MQ=255)
                                                                                                                                               ccAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGc  <  7:176930/149‑1 (MQ=255)
                                                                                                                                                cAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTTATAGCGGCACGGGCGGCGTCCGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTATCGACATCTTTGTAGTTGGCAAAGAAGTGc  >  1:66494/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGC  >  NC_002947/624920‑625209

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: