Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,653,241 T→C 8.6% K225K (AAA→AAG rspB ← putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,653,2410TC 8.6% 199.1 / 7.5 69K225K (AAA→AAGrspBputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (3/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAA  >  NC_000913/1653096‑1653386
                                                                                                                                                 |                                                                                                                                                 
cGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCtt                                                                                                                                                  <  1:35736/147‑1 (MQ=255)
   tAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTc                                                                                                                                              <  3:304952/148‑1 (MQ=255)
     aTTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTACTGTCGGCTTGATGCCTTTTTCAg                                                                                                                                            <  1:429754/148‑1 (MQ=255)
       tCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGt                                                                                                                                           >  7:107833/1‑147 (MQ=255)
       tCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGt                                                                                                                                           >  7:208970/1‑147 (MQ=255)
          ttGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAATGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCt                                                                                                                                                   <  1:301912/136‑1 (MQ=255)
          ttGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAATGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCt                                                                                                                                                   >  2:301912/1‑136 (MQ=255)
            gCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGaaa                                                                                                                                       <  2:111262/146‑1 (MQ=255)
                aaTCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCg                                                                                                                                 >  6:218081/1‑148 (MQ=255)
                aaTCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCg                                                                                                                                 <  8:107833/148‑1 (MQ=255)
                aaTCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGTGAAAAtct                                                                                                                                   <  8:208970/146‑1 (MQ=255)
                 aTCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTc                                                                                                                                              <  3:17053/134‑1 (MQ=255)
                 aTCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTc                                                                                                                                              >  4:17053/1‑134 (MQ=255)
                      ttCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCg                                                                                                                                 <  7:138820/142‑1 (MQ=255)
                      ttCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCg                                                                                                                                 >  8:138820/1‑142 (MQ=255)
                       tCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCg                                                                                                                          >  4:343237/1‑148 (MQ=255)
                         aGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGgt                                                                                                                        >  7:75677/1‑148 (MQ=255)
                            cGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTc                                                                                                                     <  6:265512/148‑1 (MQ=255)
                                      gagaACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGtta                                                                                                            <  7:156426/147‑1 (MQ=255)
                                        gaACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGttatta                                                                                                         <  5:218081/148‑1 (MQ=255)
                                                tCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGcc                                                                                                   >  6:203488/1‑146 (MQ=255)
                                                   atacaatacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTc                                                                                                                                              >  2:4367/1‑100 (MQ=255)
                                                   atacaatacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTc                                                                                                                                              <  1:4367/100‑1 (MQ=255)
                                                     acaatacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAg                                                                                            >  3:73283/1‑148 (MQ=255)
                                                       aatacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGcc                                                                                          >  7:320737/1‑148 (MQ=255)
                                                         tacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGtt                                                                                                             <  3:325806/127‑1 (MQ=255)
                                                         tacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCAGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGtt                                                                                                             >  4:325806/1‑127 (MQ=255)
                                                           cGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGc                                                                                      <  8:75677/148‑1 (MQ=255)
                                                            gTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGct                                                                                     >  8:225488/1‑148 (MQ=255)
                                                                      agaAGCCAGCGTTACGGCCTCTTTCAGGATAGACGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGctctc                                                                                  >  7:428347/1‑141 (MQ=255)
                                                                      agaAGCCAGCGTTACGGCCTCTTTCAGGATAGACGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGctctc                                                                                  <  8:428347/141‑1 (MQ=255)
                                                                          gCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCttttt                                                                       <  5:101628/148‑1 (MQ=255)
                                                                            cAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTc                                                                      >  1:70514/1‑147 (MQ=255)
                                                                            cAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCACGAGCTCTTTAGCTTTTTc                                                                      >  4:224845/1‑147 (MQ=255)
                                                                                    cGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGtt                                                              <  8:232824/147‑1 (MQ=255)
                                                                                    cGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             <  8:320737/148‑1 (MQ=255)
                                                                                      gCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtcat                                                           <  4:73283/148‑1 (MQ=255)
                                                                                        ctctTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtcatca                                                         <  7:225488/148‑1 (MQ=255)
                                                                                            ttCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGc                                                     >  8:353770/1‑148 (MQ=255)
                                                                                            ttCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGc                                                     >  8:195336/1‑148 (MQ=255)
                                                                                            ttCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGc                                                     >  3:282797/1‑148 (MQ=255)
                                                                                             tCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGc                                                     <  3:343237/147‑1 (MQ=255)
                                                                                                ggATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             <  5:349136/136‑1 (MQ=255)
                                                                                                 gATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             <  1:1649/135‑1 (MQ=255)
                                                                                                 gATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             >  2:1649/1‑135 (MQ=255)
                                                                                                 gATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCCCTCTCTTTCTCTTTTTCCAGTCGttaa                                                            >  6:349136/1‑134 (MQ=255)
                                                                                                      aaGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCACCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTAGTTCATCAATGCGATCGGcaa                                            >  6:231640/1‑147 (MQ=255)
                                                                                                       aGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGcaac                                           >  3:122585/1‑147 (MQ=255)
                                                                                                          aTGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGTGAAAATCTCGCCaa                                                                                                                             >  7:375571/1‑62 (MQ=255)
                                                                                                          aTGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGTGAAAATCTCGCCaa                                                                                                                             >  7:369267/1‑62 (MQ=255)
                                                                                                          aTGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGTGAAAATCTCGCCaa                                                                                                                             <  8:369267/62‑1 (MQ=255)
                                                                                                          aTGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCCTTTTCAGTGAAAATCTCGCCaa                                                                                                                             <  8:375571/62‑1 (MQ=255)
                                                                                                                   cGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACAttttt                               >  7:418000/1‑147 (MQ=255)
                                                                                                                        cGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             <  4:387021/112‑1 (MQ=255)
                                                                                                                        cGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTtc                                                             >  3:387021/1‑112 (MQ=255)
                                                                                                                         gATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGtt                         <  8:418000/147‑1 (MQ=255)
                                                                                                                          aTAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAAGCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTGTCGCTTTTTACAGGCGTTCATCAATGCGATCGGCAACAATCACATTGTTAACGtga                        >  5:133073/1‑145 (MQ=255)
                                                                                                                            aattaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGtt                                                                                                             <  3:221456/60‑1 (MQ=255)
                                                                                                                            aattaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGtt                                                                                                             >  4:221456/1‑60 (MQ=255)
                                                                                                                            aattaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTta                        >  6:204423/1‑145 (MQ=255)
                                                                                                                             attaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAg                     <  2:70514/147‑1 (MQ=255)
                                                                                                                             attaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAg                     <  5:204423/147‑1 (MQ=255)
                                                                                                                             attaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAg                     <  4:282797/147‑1 (MQ=255)
                                                                                                                              ttaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGAc                   >  8:24631/1‑148 (MQ=255)
                                                                                                                                    tCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCtt              <  7:195336/147‑1 (MQ=255)
                                                                                                                                        cTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTCTTATTAATCGCCCAGTCCGCCCCGGTCTCATTCGCTTATTCCAGTCGTTCATCAAAGCCATAGGCAACAATCACCTTTTTAACCTTATAGACCCCTTTTaa          >  2:248445/1‑147 (MQ=255)
                                                                                                                                             tGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATAcc      >  3:84359/1‑146 (MQ=255)
                                                                                                                                             tGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATAcc      <  4:84359/146‑1 (MQ=255)
                                                                                                                                              gCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAg                     <  6:27216/130‑1 (MQ=255)
                                                                                                                                              gCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAAACGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAg                     >  5:27216/1‑130 (MQ=255)
                                                                                                                                               ccTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGaa  >  1:292190/1‑148 (MQ=255)
                                                                                                                                               ccTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGa   <  7:24631/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAA  >  NC_000913/1653096‑1653386

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