Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,917 T→C 7.0% intergenic (+315/+299) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9170TC 7.0% 188.0 / 4.4 72intergenic (+315/+299)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (24/42);  new base C (2/3);  total (27/45)
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.70e-01

TGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTC  >  NC_000913/3269775‑3270033
                                                                                                                                               |                                                                                                                    
tGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGt                                                                                                                 <  8:298038/148‑1 (MQ=255)
  cGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGa                                                                                                              <  7:376475/149‑1 (MQ=255)
  cGTGCTTTTGTAGGCAGGATAAGGAGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGa                                                                                                              <  1:324072/149‑1 (MQ=255)
    tGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACa                                                                                                            <  4:324204/149‑1 (MQ=255)
     gCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAg                                                                                                           >  2:78955/1‑149 (MQ=255)
      cTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAg                                                                                                           <  2:371007/148‑1 (MQ=255)
         ttGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGtca                                                                                                        <  5:161354/148‑1 (MQ=255)
            tAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGtcattc                                                                                                     <  1:355335/148‑1 (MQ=255)
              ggCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                                                                  <  4:262557/149‑1 (MQ=255)
              ggCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                                                                  <  8:237875/149‑1 (MQ=255)
                 aGGATAAGGCGTTCACGCCGCATCCGGCATGAATACAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGCCCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               <  6:22907/149‑1 (MQ=255)
                        ggCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTGTGGGTATAGTGTGGTGGACAGGGATTGATCTTTCTgacgc                                                                                        >  5:119867/1‑145 (MQ=25)
                           gTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAg                                                                                                           <  7:211477/127‑1 (MQ=255)
                           gTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAg                                                                                                           >  8:211477/1‑127 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               >  1:228642/1‑138 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               <  2:228642/138‑1 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           >  8:47569/1‑142 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           <  7:47569/142‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGACCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTGCTTCTTCTTTCTGCCCCTCCaa                                                                                   >  1:373250/1‑147 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  3:140027/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  3:347251/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  7:275234/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                                                                   >  5:16226/1‑147 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTC‑TTTTGTGTATAGTGTCGTGGACAGTAATTCATATTTCTGCCACTCaaa                                                                                   >  1:314442/1‑147 (MQ=21)
                                   cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTATC‑TTATGCGAATAGTGTAGTGGAAAATAATTAATATTTCTGACCCTCAAAAAGTa                                                                              >  2:154337/1‑148 (MQ=25)
                                    gcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                            <  3:305850/149‑1 (MQ=255)
                                       tCCGGCATGAATAAAGCGCACTAAACAATCTGAAACCAGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTCTATTGGAAAGTCAGAAATCTTTCAGCCACTCACAAGTTAAAAAc                                                                         >  1:77500/1‑149 (MQ=21)
                                            cATGAATAAAGCGCACTCAACAATCTGAAACCCTCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCTCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccg                                                                      <  8:248886/147‑1 (MQ=255)
                                            cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGACCCTCCAAAAAAAAAAAcccc                                                                      >  8:306983/1‑146 (MQ=255)
                                               gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTgg                                                                                                               >  8:378789/1‑103 (MQ=255)
                                               gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTgg                                                                                                               <  7:378789/103‑1 (MQ=255)
                                                aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTAATTCAGAATTCGACACCTCAAAAAGAAAAAAAccgccggcg                                                                >  7:300532/1‑146 (MQ=18)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           <  4:322938/113‑1 (MQ=255)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           >  3:322938/1‑113 (MQ=255)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                       <  7:42818/149‑1 (MQ=255)
                                                          aCTCAACTATCTGAAACCCGGCTGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                      <  8:72578/149‑1 (MQ=255)
                                                           cTCAACAATCACAAACCCGCCGTACCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  2:80269/149‑1 (MQ=11)
                                                            tctacaaTCTGAAACCCGCCGTAACGGTTTATTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGCGTCGTGGACa                                                                                                            <  6:289125/90‑1 (MQ=17)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  6:91380/1‑149 (MQ=255)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACg                                                    >  1:80269/1‑149 (MQ=255)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTAAAGTGACGAGGACACGAATTAATAAGTCTGCCCAGAAAAAAGTAAAAAACCGCAGAACCGGGtgattaa                                                     >  7:288785/1‑142 (MQ=16)
                                                              aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  5:138092/146‑1 (MQ=255)
                                                              aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                     >  6:138092/1‑146 (MQ=255)
                                                                 aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGAACGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGGCAGTCCTTCCTCTTTCTGCCCCTCCAAAAGAAAAAACACGCCAAAGGGGGTTTTTACGTaaaa                                               >  5:322470/1‑149 (MQ=17)
                                                                  aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                              <  3:68784/149‑1 (MQ=255)
                                                                  aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTaaaa                                               >  4:68784/1‑148 (MQ=255)
                                                                     tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTAGTGGACAGTCCTTAATCTGTATGCCCCTGCATAAGTAAAAAACCGCCGCATCGGGTTTATAAATAAAACAgg                                           >  5:220814/1‑149 (MQ=17)
                                                                       aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                            <  2:409858/146‑1 (MQ=255)
                                                                        aaCCCGCCGGAACTGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                        <  7:80626/149‑1 (MQ=255)
                                                                          cctgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTAGTGGAAAGTGATTGATATTTCTGCCCTTACAACAGCAAAAACCAGACGAAGAGGGTTTTTACGTATAAAAGGTGaa                                       >  7:373362/4‑148 (MQ=14)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  6:16226/149‑1 (MQ=255)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  5:232582/149‑1 (MQ=255)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTATGGGTATAGTGCCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  1:325607/149‑1 (MQ=255)
                                                                          ccagccgGAACGGTTTATTACGTACATCATGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCTGGTTTTTACGTAAAACAGGTGaaa                                      <  1:20284/146‑1 (MQ=255)
                                                                            cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                    <  5:91380/149‑1 (MQ=255)
                                                                                  aaCGGTATATTACGTACATCAGGTACAACTGACCGATAAGCCGCTTTC‑TTGTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                              <  4:324219/149‑1 (MQ=14)
                                                                                   actgTTTATTACTAACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                             <  1:362022/146‑1 (MQ=255)
                                                                                    cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                            <  5:112132/149‑1 (MQ=255)
                                                                                         taataCGTTCATCATGTAAATCTGTCCTATATTCCTCTTTC‑TTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTATGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                        <  1:154337/145‑1 (MQ=14)
                                                                                           tttcGTACATAAGTTAATACTGACCGATAAGCCGCTTTC‑TTTTTGTTATAGTGTCGTGGACAGTGATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                     <  1:284106/146‑1 (MQ=14)
                                                                                               gTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGCGTCGTGGACa                                                                                                            >  5:109731/1‑58 (MQ=255)
                                                                                               gTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGCGTCGTGGACa                                                                                                            <  6:109731/58‑1 (MQ=255)
                                                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGAGTCGTGGACAGTCATTCATCt                                                                                                 >  4:247450/1‑65 (MQ=255)
                                                                                                    tcatttAAAACTGTCCGATAAGTAGTTTTCATTTTTTGT‑TTGTGTCGTGGACAGTCTTTCTTATATCTGCCCCTCCAAAAGTACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg            <  2:77296/144‑1 (MQ=11)
                                                                                                    tcatgcaaaaCTGTCCGATAAGCCGCTTTC‑TTTAGGGTATACTGGCTTGGACACTCTTTGCTCTTTCTCCCTCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg            <  8:162478/143‑1 (MQ=11)
                                                                                                        gTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 >  3:349542/1‑76 (MQ=34)
                                                                                                        gTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  4:349542/76‑1 (MQ=34)
                                                                                                          aataattaCCGTTAAGCTGATTGC‑TTTTGTCTATTGTGTCGACGACAGGCTTAATTCTATCTCCCACTCCAAAATTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg      <  1:299671/142‑1 (MQ=2)
                                                                                                              gagACCTATTATCCGCATTC‑ATTAGGGCATAGCGTCGTGGACACTCCTTCTTCAGTCAGCCCCTGCATCAGATATAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc  <  6:392511/147‑1 (MQ=2)
                                                                                                              ctttccGATATGTCGATTTC‑AGTTGTCTCTAGTGTCGTCGACAGCGTTTTCTCTTTCTCCCACTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc  <  7:295649/145‑1 (MQ=2)
                                                                                                                                   ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt              >  6:325840/1‑117 (MQ=17)
                                                                                                                                            taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt              >  2:210878/1‑108 (MQ=35)
                                                                                                                                               |                                                                                                                    
TGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTC  >  NC_000913/3269775‑3270033

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