Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,193 C→T 6.7% intergenic (+591/+23) 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,270,1930CT 6.7% 168.8 / 2.6 62intergenic (+591/+23)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (37/20);  new base T (2/2);  total (40/22)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTT  >  NC_000913/3270057‑3270339
                                                                                                                                         |                                                                                                                                                  
ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCTGGGAGAGGCCTTCATATTTAGGCGCCGCCAAAAACAAAAACCCGCCGAAGAgggggtg                                                                                                                                        >  8:122208/1‑145 (MQ=11)
   cTCCAAAAGCAAACACCCGCCGCACCGGGTTATTACGGAAACCCGGTGACACTCTCCGAGACGCCGCCTGCTTCTCGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgc                                                                                                                                     <  7:121955/149‑3 (MQ=2)
   cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                     >  7:349168/1‑149 (MQ=37)
    tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCAATCATCTTTCAGCACCGACAAAAGCAAAAACCCGCCTAAGCAGGTGGTTACg                                                                                                                                    >  2:24653/1‑149 (MQ=11)
       aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                                                 <  8:349168/149‑1 (MQ=37)
         aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                                                                                                                                       >  8:323657/1‑141 (MQ=255)
         aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                                                                                                                                       <  7:323657/141‑1 (MQ=255)
                aaCCCGCCGAAGCGGGTTTTTACGGAAACCAGGTGAAACTGACGGATACGCCGCTTTCTTTTGGGTATATCTGCGTGGACATTTATTCATCTTTCTTCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                                                       <  6:334441/134‑1 (MQ=11)
                aaCCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCTATAAGCAGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAACCAAAAACCCGCCGAAGAGGGgtttt                                                                                                                                       >  5:334441/1‑134 (MQ=11)
                 aCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTAGAGCGTCGTGGCCAGTCATTCATCTTTCTGCGCCGGCAATAGCAAAAACCCGCCGAAGCGGGTGTTTTAGGAAATGAGGTGa                                                                                                                        >  1:225507/1‑148 (MQ=11)
                                tttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTAGTGGAAAGTCATTCATGTTTATGCGCCTCCAAAAGCAATAACCCGCGGAAGGGGGTTTTTAAGGTAATCAGGTGAAACTGACCGATAATcc                                                                                                        >  7:194322/1‑149 (MQ=11)
                                                   aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCGGCCCCTGCAAAAGCAAAAACCCGCCGCATCGGGGTTTTACGTAAATCAGGTGCAACAGACCGATAAGCCGTGTTCTGTAGTGGACAGt                                                                                     >  1:79665/1‑149 (MQ=14)
                                                         aCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCAAGTCATTCATCGTTCGGCCCCACCAAAAGCAAAAACCCGCCGAAGAGCGGGTGTACGTAAAACAGGTGCAACTGACCCAGAAACCGGGTTCGTTCGGGGAAAGtcattc                                                                                >  6:117426/1‑148 (MQ=9)
                                                          ccGATATGGCGCTTTCTTTTGGGTATACGTGAGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattcat                                                                              <  7:58056/149‑1 (MQ=17)
                                                            gATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                                                                             >  1:8331/1‑148 (MQ=255)
                                                                aGCCGCTTTCTTTTGGGTATAGCGTCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                        <  1:203807/149‑1 (MQ=255)
                                                                 gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTATGTCGTGGACAGTCATTCATATATGc                                                                        >  2:42744/1‑148 (MQ=255)
                                                                 gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCATCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATAAGGTGAAACTGACCGATAATCAGGGTTATTTAGTGGACAGTCATTCCTCTAGGcc                                                                       >  2:203807/1‑149 (MQ=17)
                                                                     ctttctttTGGGTATAGCGTCGTGGACATTAATTCATCTTACTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATAAGCTGAAACTGACCGATAAGCCGGGTTCTGTCTTGGACAGTCATTCATCTAGGCCAGc                                                                    >  7:262759/1‑148 (MQ=18)
                                                                          ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTAAGTAAATAAGGTGAAACTGACCCATAAGCCGGGTTCTGTCGTGGACAGTCATTCCTCTAGGCCAGCAATAGc                                                              >  8:18958/1‑149 (MQ=18)
                                                                            ttGGGTATAGTGTAGTTTACATTCAGTCATCTTTCTGACCATCCAAAAGCAAAAACCCGACTAAGATGGTTTTTACGTAAATCACTTTAAACTGACCGATAAGGCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc                                                              <  4:345366/147‑1 (MQ=18)
                                                                               ggTATAGCGTCGTGGACAGTCATTCATCTTTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTAATCTAGGCCCGCAATCGCTCCCt                                                         >  4:118456/1‑149 (MQ=255)
                                                                                   taGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCCGCAATCGCTCACTGGCt                                                     >  2:114543/1‑149 (MQ=255)
                                                                                       gtcgtGGACAGTCATTCATCATTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAag                                                 <  1:42744/149‑1 (MQ=255)
                                                                                        tcgtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAag                                                 >  2:274112/1‑148 (MQ=255)
                                                                                             ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                    <  6:10151/77‑1 (MQ=34)
                                                                                             ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                    >  5:10151/1‑77 (MQ=34)
                                                                                               aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTa                                          >  7:172129/1‑149 (MQ=255)
                                                                                                 aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAcc                                        >  1:10182/1‑149 (MQ=255)
                                                                                                    tattcatCTTTCTGCCACTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCgg                                     <  3:218691/148‑1 (MQ=255)
                                                                                                     attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                    >  3:167625/1‑149 (MQ=255)
                                                                                                      ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                           >  1:351667/1‑93 (MQ=255)
                                                                                                      ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                           <  2:351667/93‑1 (MQ=255)
                                                                                                      ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCgc                                     >  4:218691/1‑146 (MQ=255)
                                                                                                      ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                   >  5:326406/1‑149 (MQ=255)
                                                                                                      gtcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                   >  7:268157/2‑149 (MQ=255)
                                                                                                        catCTTTCTGCCCCTGCAAAAGCAAAAACCCGCCGAAGCTGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTGAAGCAGCCTACCCGGGtt                                  <  8:262759/148‑1 (MQ=255)
                                                                                                         atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                  >  7:176365/1‑147 (MQ=255)
                                                                                                           cTTTCTTCCCCTCCAAAATGAAAAACCCGCCGAAGCTTGTTTTTACGTAAATGAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                                 <  2:79665/146‑1 (MQ=37)
                                                                                                           cTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                               >  2:213366/1‑148 (MQ=255)
                                                                                                              tcgccctcTCCAAACGCAAAACCCCGCCGAAGCGGGTTTTTACGTAAATCCGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                <  1:208624/137‑1 (MQ=37)
                                                                                                                 gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGc                        >  3:25345/1‑149 (MQ=255)
                                                                                                                   cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                >  8:258247/1‑139 (MQ=255)
                                                                                                                   cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                <  7:258247/139‑1 (MQ=255)
                                                                                                                   cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCg                      >  7:243976/1‑149 (MQ=255)
                                                                                                                    ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                           <  6:238856/79‑1 (MQ=255)
                                                                                                                    ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                           >  5:238856/1‑79 (MQ=255)
                                                                                                                     cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACCGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTa                    >  8:210741/1‑149 (MQ=255)
                                                                                                                      tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAc                   >  6:80987/1‑149 (MQ=255)
                                                                                                                        cAAAAGCAAAAACGCGCCGAAGCGGGTTTTTACGTAAATCAGGTGATACTGACTGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAcc                  <  2:110345/148‑1 (MQ=255)
                                                                                                                        cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                 >  4:298318/1‑149 (MQ=255)
                                                                                                                        cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                 >  3:188341/1‑149 (MQ=255)
                                                                                                                        cAAAAGCAAAAACCCGCCGAAGCGGGATTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCAGTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                 >  6:284039/1‑149 (MQ=255)
                                                                                                                         aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCtt                <  2:8331/149‑1 (MQ=255)
                                                                                                                           aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt              <  1:213366/149‑1 (MQ=255)
                                                                                                                           aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTGCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt              <  1:57530/149‑1 (MQ=255)
                                                                                                                             gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGa            >  3:144769/1‑149 (MQ=255)
                                                                                                                             gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGACGTACCTTATg             >  6:48617/1‑148 (MQ=255)
                                                                                                                              taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa           >  3:73800/2‑149 (MQ=255)
                                                                                                                              cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGATCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa           >  2:3220/1‑149 (MQ=255)
                                                                                                                                 aaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                                                                                                       <  4:273358/54‑1 (MQ=35)
                                                                                                                                  aaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCTATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAcccc       <  3:305286/149‑1 (MQ=255)
                                                                                                                                     ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGTACAGTAATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCATACCCGGGTTCAGTACCGGCCGTCCCGTATGCACCCCTAt    >  7:346625/1‑149 (MQ=38)
                                                                                                                                      cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAtt   >  1:50426/1‑149 (MQ=255)
                                                                                                                                       gccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAttt  >  3:296376/1‑149 (MQ=255)
                                                                                                                                       gccgAAGCGGGTTTTTACGTAAAACAGGAGAACCTGACCGATAAGCCGGGTGCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAttt  <  5:117426/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTT  >  NC_000913/3270057‑3270339

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: