Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,195 A→C 7.3% intergenic (+593/+21) 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,1950AC 7.3% 155.9 / 6.6 66intergenic (+593/+21)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (40/21);  new base C (3/2);  total (43/23)
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.23e-01

CCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270057‑3270343
                                                                                                                                           |                                                                                                                                                    
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)
                                                                                                                                           aaGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGGTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGcc  <  7:75120/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270057‑3270343

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: