Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 807,261 G→A 5.1% S7S (TCC→TCT)  ybhC ← acyl‑CoA thioesterase, lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,2610GA 5.1% 243.8 / 2.3 79S7S (TCC→TCT) ybhCacyl‑CoA thioesterase, lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (43/30);  new base A (2/2);  total (46/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 = 6.63e-01

GCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGAT  >  NC_000913/807118‑807407
                                                                                                                                               |                                                                                                                                                  
gCATCGAAATTCTGCGCGTCTTTTGCCGACAGAATCGGGCGCGAAGAGGTACCACGCGCGGTTTGATCTGAAGGACGTTGATCGGGCGGGGTTGAGCTAGAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAAc                                                                                                                                               <  2:370172/149‑1 (MQ=255)
     gAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGACAATGCCAGCGCCAGACGGGAAACTGaa                                                                                                                                           >  7:344788/1‑148 (MQ=255)
       aaTTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAAt                                                                                                                                        <  3:303201/149‑1 (MQ=255)
                    ttttGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGtt                                                                                                                                   >  1:132153/1‑141 (MQ=255)
                    ttttGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGtt                                                                                                                                   <  2:132153/141‑1 (MQ=255)
                         ccGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGTGAAACTGAAAATGTGTTCACAGGTTGCTcc                                                                                                                      >  5:320084/1‑149 (MQ=255)
                         ccGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTcc                                                                                                                      >  3:141661/1‑149 (MQ=255)
                           gACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCgg                                                                                                                    >  8:165604/1‑149 (MQ=255)
                           gACAGAATAGGGCGAGAAGAGGTACCAGGAGCGGTTTGATAAGAAGGCCGTTGATCGGGCGGGGTTGAGCTACAGGAGGTCAGCGTCACGCCAAAAGCCAATTCCAGCGCCAGAAGGGAAACTGAAAATGAGGTCACAGGATGCTCCgg                                                                                                                    >  4:294925/1‑149 (MQ=255)
                            aCAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCggg                                                                                                                   >  7:151903/1‑149 (MQ=255)
                            aCAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCggg                                                                                                                   >  4:195186/1‑149 (MQ=255)
                            aCAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCggg                                                                                                                   >  3:150577/1‑149 (MQ=255)
                            aCAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGAGAAACTGAAAATGTGTTCACAGGTTGCTCCggg                                                                                                                   >  7:37378/1‑149 (MQ=255)
                            aCAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGAGAAACTGAAAATGTGTTCACAGGTTGCTCCggg                                                                                                                   >  7:30764/1‑149 (MQ=255)
                                aaTCGGGCGAGAATAGGTACCAGGCGCGGTTTGAGCAGTATGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGGCAGCGTCTCGCCAAAAGCCAATGCGAGCGCCAGACGGGAACCTGAAACTGTGGTCACAGGTTGCTCCGGGCTAt                                                                                                               <  2:404010/149‑1 (MQ=255)
                                aaTCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTAt                                                                                                               >  6:132426/1‑149 (MQ=255)
                                 aTCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATg                                                                                                              <  4:141661/149‑1 (MQ=255)
                                  tCGGGCGAGAAGAGGTACCAGGCGNGGTNTGATCAGAAGGACGTTGAGCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGa                                                                                                             <  6:320084/149‑1 (MQ=255)
                                     ggCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAAt                                                                                                          <  4:150577/149‑1 (MQ=255)
                                      gCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATa                                                                                                         <  3:195186/149‑1 (MQ=255)
                                          gaagaGGTACCAGGCGCGGTTTGAGCAGAAGGACGTTGATCGGGCGGGGTGGAGCTACAGGCGGTCAGCGTCACGCAAAAAGCCAATGCCAGCTACAGACGTGAAACTTAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGaaa                                                                                                     <  3:294925/149‑1 (MQ=255)
                                               ggTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCCCGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGa                                                                                                >  2:47320/1‑149 (MQ=255)
                                                gTACCAGGCGCGGTTTGATCAGAAGGACGTAGATCGGGCGGGGTAGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTCAAAATGTGTTCACAGCTTGCTCCGGTCTATGAAATAGAAAAATGaa                                                                                               >  2:49427/1‑149 (MQ=255)
                                                   ccAGGCGCGGTTTGATCAGAAGGACTTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATcc                                                                                            >  3:187960/1‑149 (MQ=255)
                                                    cAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCg                                                                                           >  5:1450/1‑149 (MQ=255)
                                                      ggCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACTCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGtt                                                                                         >  3:58385/1‑149 (MQ=255)
                                                      ggCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGtt                                                                                         >  6:68073/1‑149 (MQ=255)
                                                             tttGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCt                                                                                  >  1:237997/1‑149 (MQ=255)
                                                                  tCAGAAGGACGTTGATCGGGCGGGCTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGGTCACAGGTTGCTCGGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCttt                                                                             >  8:214073/1‑149 (MQ=255)
                                                                   cAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGACTCCGTTGAAGCCTGCtttt                                                                            >  2:410767/1‑149 (MQ=255)
                                                                      acggACGTTGATCGGGCGTGGTGGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCAAGACGAGAAACTGAAAATGTTTTCACAGGTTGCTCCGGGCTATGAGATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTta                                                                         <  8:37378/147‑1 (MQ=255)
                                                                      aaGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTta                                                                         <  5:132426/149‑1 (MQ=255)
                                                                      aaGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGAGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGATATAGAAAAATGAATCCGGTGAAGCCTGCTTTTTta                                                                         <  8:30764/149‑1 (MQ=255)
                                                                         gACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATAc                                                                      <  8:59212/149‑1 (MQ=255)
                                                                             ttGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAc                                                                  <  6:1450/149‑1 (MQ=255)
                                                                             ttGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAc                                                                  <  2:132277/149‑1 (MQ=255)
                                                                                aTCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTTAAAATGGGTTCACAGGTAGCTCCGGGCTATGAAATAGAAAAATGAATCCGATGAAGCCTGCTTTTTTATACTAACtt                                                                >  6:125347/1‑148 (MQ=255)
                                                                                 tCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGa                                                              <  1:47320/149‑1 (MQ=255)
                                                                                   ggggggggTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGAGGGGAAACTGTAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCGTGCTTTTTTATACTAACTTGAg                                                             <  1:49427/148‑1 (MQ=255)
                                                                                     gcgggGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTATTTTATACTAACTTGAgag                                                           >  6:16069/1‑146 (MQ=255)
                                                                                         tgTTTAGCTACAGGCGTTCATCGTCACGCCAAAAGCCAATGCTTGCGCCAGACTGGAAACTTAAAATGTGTTCACAGGTTTCTCCGGGCTATGAAATAGAAAAATGAATCCTTTGAAGCCTGCTTTTTTATACTAACTTGATCGAAACg                                                      <  5:16069/148‑1 (MQ=255)
                                                                                          gTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACgg                                                     >  4:219213/1‑149 (MQ=255)
                                                                                             gAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACggg                                                    >  3:306675/1‑147 (MQ=255)
                                                                                              aGCTACAGGCGGTCAGCGTCACGCCAAAGGCCAATGCCAGCTCCAGACGGGTAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGaa                                                  <  1:281659/148‑1 (MQ=255)
                                                                                               gCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACgg                                                     >  1:234314/1‑144 (MQ=255)
                                                                                               gCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACgg                                                     <  2:234314/144‑1 (MQ=255)
                                                                                               gCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAg                                                 >  7:14991/1‑148 (MQ=255)
                                                                                               gCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAg                                                 <  8:14991/148‑1 (MQ=255)
                                                                                                 tACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGt                                               >  4:281194/1‑148 (MQ=255)
                                                                                                 tACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGt                                               >  8:331731/1‑148 (MQ=255)
                                                                                                 tACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAAGAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGt                                               >  5:162260/1‑148 (MQ=255)
                                                                                                    aGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTaaa                                            >  7:22713/1‑148 (MQ=255)
                                                                                                     ggCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTaaaaa                                          <  4:303847/149‑1 (MQ=255)
                                                                                                      gCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTaaa                                            >  5:123720/1‑146 (MQ=255)
                                                                                                        gggcAGCGTCACGCCAAAAGCCCATGCCAGCGCCACACGGCACCCTGAAAATGTGCTCACAGGGTGCTCCGGGCTATGAAAAAGAAAAAAGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGCGAACGTAAAAacac                                       <  6:162260/146‑1 (MQ=255)
                                                                                                         gTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACa                                      <  7:165604/149‑1 (MQ=255)
                                                                                                          tCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACa                                      >  7:121344/1‑148 (MQ=255)
                                                                                                             gCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATCAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACaaaaa                                  <  7:331731/149‑1 (MQ=255)
                                                                                                              cGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAg                                 >  7:378204/1‑149 (MQ=255)
                                                                                                                tCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAgtt                               <  6:123720/149‑1 (MQ=255)
                                                                                                                tCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAgtt                               >  5:163850/1‑149 (MQ=255)
                                                                                                                tCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAgtt                               >  7:227581/1‑149 (MQ=255)
                                                                                                                tCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTCAAGCCTGCATTTTTATACTAACTTGAGAGAAACGGGAAGGTAAAAAGACAAAAAgtt                               >  6:93125/1‑149 (MQ=255)
                                                                                                                  aCGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTtata                                                                       >  2:424411/1‑107 (MQ=255)
                                                                                                                  aCGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTtata                                                                       <  1:424411/107‑1 (MQ=255)
                                                                                                                        aaaGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGa                                                              >  8:248416/1‑110 (MQ=255)
                                                                                                                        aaaGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGa                                                              <  7:248416/110‑1 (MQ=255)
                                                                                                                          aGCCAATGCCACCGCCCTCCGGGTAACAGCAAATGTGTTCACCGGTTTCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATa                     <  8:7476/149‑1 (MQ=255)
                                                                                                                              aaTGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCCCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCtt                 >  1:269569/1‑149 (MQ=255)
                                                                                                                                  ccAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg             >  8:27337/1‑149 (MQ=255)
                                                                                                                                    aGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGt            >  7:154131/1‑148 (MQ=255)
                                                                                                                                    aGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGt            >  2:92115/1‑148 (MQ=255)
                                                                                                                                     gcgcCAGACGGGAAACTGAAAATGTGTTCCCAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAATCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAACTGa          <  8:344788/149‑1 (MQ=255)
                                                                                                                                        ccAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAc       <  8:151903/149‑1 (MQ=255)
                                                                                                                                          aGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAc       >  7:278612/1‑147 (MQ=255)
                                                                                                                                          aGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAc       >  4:157822/1‑147 (MQ=255)
                                                                                                                                          aGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCa     <  3:219213/149‑1 (MQ=255)
                                                                                                                                            aCGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAc       >  5:73882/1‑145 (MQ=255)
                                                                                                                                             cGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGAt  <  6:163850/149‑1 (MQ=255)
                                                                                                                                               ggAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGtttt                          >  6:252881/1‑123 (MQ=255)
                                                                                                                                               ggAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGtttt                          <  5:252881/123‑1 (MQ=255)
                                                                                                                                               ggAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGa   <  6:372392/146‑1 (MQ=255)
                                                                                                                                               ggAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGa   >  5:372392/1‑146 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGAT  >  NC_000913/807118‑807407

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