Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,382 G→A A147A (GCG→GCA yjtD → putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,641,3820GA98.1% 162.5 / NA 53A147A (GCG→GCAyjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/32);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACC  >  NC_000913/4641136‑4641622
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
tGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTg                                                                                                                                                                                                                                               >  1:114729/1‑250 (MQ=255)
   aTCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATgg                                                                                                                                                                                                                                           <  2:537279/251‑1 (MQ=255)
     cGTTACACGATGTCGATTTCACTGTCGCCACCACGTCGCGCAGTCGGGCGAAATATCATTACGACGTTACGCCAGGAGAACTGGTACCGCTGTGATAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATGCCGGGTGGAGTAACGAAGAGTTAGCGTGGGCTGACGTGCTTACTGGGGTGGCGAGGGGGGAAGATTATCCGTCGCTCAATTTGGTGCAGGCAGCGATGGt                                                                                                                                                                                                                                          <  2:403553/250‑1 (MQ=255)
      gTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCt                                                                                                                                                                                                                                        <  1:174797/251‑1 (MQ=255)
       ttACACGATGCCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCt                                                                                                                                                                                                                                        >  1:196295/1‑250 (MQ=255)
        tACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTa                                                                                                                                                                                                                                       <  2:403144/250‑1 (MQ=255)
                  cGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCGGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTACCGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAAt                                                                                                                                                                                                                            >  2:610400/1‑251 (MQ=255)
                   gATTTCACGGTCGCCACCAATGCGCGCAGTCGGGTGAAATATCATTACTACGCTACGCCAGTAGAATTGGTACCGCTGTTAGAGGAAAAATCTTCAGGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGTGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATTCTTCGCTCAATCTGGGGCAGGCAATGATGGGCGATGGCTATCAAtg                                                                                                                                                                                                                           <  2:390440/251‑2 (MQ=255)
                        cacggtcgactccacTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGGTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGGCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCTGATGGTGGCAGATTCTCTTTCGCTCAATCTGGGGCAGGCGGTGATGGGCTATTGCGATCAATTAGCaa                                                                                                                                                                                                                      <  2:110734/240‑1 (MQ=255)
                                     acTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATacaaca                                                                                                                                                                                                         <  1:162424/251‑1 (MQ=255)
                                           cgcAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCgac                                                                                                                                                                                                   <  1:3815/251‑3 (MQ=255)
                                           cgaaGTCGGGAGAAATATCATTAATACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGATGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTATGCCGATGGTGGCAGATTATCCTTCGCTCAATCTAGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCACCATTAATACAACAACCgac                                                                                                                                                                                                   >  2:554643/4‑249 (MQ=255)
                                                     cGAAATATCATTACTACGCTACGCCAGTAGAACTGGTCCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCAGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAACGTGac                                                                                                                                                                                         >  2:566946/1‑250 (MQ=255)
                                                         atatCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGacacaa                                                                                                                                                                                     <  2:438484/251‑6 (MQ=255)
                                                           atCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGacacaac                                                                                                                                                                                    <  1:560898/250‑7 (MQ=255)
                                                           atCATTAATACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACAGACGAAAAGTGacacaacg                                                                                                                                                                                   >  2:448577/1‑244 (MQ=255)
                                                                                  gAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGc                                                                                                                                                            <  1:562557/251‑1 (MQ=255)
                                                                                   aaCTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTg                                                                                                                                                                                                                                    <  1:410061/178‑1 (MQ=255)
                                                                                   aaCTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTg                                                                                                                                                                                                                                    >  2:410061/1‑178 (MQ=255)
                                                                                             gCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTg                                                                                                                                                                                                                                               >  2:482110/1‑157 (MQ=255)
                                                                                             gCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTg                                                                                                                                                                                                                                               <  1:482110/157‑1 (MQ=255)
                                                                                                gTTAGAGTAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGCCCCAACGGCAGCCCAACATCAACTGCAAGCTTTACGCGAAcc                                                                                                                                                >  2:26681/1‑248 (MQ=255)
                                                                                                    gagGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCAt                                                                                                                                          <  1:129547/251‑1 (MQ=255)
                                                                                                          aaaTCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCAtt                                                                                                                                    <  1:103095/251‑1 (MQ=255)
                                                                                                                tCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTga                                                                                                                               <  1:297865/250‑1 (MQ=255)
                                                                                                                tCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGGCATGGCATTGc                                                                                                                                  <  1:130806/247‑1 (MQ=255)
                                                                                                                    ggatgcaccATGCCGGGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGTCTAACGAAGAGTTAGCGTTGGCTGACGTTCGTACTGGTGTGCCGATTGTGGCAGATTATCCTTCGCACAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCTTTGCTgacgac                                                                                                                           <  1:446090/243‑1 (MQ=255)
                                                                                                                       tGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTAGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAAAAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCAGGGCATTGCTGACGACTCtg                                                                                                                       >  2:482207/1‑251 (MQ=255)
                                                                                                                          gCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCtggca                                                                                                                    <  2:335605/251‑1 (MQ=255)
                                                                                                                                    gcTGGTGTTTGGTCGAGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGATCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATg                                                                                                          >  1:446625/1‑251 (MQ=255)
                                                                                                                                       ggTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGTGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACa                                                                                                       <  2:148497/251‑1 (MQ=255)
                                                                                                                                               gTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTg                                                                                               <  2:613308/251‑1 (MQ=255)
                                                                                                                                                 cgcgAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGt                                                                                             >  1:394193/1‑251 (MQ=255)
                                                                                                                                                    gAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCACTAAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTCGCAACATTAATACAACAACAGCCGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTCAGCGAACGAGTCATGGCATTGCTGACGAATCTGGCAGTGGCAGATGACATAAAACTGGTCg                                                                                           >  1:605152/1‑250 (MQ=255)
                                                                                                                                                      aGATTCCGGGCTGACGAACGAAGGGTTGGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACt                                                                                        <  1:334330/251‑1 (MQ=255)
                                                                                                                                                                   aCTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACg                                                                           <  1:448577/251‑1 (MQ=255)
                                                                                                                                                                      aaCGAAGAGTTAGCGGTGGCTGACGTTCTTACGGGGGTGCTGATGGTGGCAGATTATACTGCGCTCAATCTGGCGCAGGCAGTGATGGGCTAGTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATGGCTGACGACTCTGGCAGTGGCAGTTGACATAAAACTGGTCGACTGGTTACAGCAACGCCt                                                                        <  2:105790/251‑1 (MQ=255)
                                                                                                                                                                       aCGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCCTTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGtt                                                                                    >  2:540113/1‑238 (MQ=255)
                                                                                                                                                                       aCGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGtt                                                                                    <  1:540113/238‑1 (MQ=255)
                                                                                                                                                                                tAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGAGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGCCATAAAACTGGTCGCCTGGTTACAGCAACGCCTGGGGCtttt                                                               >  2:538925/1‑250 (MQ=255)
                                                                                                                                                                                       ggCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCaa                                                        <  2:209162/250‑1 (MQ=255)
                                                                                                                                                                                             cGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGaca                                                 <  1:390702/251‑1 (MQ=255)
                                                                                                                                                                                                tCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACgg                                              >  1:163746/1‑251 (MQ=255)
                                                                                                                                                                                                tCTTACTGGAGTGCCGATCGAGGCAGATTATCCTTCGCTCCATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGAGACATTAATAAAGCAACCGACGAAAAGTGACCCAACGTCAGACCAACATCAACTGCACGCTTTACGCGAACGAGTCATGGCATTGCTGCCGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACgg                                              >  2:278939/1‑251 (MQ=255)
                                                                                                                                                                                                       ggTGTGCCGATGGTGGCATATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGtt                                       <  1:610400/251‑1 (MQ=255)
                                                                                                                                                                                                       ggTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAgtc                                                                                                                                            >  1:56096/1‑148 (MQ=255)
                                                                                                                                                                                                       ggTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAgtc                                                                                                                                            <  2:56096/150‑3 (MQ=255)
                                                                                                                                                                                                             ccGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCg                                 <  1:437795/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggtggCAGATTATCCTTCGCTAAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgc                            >  1:141822/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cGCTCAAGTTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGGCATGGCAGTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGaaaa            <  2:458450/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAata       <  1:521635/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCAcc  >  2:485060/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggggCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATc     <  1:602121/247‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
TGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACC  >  NC_000913/4641136‑4641622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: