Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,196 A→G 6.7% intergenic (+594/+20) 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,1960AG 6.7% 254.0 / 7.5 91intergenic (+594/+20)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (50/34);  new base G (4/2);  total (54/37)
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.32e-01

CCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270058‑3270343
                                                                                                                                          |                                                                                                                                                   
cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTTCCCCACCAAAAGCAAAAACCCGCTTAGGCGGGgtttt                                                                                                                                           >  2:138091/1‑149 (MQ=11)
   tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgag                                                                                                                                        >  7:353873/1‑146 (MQ=11)
   aacAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGCGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATCGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                                                                                                                                           <  8:353873/144‑1 (MQ=11)
    ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCATTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                        >  7:315799/1‑148 (MQ=37)
     cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGCCCGATAAGCCGCTTTCTTTTGGGTATAGCGTAGTGTAAAGGCATTCAACTTTCGGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGt                                                                                                                                       >  8:255944/1‑148 (MQ=11)
      aaacGCAAACACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAACCTGACCGATACGCCGCTTTCATTTGGGTATAGCGGCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgctttt                                                                                                                                     <  1:419791/149‑7 (MQ=11)
                                      gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGTGTATAGCGTCTAGGACAGTCAATCATAATAATCACCCAACACAAGCAAAAACCCCCAGAAGCGGGAATTTACTGTAATAATCTGTAAATGACCGATAATAAGGGTTCt                                                                                                     >  2:95219/1‑149 (MQ=2)
                                       gaacccAGGAGAATCTGACCTATATCCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTg                                                                                                    <  7:415956/145‑1 (MQ=14)
                                        aaaCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                              <  6:247005/122‑1 (MQ=38)
                                         aaCCAGGTGAAACTGACCGATAAGCCGCGATCTTATCTTGCTAGCGGCGTGGACAGTCATTCCTCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTc                                                                                                  <  5:262614/149‑1 (MQ=14)
                                            cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCATGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                                                                                                                     >  6:70431/1‑127 (MQ=25)
                                             attttaaaCTGACCTATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTgg                                                                                              <  1:247283/144‑1 (MQ=255)
                                             aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                              <  6:394131/117‑1 (MQ=35)
                                             aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                              >  5:394131/1‑117 (MQ=35)
                                                tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                                                                  <  7:114477/94‑1 (MQ=33)
                                                   aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCTAGGAAAGTCATTAATATTTAAGACCCTACAAAAACAAAAACCCCCCCAAGCCGGATTTTAAGTAAATCACGTTCAACTGACAAAAAAGCAGGGTTCTGGAGTGGAAAGt                                                                                         >  8:439538/1‑148 (MQ=9)
                                                    aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCATTAATCTTTCAGCCCCACCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTACGTAAATCAGGTGAAACTGACCCATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                       >  7:34472/1‑149 (MQ=17)
                                                     cttaCCGATAAGCCGCTTTCTTTTGGGTAACCCTTCGTGGACTGTCTTTCATCTTTCTGCCCCTCCTAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcat                                                                                      <  7:183907/146‑1 (MQ=17)
                                                     attaactaTAAGACGATTTCTTTAGCGAATAGCGTCTGTGACAGTCATTCATCATTCTGCCCCTCCAAAAGCAAAAACCCTCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcat                                                                                      <  7:176835/142‑1 (MQ=17)
                                                          cGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                       <  5:224417/143‑1 (MQ=255)
                                                          cGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTgggccggca                                                                                       >  6:224417/1‑136 (MQ=255)
                                                            aTAAGCCGCCTTCATTTGGGTATAGCGTCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTa                                                                               <  7:95619/149‑1 (MQ=255)
                                                                  cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                       >  4:64295/1‑119 (MQ=40)
                                                                  cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                       <  3:64295/119‑1 (MQ=40)
                                                                     tttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaa                                                                      >  7:396489/1‑149 (MQ=255)
                                                                       tctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                                                                            <  4:302753/109‑1 (MQ=35)
                                                                       tctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                                                                            >  3:302753/1‑109 (MQ=35)
                                                                        ctttTGGGTATAGAGTCGTGGACATTCATTAATATATCTAACCCAACAAAAGCAAAAACACGCCGAAGCGCGTTTTTACGTAAATCAGGTGAAACTTACCGATAAGACGGGTTCTGTCGTGGCCAGTCATTCATATATGCCAGccaaag                                                                   >  3:369628/1‑144 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTTGACAGTCATTCTTCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGTGTTTTTACGTAACTCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                          <  5:363821/141‑1 (MQ=18)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                                                                    <  2:394000/131‑1 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                                                                    >  1:394000/1‑131 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc                                                                  >  3:342986/1‑149 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGCGTTTTTACGTAAATCAGGTGAAACTGACCGATAACCCGGGTTCTGTCGGGGACAATCATTCATCTTGGCCAGCAATcc                                                                   >  8:96270/1‑147 (MQ=18)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCAGACCCAACAAAAGCAAAAACCCGCCGAAAAGGGTGTTTAAGTAAATAAGGAGAAAATGACCGATAAGACGGGTTCAGTCGTAGACAGTGAGTGCTCTTGGACAGCGATCg                                                                   >  6:363821/1‑148 (MQ=9)
                                                                         atttGTGTATAGCGTCGTGGACAGTCATTCATCTATCTGCCCCTCCAATTGCAAAAACCCGAGGTAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCg                                                                   <  1:95219/147‑1 (MQ=21)
                                                                             gggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACTGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCa                                                               <  8:315799/148‑1 (MQ=255)
                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggtt                                                                                                        >  3:329100/1‑103 (MQ=38)
                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggtt                                                                                                        <  4:329100/105‑3 (MQ=38)
                                                                                   aGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCt                                                         >  3:335292/1‑148 (MQ=255)
                                                                                     tgtcgtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                                                                          <  7:191231/96‑1 (MQ=18)
                                                                                     cttaggtgACTGTAATTCATCTTTCTGCACGTCCTATAGCAAAAACCCGCCTAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCaa                                                      <  7:188361/142‑1 (MQ=21)
                                                                                      gtcgtGGACAGTCATTCATCTTTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCCGCTTTTTACGTAAATCAGGTGAAACTGACCGAGAAGCAGGGTTCTGTCGTGGAAAGTCATTAATCTAGGCCAGCAATCGCTAACTGGATCAag                                                     >  2:30183/1‑149 (MQ=18)
                                                                                         gtgtACAGTCATTCATATTTCTGCCCCTCCATAAGCTAAAACCGGCCGATGCTGGTTTTTACGTATATCAGGTGAAACTGACCGATAAGCGGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagca                                                   <  1:260455/148‑1 (MQ=25)
                                                                                              cAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAc                                             >  1:399148/1‑149 (MQ=255)
                                                                                               aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                               <  7:60954/98‑1 (MQ=255)
                                                                                               aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACACGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                               >  8:60954/1‑98 (MQ=255)
                                                                                                   attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                        >  8:246076/1‑149 (MQ=255)
                                                                                                   attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                        >  3:287222/1‑149 (MQ=255)
                                                                                                    ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                       >  5:339221/1‑149 (MQ=255)
                                                                                                     tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                            >  3:124558/1‑111 (MQ=255)
                                                                                                     tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                            <  4:124558/111‑1 (MQ=255)
                                                                                                     tcatCTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                      >  3:75110/1‑149 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                             <  5:14748/125‑1 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                             >  6:14748/1‑125 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                      >  8:40864/1‑148 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                      >  5:125933/1‑148 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                              >  5:52959/1‑59 (MQ=35)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                              <  6:52959/59‑1 (MQ=35)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt                                                                                                     <  4:219107/84‑1 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt                                                                                                     >  3:219107/1‑84 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                             <  7:133311/124‑1 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                             >  8:133311/1‑124 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCaa                                                      >  4:352284/1‑131 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCaa                                                      <  3:352284/131‑1 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                      >  5:68609/1‑147 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                      >  7:387498/1‑147 (MQ=255)
                                                                                                         cTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                          <  6:11263/109‑1 (MQ=255)
                                                                                                         cTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                          >  5:11263/1‑109 (MQ=255)
                                                                                                                ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACgg                              <  7:506911/146‑1 (MQ=255)
                                                                                                                 cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTc                                                        >  2:118111/1‑119 (MQ=255)
                                                                                                                 cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTc                                                        <  1:118111/119‑1 (MQ=255)
                                                                                                                 cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACgg                              >  8:444774/1‑145 (MQ=255)
                                                                                                                 cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACgg                              <  7:444774/145‑1 (MQ=255)
                                                                                                                 cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCg                          >  2:7140/1‑149 (MQ=255)
                                                                                                                  ccTCCAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCCCTGGCTCAAGCAGCATACCCGGGTTCATTACGGGCCg                          >  2:150400/1‑148 (MQ=255)
                                                                                                                  ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCg                          >  8:417520/1‑148 (MQ=255)
                                                                                                                   cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAACCTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCACTCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTa                        >  1:75522/1‑149 (MQ=255)
                                                                                                                     ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAGGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAcc                      >  1:198981/1‑149 (MQ=255)
                                                                                                                       aaaaTCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCtt                    <  8:396489/149‑1 (MQ=255)
                                                                                                                       aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                     >  5:243191/1‑148 (MQ=255)
                                                                                                                       aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                     >  8:4468/1‑148 (MQ=255)
                                                                                                                           gtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg                 >  4:152927/3‑148 (MQ=255)
                                                                                                                           gtaaaaaCCCGCTTCGGCGGGTTCTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg                 >  2:397800/3‑148 (MQ=39)
                                                                                                                                  cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTa         >  2:151297/1‑149 (MQ=255)
                                                                                                                                   ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt        <  7:4468/149‑1 (MQ=255)
                                                                                                                                    cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAtt       >  4:225617/1‑149 (MQ=255)
                                                                                                                                     gccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAccc            <  4:318884/143‑1 (MQ=255)
                                                                                                                                     gccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAccc            >  3:318884/1‑143 (MQ=255)
                                                                                                                                     gccgAAGCGGGTTTTTACGTAAATCAGCTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATATAGGCCAGCAAACGCTCACTGGCTAAAGCAGCCTACCCCGGTTCAGTACGGGCCCTACCATATCAACCCCTAttt      >  3:117854/1‑149 (MQ=38)
                                                                                                                                        gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGc   <  5:275432/149‑1 (MQ=255)
                                                                                                                                        gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGc   >  1:274358/1‑149 (MQ=255)
                                                                                                                                         aaGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGcc  <  6:243191/149‑1 (MQ=255)
                                                                                                                                         aaGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGcc  <  4:287222/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270058‑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: