Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,192 C→G 6.2% S175C (TCC→TGC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,1920CG 6.2% 168.8 / 2.6 64S175C (TCC→TGC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base G (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAA  >  NC_000913/4262046‑4262335
                                                                                                                                                  |                                                                                                                                               
gtgtCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                                                                                                                                               >  1:149336/1‑149 (MQ=255)
  gtCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCgg                                                                                                                                              >  8:168821/1‑148 (MQ=255)
     tGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCa                                                                                                                                            >  2:273304/1‑147 (MQ=255)
          ggTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                                                                                                       >  4:101225/1‑147 (MQ=255)
          ggTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGa                                                                                                                                     <  3:357307/149‑1 (MQ=255)
           gTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAg                                                                                                                                    >  5:209035/1‑149 (MQ=255)
            tAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCTATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                                                                                                                   <  1:108750/149‑1 (MQ=255)
             aaTGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCAGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgtg                                                                                                                                  >  2:28758/1‑149 (MQ=255)
                 cgcAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgtgt                                                                                                                                 >  1:44635/1‑146 (MQ=255)
                 cgcAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                                                                                                               <  6:173799/148‑1 (MQ=255)
                     gCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGt                                                                                                                          <  6:158439/149‑1 (MQ=255)
                         gTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACTGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtc                                                                                                                        <  3:134125/147‑1 (MQ=255)
                         gTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtca                                                                                                                       <  1:273304/148‑1 (MQ=255)
                         gTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtc                                                                                                                        <  6:209035/147‑1 (MQ=255)
                            gCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtcatca                                                                                                                    <  2:149336/148‑1 (MQ=255)
                             ccGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                                                                                                                  >  3:8550/1‑149 (MQ=255)
                              cGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                                                                                                                  >  3:232026/1‑148 (MQ=255)
                                aCTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcac                                                                                                               <  4:8550/149‑1 (MQ=255)
                                aCTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcac                                                                                                               <  7:218852/149‑1 (MQ=255)
                                aCTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcac                                                                                                               <  7:224240/149‑1 (MQ=255)
                                  tGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacg                                                                                                              >  7:284257/1‑148 (MQ=255)
                                     gTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                                                                                                          >  4:186456/1‑149 (MQ=255)
                                        aaaGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTa                                                                                                        >  3:229933/1‑148 (MQ=255)
                                             gATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGCCAAAGGCGAGTGTGAGATCTTCATCATCCACGCACGTAAAGcat                                                                                                  >  6:1908/1‑147 (MQ=255)
                                              aTGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTg                                                                                                 <  3:101225/149‑1 (MQ=255)
                                                               gATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCt                                                                                              >  3:345964/1‑135 (MQ=255)
                                                               gATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCt                                                                                              <  4:345964/135‑1 (MQ=255)
                                                                aTTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                               >  5:35928/1‑149 (MQ=255)
                                                                 ttCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTTTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCt                                                                                                  <  5:222894/129‑1 (MQ=255)
                                                                 ttCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTTTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCt                                                                                                  >  6:222894/1‑129 (MQ=255)
                                                                 ttCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                               >  8:264541/1‑148 (MQ=255)
                                                                           ggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                    >  7:43406/1‑149 (MQ=255)
                                                                            gtgaAAACGCGTATTGGCATCGACTACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                   <  7:20553/149‑1 (MQ=255)
                                                                               aaaaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                <  7:168821/149‑1 (MQ=255)
                                                                                     cgTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCACTACGGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                          <  8:172266/149‑1 (MQ=255)
                                                                                              aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGa                                                 >  6:380291/1‑149 (MQ=255)
                                                                                                cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGt                                                                  <  7:205574/130‑1 (MQ=255)
                                                                                                cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGt                                                                  >  8:205574/1‑130 (MQ=255)
                                                                                                    gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATcc                                                           <  2:314715/133‑1 (MQ=255)
                                                                                                       cAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                        >  5:205020/1‑149 (MQ=255)
                                                                                                        aGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgt                                       <  2:44635/149‑1 (MQ=255)
                                                                                                        aGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgt                                       <  8:43406/149‑1 (MQ=255)
                                                                                                        aGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATC‑‑CGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                     >  1:246654/1‑149 (MQ=255)
                                                                                                          gACAGCTTTGAATTTCTCTGCGATTTCATCAATACAGTTTCCGGCAAAGGCGCGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGTCTTTCGGGGTTAATCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                     <  1:351242/149‑1 (MQ=255)
                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGTTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                               <  2:115782/123‑1 (MQ=255)
                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGGTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                               >  1:115782/1‑123 (MQ=255)
                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                     <  1:28758/149‑1 (MQ=255)
                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                     <  1:345867/149‑1 (MQ=255)
                                                                                                               cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTAt                                 >  4:76344/1‑148 (MQ=255)
                                                                                                                            tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                   >  7:89784/1‑149 (MQ=255)
                                                                                                                             gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                  >  3:353280/1‑149 (MQ=255)
                                                                                                                              cGATTTCATCAATACCGTTTGCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                    <  5:339979/114‑1 (MQ=255)
                                                                                                                              cGATTTCATCAATACCGTTTGCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                    <  5:345626/114‑1 (MQ=255)
                                                                                                                              cGATTTCATCAATACCGTTTGCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                    >  6:339979/1‑114 (MQ=255)
                                                                                                                              cGATTTCATCAATACCGTTTGCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCACGAAAGAAAACCGTGAAATCCCTCCGCt                                                    >  6:345626/1‑114 (MQ=255)
                                                                                                                              agATTTCATCAATACCGTTTCCGGCATAGGCGAGTGTGAGATGTTCATCATGCACGCACGTAAAGGCTGGATTTCGGGGTTAAGCCCGATAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                  <  5:1908/147‑1 (MQ=255)
                                                                                                                               gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt                >  7:333244/1‑149 (MQ=255)
                                                                                                                                    catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg            >  7:337832/1‑148 (MQ=255)
                                                                                                                                    catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg            >  2:305182/1‑148 (MQ=255)
                                                                                                                                     atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa          >  8:338157/1‑149 (MQ=255)
                                                                                                                                        aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc        <  4:353280/148‑1 (MQ=255)
                                                                                                                                        aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGAGTGTATCAACTGAAGCGTGACTTTCCGCATCt       >  1:67292/1‑149 (MQ=255)
                                                                                                                                            ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATC‑‑CGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa  <  2:246654/148‑1 (MQ=255)
                                                                                                                                             cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa  <  8:284257/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAA  >  NC_000913/4262046‑4262335

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: ‑