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

TGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTC  >  NC_000913/3269775‑3270033
                                                                                                                                               |                                                                                                                    
tGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                                                 <  8:298038/148‑1 (MQ=255)
  cGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                                              <  7:376475/149‑1 (MQ=255)
  cGTGCTTTTGTAGGCAGGATAAGGAGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                                              <  1:324072/149‑1 (MQ=255)
    tGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                                                            <  4:324204/149‑1 (MQ=255)
     gCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAg                                                                                                           >  2:78955/1‑149 (MQ=255)
      cTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAg                                                                                                           <  2:371007/148‑1 (MQ=255)
         ttGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtca                                                                                                        <  5:161354/148‑1 (MQ=255)
            tAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                                                                     <  1:355335/148‑1 (MQ=255)
              ggCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                                                                  <  4:262557/149‑1 (MQ=255)
              ggCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                                                                  <  8:237875/149‑1 (MQ=255)
                 aGGATAAGGCGTTCACGCCGCATCCGGCATGAATACAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGCCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               <  6:22907/149‑1 (MQ=255)
                        ggCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGTGTGGTGGACAGGGATTGATCTTTCTgacgc                                                                                        >  5:119867/1‑145 (MQ=25)
                           gTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAg                                                                                                           <  7:211477/127‑1 (MQ=255)
                           gTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAg                                                                                                           >  8:211477/1‑127 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               >  1:228642/1‑138 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                               <  2:228642/138‑1 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           >  8:47569/1‑142 (MQ=255)
                            ttCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           <  7:47569/142‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGACCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGCTTCTTCTTTCTGCCCCTCCaa                                                                                   >  1:373250/1‑147 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  3:140027/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  3:347251/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  7:275234/149‑1 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                                                                   >  5:16226/1‑147 (MQ=255)
                               aCGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTAATTCATATTTCTGCCACTCaaa                                                                                   >  1:314442/1‑147 (MQ=21)
                                   cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTATCTTATGCGAATAGTGTAGTGGAAAATAATTAATATTTCTGACCCTCAAAAAGTa                                                                              >  2:154337/1‑148 (MQ=25)
                                    gcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                            <  3:305850/149‑1 (MQ=255)
                                       tCCGGCATGAATAAAGCGCACTAAACAATCTGAAACCAGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTCTATTGGAAAGTCAGAAATCTTTCAGCCACTCACAAGTTAAAAAc                                                                         >  1:77500/1‑149 (MQ=21)
                                            cATGAATAAAGCGCACTCAACAATCTGAAACCCTCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccg                                                                      <  8:248886/147‑1 (MQ=255)
                                            cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGACCCTCCAAAAAAAAAAAcccc                                                                      >  8:306983/1‑146 (MQ=255)
                                               gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg                                                                                                               >  8:378789/1‑103 (MQ=255)
                                               gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg                                                                                                               <  7:378789/103‑1 (MQ=255)
                                                aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTCAGAATTCGACACCTCAAAAAGAAAAAAAccgccggcg                                                                >  7:300532/1‑146 (MQ=18)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           <  4:322938/113‑1 (MQ=255)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                                                           >  3:322938/1‑113 (MQ=255)
                                                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                       <  7:42818/149‑1 (MQ=255)
                                                          aCTCAACTATCTGAAACCCGGCTGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                      <  8:72578/149‑1 (MQ=255)
                                                           cTCAACAATCACAAACCCGCCGTACCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  2:80269/149‑1 (MQ=11)
                                                            tctacaaTCTGAAACCCGCCGTAACGGTTTATTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                                                            <  6:289125/90‑1 (MQ=17)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                    >  6:91380/1‑149 (MQ=255)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACg                                                    >  1:80269/1‑149 (MQ=255)
                                                            tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGTGACGAGGACACGAATTAATAAGTCTGCCCAGAAAAAAGTAAAAAACCGCAGAACCGGGtgattaa                                                     >  7:288785/1‑142 (MQ=16)
                                                              aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                     <  5:138092/146‑1 (MQ=255)
                                                              aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                     >  6:138092/1‑146 (MQ=255)
                                                                 aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGAACGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCCTTCCTCTTTCTGCCCCTCCAAAAGAAAAAACACGCCAAAGGGGGTTTTTACGTaaaa                                               >  5:322470/1‑149 (MQ=17)
                                                                  aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                              <  3:68784/149‑1 (MQ=255)
                                                                  aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTaaaa                                               >  4:68784/1‑148 (MQ=255)
                                                                     tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGTGGACAGTCCTTAATCTGTATGCCCCTGCATAAGTAAAAAACCGCCGCATCGGGTTTATAAATAAAACAgg                                           >  5:220814/1‑149 (MQ=17)
                                                                       aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                            <  2:409858/146‑1 (MQ=255)
                                                                        aaCCCGCCGGAACTGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                        <  7:80626/149‑1 (MQ=255)
                                                                          cctgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGTGGAAAGTGATTGATATTTCTGCCCTTACAACAGCAAAAACCAGACGAAGAGGGTTTTTACGTATAAAAGGTGaa                                       >  7:373362/4‑148 (MQ=14)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  6:16226/149‑1 (MQ=255)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  5:232582/149‑1 (MQ=255)
                                                                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTATGGGTATAGTGCCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                      <  1:325607/149‑1 (MQ=255)
                                                                          ccagccgGAACGGTTTATTACGTACATCATGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCTGGTTTTTACGTAAAACAGGTGaaa                                      <  1:20284/146‑1 (MQ=255)
                                                                            cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                    <  5:91380/149‑1 (MQ=255)
                                                                                  aaCGGTATATTACGTACATCAGGTACAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                              <  4:324219/149‑1 (MQ=14)
                                                                                   actgTTTATTACTAACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                             <  1:362022/146‑1 (MQ=255)
                                                                                    cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                            <  5:112132/149‑1 (MQ=255)
                                                                                         taataCGTTCATCATGTAAATCTGTCCTATATTCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTATGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                        <  1:154337/145‑1 (MQ=14)
                                                                                           tttcGTACATAAGTTAATACTGACCGATAAGCCGCTTTCTTTTTGTTATAGTGTCGTGGACAGTGATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                     <  1:284106/146‑1 (MQ=14)
                                                                                               gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                                                            >  5:109731/1‑58 (MQ=255)
                                                                                               gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                                                            <  6:109731/58‑1 (MQ=255)
                                                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCt                                                                                                 >  4:247450/1‑65 (MQ=255)
                                                                                                    tcatttAAAACTGTCCGATAAGTAGTTTTCATTTTTTGTTTGTGTCGTGGACAGTCTTTCTTATATCTGCCCCTCCAAAAGTACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg            <  2:77296/144‑1 (MQ=11)
                                                                                                    tcatgcaaaaCTGTCCGATAAGCCGCTTTCTTTAGGGTATACTGGCTTGGACACTCTTTGCTCTTTCTCCCTCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg            <  8:162478/143‑1 (MQ=11)
                                                                                                        gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 >  3:349542/1‑76 (MQ=34)
                                                                                                        gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                 <  4:349542/76‑1 (MQ=34)
                                                                                                          aataattaCCGTTAAGCTGATTGCTTTTGTCTATTGTGTCGACGACAGGCTTAATTCTATCTCCCACTCCAAAATTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg      <  1:299671/142‑1 (MQ=2)
                                                                                                              gagACCTATTATCCGCATTCATTAGGGCATAGCGTCGTGGACACTCCTTCTTCAGTCAGCCCCTGCATCAGATATAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc  <  6:392511/147‑1 (MQ=2)
                                                                                                              ctttccGATATGTCGATTTCAGTTGTCTCTAGTGTCGTCGACAGCGTTTTCTCTTTCTCCCACTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc  <  7:295649/145‑1 (MQ=2)
                                                                                                                                   ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt              >  6:325840/1‑117 (MQ=17)
                                                                                                                                            taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt              >  2:210878/1‑108 (MQ=35)
                                                                                                                                               |                                                                                                                    
TGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTC  >  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: