Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,194 G→T 7.5% intergenic (+592/+22) 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,1940GT 7.5% 248.5 / 4.6 84intergenic (+592/+22)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (45/31);  new base T (4/3);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270058‑3270342
                                                                                                                                        |                                                                                                                                                    
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                                                                                                      <  3:64295/119‑1 (MQ=40)
                                                                  cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                      >  4:64295/1‑119 (MQ=40)
                                                                     tttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaa                                                                     >  7:396489/1‑149 (MQ=255)
                                                                       tctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                                                                           >  3:302753/1‑109 (MQ=35)
                                                                       tctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                                                                           <  4:302753/109‑1 (MQ=35)
                                                                        ctttTGGGTATAGAGTCGTGGACATTCATTAATATATCTAACCCAACAAAAGCAAAAACACGCCGAAGCGCGTTTTTACGTAAATCAGGTGAAACTTACCGATAAGACGGGTTCTGTCGTGGCCAGTCATTCATATATGCCAGccaaag                                                                  >  3:369628/1‑144 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTTGACAGTCATTCTTCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGTGTTTTTACGTAACTCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                         <  5:363821/141‑1 (MQ=18)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                                                                   >  1:394000/1‑131 (MQ=255)
                                                                         ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                                                                   <  2:394000/131‑1 (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                                                                                                       <  4:329100/105‑3 (MQ=38)
                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggtt                                                                                                       >  3:329100/1‑103 (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                                       >  3:287222/1‑149 (MQ=255)
                                                                                                   attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                       >  8:246076/1‑149 (MQ=255)
                                                                                                    ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                      >  5:339221/1‑149 (MQ=255)
                                                                                                     tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                           <  4:124558/111‑1 (MQ=255)
                                                                                                     tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                           >  3:124558/1‑111 (MQ=255)
                                                                                                     tcatCTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  3:75110/1‑149 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                            <  5:14748/125‑1 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                            >  6:14748/1‑125 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  5:125933/1‑148 (MQ=255)
                                                                                                      catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  8:40864/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                                                            >  8:133311/1‑124 (MQ=255)
                                                                                                       atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                            <  7:133311/124‑1 (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)
                                                                                                                                        |                                                                                                                                                    
CCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270058‑3270342

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: