Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,642 G→T 6.2% T165N (ACT→AAT)  chbG ← chito‑oligosaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,816,6420GT 6.2% 214.2 / 3.6 64T165N (ACT→AAT) chbGchito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (37/23);  new base T (2/2);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

GGATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTC  >  NC_000913/1816502‑1816789
                                                                                                                                            |                                                                                                                                                   
ggATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAgg                                                                                                                                             <  8:337351/149‑1 (MQ=255)
 gATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTTCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGt                                                                                                                                            <  6:311736/149‑1 (MQ=255)
 gATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGt                                                                                                                                            <  3:71959/149‑1 (MQ=255)
        ttatCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCgg                                                                                                                                     >  8:236992/1‑149 (MQ=255)
          atCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGc                                                                                                                                    >  6:371644/1‑148 (MQ=255)
                       ggATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCAt                                                                                                                      <  2:237575/149‑1 (MQ=255)
                          tGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGCCCGGCAAATCAAACGCCATCt                                                                                                                    >  1:17609/1‑148 (MQ=255)
                          tGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCt                                                                                                                    >  2:340448/1‑148 (MQ=255)
                            acacATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGa                                                                                                                  <  1:340448/148‑1 (MQ=255)
                              acaTCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGa                                                                                                               <  7:279732/149‑1 (MQ=255)
                                 tCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATcc                                                                                                            >  1:85417/1‑149 (MQ=255)
                                     cTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgccc                                                                                                        >  4:227761/1‑149 (MQ=255)
                                     cTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgccc                                                                                                        >  1:135343/1‑149 (MQ=255)
                                      tCCAGTGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGGCCGGCAAATCAAACGCCATCTGACGATccgcccg                                                                                                       >  1:306441/1‑149 (MQ=255)
                                        cAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgcccg                                                                                                       <  3:227761/147‑1 (MQ=255)
                                              aCGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCcgcaacgcaa                                                                                               <  7:236992/149‑1 (MQ=255)
                                              aCGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCcgcaacgcaa                                                                                               <  5:371644/149‑1 (MQ=255)
                                                    aCCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATACACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCt                                                                                         >  7:297235/1‑149 (MQ=255)
                                                         gATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGctct                                                                                    >  4:363461/1‑149 (MQ=255)
                                                         gATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGctct                                                                                    >  3:286219/1‑149 (MQ=255)
                                                               cTGCGTCATCCAGCACTTGCAGGCACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGATGAATCCCTGGGTAGTGCGCAGTTTTACCGGCAAATCAAACGCCATCTGACGCTCCGCCCGCAACGCAATGCCCTGCTCTGCCGc                                                                               >  7:287427/1‑148 (MQ=255)
                                                               cTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGTTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGc                                                                               >  4:120589/1‑148 (MQ=255)
                                                                       tCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  2:351473/1‑148 (MQ=255)
                                                                       tCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  3:269690/1‑148 (MQ=255)
                                                                        ccAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGc                                                                      >  8:56642/1‑148 (MQ=255)
                                                                               ttGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCg                                                              >  6:138295/1‑149 (MQ=255)
                                                                               ttGCAGGAACAGCGACTCACAGATCTCTTCACCGTAGACTGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTCGCAACAATCg                                                              >  6:248796/1‑149 (MQ=255)
                                                                                  cAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCggg                                                            <  7:56642/148‑1 (MQ=255)
                                                                                           cGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCAGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCg                                                  >  4:120891/1‑149 (MQ=255)
                                                                                               tCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGaa                                              >  2:234943/1‑149 (MQ=255)
                                                                                                 aCTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCCACGCAATGCCCTGCTCTGCCCCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACa                                            >  4:279581/1‑149 (MQ=255)
                                                                                                 aCTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACa                                            <  3:363461/149‑1 (MQ=255)
                                                                                                  cTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACat                                           <  2:85417/149‑1 (MQ=255)
                                                                                                   tGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACata                                          >  1:301071/1‑149 (MQ=255)
                                                                                                   tGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACata                                          >  2:259843/1‑149 (MQ=255)
                                                                                                    gATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACatat                                         >  3:122293/1‑149 (MQ=255)
                                                                                                     aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       <  3:78551/118‑1 (MQ=255)
                                                                                                     aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  4:78551/1‑118 (MQ=255)
                                                                                                     aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  4:83850/1‑118 (MQ=255)
                                                                                                     aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       <  3:83850/118‑1 (MQ=255)
                                                                                                     aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATg                                        >  7:133624/1‑149 (MQ=255)
                                                                                                            cACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatga                                 >  2:230412/1‑149 (MQ=255)
                                                                                                             aCCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGGCCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatga                                 <  2:306441/148‑1 (MQ=255)
                                                                                                               cGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACattatga                              <  4:122293/149‑1 (MQ=255)
                                                                                                               cGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatga                              <  5:248796/149‑1 (MQ=255)
                                                                                                                gTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgat                             <  4:269690/149‑1 (MQ=255)
                                                                                                                gTAGAATGCACTGCTGAATCACTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgat                             <  3:6316/149‑1 (MQ=255)
                                                                                                                 tAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgatg                            <  1:351473/149‑1 (MQ=255)
                                                                                                                    aaTGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGc                          >  6:166421/1‑148 (MQ=255)
                                                                                                                      tGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTAt                       <  8:287427/149‑1 (MQ=255)
                                                                                                                         aCTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCCATGCCCTGCTCTGCCGCCAAACTGGCAACAATCGGGATAATCTGCGGGAACATTTGCACATGATGATGGCTATc                      >  4:122303/1‑147 (MQ=255)
                                                                                                                         aCTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATc                      >  7:308741/1‑147 (MQ=255)
                                                                                                                             ctgAATCCCTGGGGAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGCTCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATg                <  7:305766/149‑1 (MQ=255)
                                                                                                                              tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  1:332205/1‑149 (MQ=255)
                                                                                                                              tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  5:292880/1‑149 (MQ=255)
                                                                                                                              tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  7:281653/1‑149 (MQ=255)
                                                                                                                              tgAATCCCTGGGTAGTGCGCAGGTTCTCCGGCAAATCAACCGCCATCTGACCATCCGCCCGCAACGCAAAGCCCTGCACTGACGCAAACCTGCCAACAAGCGGGAAAAACTGCGGGACCATATGCACCTGATGATGGCTATCAACAt                 >  3:35305/1‑147 (MQ=255)
                                                                                                                                aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGTTCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGt             <  2:135343/149‑1 (MQ=255)
                                                                                                                                aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               <  4:286219/147‑1 (MQ=255)
                                                                                                                                aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCg              >  3:278978/1‑148 (MQ=255)
                                                                                                                                aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCg              >  8:377873/1‑148 (MQ=255)
                                                                                                                                     cTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGtt        >  3:257127/1‑149 (MQ=255)
                                                                                                                                        ggTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGc     <  5:166421/149‑1 (MQ=255)
                                                                                                                                           aGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGACCATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTc  >  4:390720/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTC  >  NC_000913/1816502‑1816789

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: