Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 837,636 C→T 5.1% intergenic (+200/+29) ybiC → / ← ybiJ putative dehydrogenase/DUF1471 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913837,6360CT 5.1% 255.1 / 2.4 78intergenic (+200/+29)ybiC/ybiJputative dehydrogenase/DUF1471 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base T (2/2);  total (35/43)
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.83e-01

TGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGC  >  NC_000913/837488‑837778
                                                                                                                                                    |                                                                                                                                              
tGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATgc                                                                                                                                                >  5:362072/1‑149 (MQ=255)
tGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATgc                                                                                                                                                >  4:126883/1‑149 (MQ=255)
 gAGGGTTAGAGCATAAGCGTCAGTCGGCAAACAGACAGGGAAATACTTTTGCTGGAAGTAGCGTAAACGCCTGAACCGTCCAACCGTACGGATTTAGAAGGCCTGATAACACGCGTTAGAGACGCATCAGGCCTCTGTTGCCGGATgcc                                                                                                                                               >  2:69323/1‑148 (MQ=255)
  aGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATgcgg                                                                                                                                              <  3:243522/149‑1 (MQ=255)
   gggTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTATTGCCGGATgcggc                                                                                                                                             >  4:195325/1‑149 (MQ=255)
         gagCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAAc                                                                                                                                       >  8:341556/1‑149 (MQ=255)
           gCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                     >  5:168662/1‑149 (MQ=255)
             atatGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                   <  4:121069/149‑1 (MQ=255)
              tatGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                   >  5:378069/1‑148 (MQ=255)
                tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGGTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                >  1:184889/1‑149 (MQ=255)
                  cGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                              >  3:158617/1‑149 (MQ=255)
                    tCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGt                                                                                                                            <  7:85440/149‑1 (MQ=255)
                    tCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGt                                                                                                                            <  7:115190/149‑1 (MQ=255)
                        tCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTa                                                                                                                        >  2:286699/1‑149 (MQ=255)
                         cGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTAc                                                                                                                       <  5:116838/149‑1 (MQ=255)
                           gCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACg                                                                                                                      >  5:7034/1‑148 (MQ=255)
                           gCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGAAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACg                                                                                                                      >  6:239949/1‑148 (MQ=255)
                            cAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                    <  8:4740/149‑1 (MQ=255)
                               acagacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTa                                                                                                                 <  4:331103/149‑1 (MQ=255)
                                 agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                               >  8:26677/1‑133 (MQ=255)
                                 agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                               <  7:26677/133‑1 (MQ=255)
                                 agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACt                                                                                                               <  4:158617/149‑1 (MQ=255)
                                          aaTACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTa                                                                                                      <  5:239949/149‑1 (MQ=255)
                                           aTACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTAc                                                                                                     <  3:126883/149‑1 (MQ=255)
                                             aCTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCg                                                                                                   >  7:222747/1‑149 (MQ=255)
                                               ttGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCg                                                                                                   >  8:52881/1‑147 (MQ=255)
                                                 gtgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                 <  3:158351/147‑1 (MQ=255)
                                                 gtgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                 <  3:152903/147‑1 (MQ=255)
                                                       gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCg                                                                                           <  3:74584/147‑1 (MQ=255)
                                                           tAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAAt                                                                                     <  1:286699/149‑1 (MQ=255)
                                                                  aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttg                                                                               <  6:362072/148‑1 (MQ=255)
                                                                   aCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgtt                                                                             >  4:226645/1‑149 (MQ=255)
                                                                   aCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgtt                                                                             >  4:221793/1‑149 (MQ=255)
                                                                   aCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgtt                                                                             >  2:21739/1‑149 (MQ=255)
                                                                     gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgt                                                                           >  7:355844/1‑149 (MQ=255)
                                                                       cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATaa                                                                                                          <  6:35802/116‑1 (MQ=255)
                                                                       cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgttg                                                                         <  2:336029/149‑1 (MQ=255)
                                                                         gATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTaa                                                                                      <  4:9951/134‑1 (MQ=255)
                                                                         gATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTaa                                                                                      >  3:9951/1‑134 (MQ=255)
                                                                          aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                      >  8:229653/1‑149 (MQ=255)
                                                                          aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                      >  4:70007/1‑149 (MQ=255)
                                                                          aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtgg                                                                        <  8:222747/147‑1 (MQ=255)
                                                                          aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtgg                                                                        <  4:118006/147‑1 (MQ=255)
                                                                              gTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCgg                                                                  <  7:52881/149‑1 (MQ=255)
                                                                               tCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGTGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCgg                                                                  >  1:190725/1‑148 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGc                                                                                          <  4:43188/122‑1 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGc                                                                                          >  3:43188/1‑122 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGaa                                                                <  6:168662/148‑1 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGaa                                                                <  2:273158/148‑1 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                               <  3:226645/149‑1 (MQ=255)
                                                                                 cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                               <  3:221793/149‑1 (MQ=255)
                                                                                       ttCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGc                                                         >  3:204739/1‑149 (MQ=255)
                                                                                       ttCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGc                                                         <  5:240002/149‑1 (MQ=255)
                                                                                        tCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCt                                                        >  4:344300/1‑149 (MQ=255)
                                                                                        tCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCt                                                        >  3:13551/1‑149 (MQ=255)
                                                                                         cGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCt                                                        <  7:341556/148‑1 (MQ=255)
                                                                                           gATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTATTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTa                                                     <  3:195325/149‑1 (MQ=255)
                                                                                                 gTAGGCCTGATAAGACGCGGTAGCTTCGCATCAGGCATCTGTTGCCGGATGTGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCAc                                                <  2:190725/148‑1 (MQ=255)
                                                                                                 gTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACt                                               >  2:229383/1‑149 (MQ=255)
                                                                                                  tAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTg                                              <  6:7034/149‑1 (MQ=255)
                                                                                                      ccTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                  >  2:388247/1‑61 (MQ=255)
                                                                                                      ccTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                  <  1:388247/61‑1 (MQ=255)
                                                                                                       cTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCAcc                                         >  7:177782/1‑149 (MQ=255)
                                                                                                        tGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCg                                        >  7:342527/1‑149 (MQ=255)
                                                                                                          aTAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTAt                                                    <  5:196413/135‑1 (MQ=255)
                                                                                                          aTAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTAt                                                    >  6:196413/1‑135 (MQ=255)
                                                                                                           tAAGTCGCGGTAGCGTCGCATCTGGCATCTGTTGCCGGATGCGGCGTTAACGCCTTATCCGTCCTTCGAATTTCTTATAATTTACCGCAGTTCCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                     <  1:69323/149‑1 (MQ=255)
                                                                                                                 gcgGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCcgcg                               >  1:332072/1‑149 (MQ=255)
                                                                                                                            gCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGc                    <  2:332072/149‑1 (MQ=255)
                                                                                                                             cATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtgg                                                                     <  3:17444/99‑1 (MQ=255)
                                                                                                                             cATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtgg                                                                     >  4:17444/1‑99 (MQ=255)
                                                                                                                             cATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGcc                   <  6:378069/149‑1 (MQ=255)
                                                                                                                                  ggCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGttt              <  8:342527/149‑1 (MQ=255)
                                                                                                                                   gCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTc             >  4:101407/1‑149 (MQ=255)
                                                                                                                                      tCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGCAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCGTTCTCAGCCAGTTTCGCt          >  5:122647/1‑149 (MQ=255)
                                                                                                                                        tGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAAtt                                                                                                       <  2:296202/54‑1 (MQ=255)
                                                                                                                                          ttGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTc        >  3:170042/1‑147 (MQ=255)
                                                                                                                                          ttGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTc        >  4:355301/1‑147 (MQ=255)
                                                                                                                                          ttGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCa      <  8:84513/149‑1 (MQ=255)
                                                                                                                                             ccGGATGTGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCTCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGAt                                                           >  7:347519/1‑93 (MQ=255)
                                                                                                                                             ccGGATGTGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCTCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGAt                                                           <  8:347519/93‑1 (MQ=255)
                                                                                                                                              cGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGc  <  8:250345/149‑1 (MQ=255)
                                                                                                                                              cGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGc  <  1:21739/149‑1 (MQ=255)
                                                                                                                                                   gcggcgTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCa                      <  8:127787/124‑1 (MQ=255)
                                                                                                                                                   gcggcgTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCa                      >  7:127787/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGC  >  NC_000913/837488‑837778

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: