Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,035,378 G→A 5.9% V105M (GTG→ATG)  hyaC → hydrogenase 1, b‑type cytochrome subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,035,3780GA 5.9% 212.0 / 2.4 67V105M (GTG→ATG) hyaChydrogenase 1, b‑type cytochrome subunit
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

TGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGCGAA  >  NC_000913/1035232‑1035524
                                                                                                                                                  |                                                                                                                                                  
tGGGCTACAGCAGGTGAATTCACTTCAGCGCCGGGATGGGTGTTACCGTGGTGTTGCTGAGGCGGATCTCCTGGGCTTTTGTTGGCAAGCGATACTCCCGCGCGCTGTGTATCGTGCCGGGATGGCGTAAAAGCGGGGGGCAGGGCgtg                                                                                                                                                  <  5:26874/149‑1 (MQ=255)
   gCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTgg                                                                                                                                               >  8:276447/1‑149 (MQ=255)
    cTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTTCCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGt                                                                                                                                              <  6:40522/149‑1 (MQ=255)
    cTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGt                                                                                                                                              <  8:156435/149‑1 (MQ=255)
        aTCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATg                                                                                                                                           >  5:46767/1‑148 (MQ=255)
             gTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCAGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATcc                                                                                                                                     >  4:358400/1‑149 (MQ=255)
                 aTTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGTCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTg                                                                                                                                 >  6:244575/1‑149 (MQ=255)
                 aTTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCt                                                                                                                                  <  2:136032/148‑1 (MQ=255)
                 aTTCACTTCAGCGCAGGGATGGTTTTTACCTTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTg                                                                                                                                 <  8:125830/149‑1 (MQ=255)
                  ttCACTGCACCGCCGGGATCGTTTTTACCGTGGATGTGCTGATGCGTATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGCGATGTGTTTCGTGCCGGTATGGCGTAAAAGCTTGTGGCAGGGCGTGTGGTTTGAAATCCGCTgg                                                                                                                                <  6:215710/149‑1 (MQ=255)
                     aCTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTa                                                                                                                              <  3:168587/148‑1 (MQ=255)
                     aCTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTAt                                                                                                                             >  6:258365/1‑149 (MQ=255)
                       ttCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCt                                                                                                                           >  8:278029/1‑149 (MQ=255)
                       ttCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCt                                                                                                                           >  4:114362/1‑149 (MQ=255)
                                   aTGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCaa                                                                                                                >  3:17656/1‑148 (MQ=255)
                                       tttttACCGTGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCaaaa                                                                                                              >  2:344348/1‑145 (MQ=255)
                                       tttttACCGTGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCaaaa                                                                                                              <  1:344348/145‑1 (MQ=255)
                                         tttACCGTGGTTTTGCTGATGCTGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTc                                                                                                          <  7:276447/148‑1 (MQ=255)
                                         tttACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTc                                                                                                          <  7:278029/148‑1 (MQ=255)
                                                  gTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGa                                                                                                <  8:101372/149‑1 (MQ=255)
                                                         tGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGc                                                                                         >  8:215247/1‑149 (MQ=255)
                                                           aTGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGAGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGc                                                                                         <  4:59587/147‑1 (MQ=255)
                                                             gCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCATGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGc                                                                                                   >  4:186065/1‑135 (MQ=255)
                                                             gCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCATGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGc                                                                                                   >  4:179969/1‑135 (MQ=255)
                                                             gCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCATGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGc                                                                                                   <  3:179969/135‑1 (MQ=255)
                                                             gCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCATGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGc                                                                                                   <  3:186065/135‑1 (MQ=255)
                                                               ggATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAAt                                                                                   >  3:178239/1‑149 (MQ=255)
                                                                gATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATc                                                                                  >  2:252160/1‑149 (MQ=255)
                                                                  tCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATCGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATcc                                                                                 >  8:198691/1‑148 (MQ=255)
                                                                      cTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCg                                                                            >  7:68645/1‑149 (MQ=255)
                                                                       tGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGc                                                                           >  8:366532/1‑149 (MQ=255)
                                                                       tGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGc                                                                           >  3:256956/1‑149 (MQ=255)
                                                                           cTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAg                                                                       <  4:17656/149‑1 (MQ=255)
                                                                                 ttgGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGGGTGTGGTTTGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGCGTGCAGATATAGGCCATAATCCAATCGCCAAggcg                                                                    >  6:365827/1‑146 (MQ=255)
                                                                                 ttgGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAggcggc                                                                  >  1:111157/1‑148 (MQ=255)
                                                                                    gCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATg                                                              >  3:359748/1‑149 (MQ=255)
                                                                                      aaTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGtt                                                            >  8:58242/1‑149 (MQ=255)
                                                                                       aTCGATACTCGCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTc                                                           <  6:46767/149‑1 (MQ=255)
                                                                                       aTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTTGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTc                                                           <  8:321901/149‑1 (MQ=255)
                                                                                       aTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTc                                                           >  8:83903/1‑149 (MQ=255)
                                                                                       aTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTc                                                           >  7:403029/1‑149 (MQ=255)
                                                                                        tCGATACTCCCGCTATCTGTTTATGGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGGGTGTTATGATATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGtt                                                            <  2:111240/147‑1 (MQ=255)
                                                                                          gATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCgg                                                         <  1:252160/148‑1 (MQ=255)
                                                                                            tACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGAGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTa                                                      >  2:140188/1‑149 (MQ=255)
                                                                                                 ccGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTcc                                                 >  2:20091/1‑149 (MQ=255)
                                                                                                 ccGCGAGCTGGTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGTTATTTcc                                                 <  8:403029/149‑1 (MQ=255)
                                                                                                  cgcgAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCt                                                <  4:359748/149‑1 (MQ=255)
                                                                                                           tttATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGt                                       >  4:331310/1‑149 (MQ=255)
                                                                                                             tATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGGGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTc                                      <  1:140188/148‑1 (MQ=255)
                                                                                                              aTCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGGCt                                     <  7:83903/148‑1 (MQ=255)
                                                                                                               tCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATACGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCCGTCtt                                    >  7:252617/1‑148 (MQ=255)
                                                                                                                 gTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTAt                                 >  5:243124/1‑149 (MQ=255)
                                                                                                                   gCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGa                               <  4:256956/149‑1 (MQ=255)
                                                                                                                    ccGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAAGTTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGTTATTTCCTGATGTCGGTCTTTATGat                              <  5:85510/149‑1 (MQ=255)
                                                                                                                    ccGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGa                               <  3:358400/148‑1 (MQ=255)
                                                                                                                      ggTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGatca                            <  2:111157/149‑1 (MQ=255)
                                                                                                                        tATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGatcatc                          <  3:114362/149‑1 (MQ=255)
                                                                                                                          tGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCAc                        <  5:258365/149‑1 (MQ=255)
                                                                                                                           ggCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACt                       >  7:141093/1‑149 (MQ=255)
                                                                                                                            gCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACATCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTg                      <  5:244575/149‑1 (MQ=255)
                                                                                                                             cGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGACGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCCCTgg                     >  6:277448/1‑149 (MQ=255)
                                                                                                                                      tggtggCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTcc                                                 <  7:312307/112‑1 (MQ=255)
                                                                                                                                      tggtggCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTcc                                                 >  8:312307/1‑112 (MQ=255)
                                                                                                                                            cAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAg      >  8:6851/1‑149 (MQ=255)
                                                                                                                                             aGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGc     <  7:215247/149‑1 (MQ=255)
                                                                                                                                               ggCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGCg    >  2:201983/1‑148 (MQ=255)
                                                                                                                                                gCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGTTGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGCGaa  <  1:139846/149‑1 (MQ=255)
                                                                                                                                                gCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGCGaa  <  3:250941/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGATGTCGGTCTTTATGATCATCACTGGTTTTGCGCTGTACAGCGAA  >  NC_000913/1035232‑1035524

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: