Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,296 T→C L119L (TTG→CTG)  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,2960TC98.0% 148.9 / ‑4.3 50L119L (TTG→CTG) yjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTA  >  NC_000913/4641056‑4641541
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
caacTGGAGCCAGCCAACCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGTAATATCATTACTACGCTACGCCAGTAGAACTGGTAAAGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACg                                                                                                                                                                                                                                               >  1:179066/4‑250 (MQ=255)
               aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGc                                                                                                                                                                                                                               >  1:232419/1‑251 (MQ=255)
               aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGAGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAACATCTTCATGGATGAGACATGCCGCGCTGGTGTTTGGTCGCGAAGCTTCCGAGCTGACTAACGAAGAGTTAGCGTTGGc                                                                                                                                                                                                                               >  2:333682/1‑251 (MQ=255)
                   gtgggggcGCACAGGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACCCTGGCTGAATCGTTACACGATGTCGATTTCCCTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGAc                                                                                                                                                                                                                           <  2:491701/244‑1 (MQ=255)
                            cacaTGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACt                                                                                                                                                                                                                  <  1:3880/251‑1 (MQ=255)
                            cacaTGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACt                                                                                                                                                                                                                  >  1:337024/1‑251 (MQ=255)
                                 ggATCGGGTGATATTAGTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGGCGCCACCACGGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGGCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGgtgt                                                                                                                                                                                                             <  2:70940/251‑1 (MQ=255)
                                       gatgaatttattGATAACTATAACGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTGCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACTCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAt                                                                                                                                                                                                       <  2:94522/244‑1 (MQ=255)
                                               tattGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGca                                                                                                                                                                                                <  1:135665/250‑2 (MQ=255)
                                               tattGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGGTTGGTCGCGAAGATTCCGGGCTGACTCACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAtggttgacg                                                                                                                                                                                               >  2:123695/1‑246 (MQ=255)
                                                attGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc                                                                                                                                                                                                 >  1:108888/1‑248 (MQ=255)
                                                 ttGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc                                                                                                                                                                                                 >  1:94365/1‑247 (MQ=255)
                                                  tGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGAt                                                                                                                                                                                             <  1:257135/250‑1 (MQ=255)
                                                              aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTCCTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCAAGGATGAGCCATGCCGCGCTGGTGTTTAGACGCGAAGATTACGGGCAGACTAACGAAGAGTTAGCGTTGGCTGACGATCTTACTGGTGTGCCACTGGGGGAAGATTATCCTTCACTc                                                                                                                                                                                 >  2:329035/1‑251 (MQ=255)
                                                              aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGa                                                                                                                                                                                                                                                     >  1:447721/1‑183 (MQ=255)
                                                              aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTc                                                                                                                                                                                 >  1:486673/1‑251 (MQ=255)
                                                              aaGTGTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGa                                                                                                                                                                                                                                                     <  2:447721/183‑1 (MQ=255)
                                                                gTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTACGGGCTGAATAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCa                                                                                                                                                                                >  2:113151/1‑250 (MQ=255)
                                                                 ttttCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGCCTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGCTGGTGGCAGATTATCCTTCGCTCAAt                                                                                                                                                                              >  2:214402/1‑251 (MQ=255)
                                                                    ccccGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTg                                                                                                                                                                           <  2:87157/250‑1 (MQ=255)
                                                                              ggCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAAGCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGca                                                                                                                                                                  <  2:303334/250‑2 (MQ=255)
                                                                               gCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAACCTGGGGCAGGcagt                                                                                                                                                                <  2:158460/251‑4 (MQ=255)
                                                                                  gAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTCTCCTTCGCTCAATCTGGGGCAGGCAGTGAt                                                                                                                                                             >  1:5750/1‑251 (MQ=255)
                                                                                   aaTCGTTACACGATGTCGATTTCCCTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTCCGCCAGTAGAACTGGTCCCGCTGTTAGAGGAAACATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGCGTTAGCGTTGGCCGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTgctg                                                                                                                                                            >  2:240347/1‑248 (MQ=255)
                                                                                          acacGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTAtt                                                                                                                                                     <  1:30366/251‑1 (MQ=255)
                                                                                           cacgagtacgagttCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTGAGAGGAAAAAGCTTCATGGATGAGCTATGCCGCGCTGGTGTTTGGTCGCGAAGATGCCGGGCTGACGAACGAAGAGTTAGCGTTGGCTGACGTTCTTACGGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGTGCAGGCAGTGATGGTCTATTg                                                                                                                                                    <  2:135974/239‑1 (MQ=255)
                                                                                             cGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCt                                                                                                                                                  <  2:108888/251‑1 (MQ=255)
                                                                                                    gATTTCACTGTCGTCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAag                                                                                                                                                                                                                                            <  1:279508/154‑1 (MQ=255)
                                                                                                        tCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGAAGTGATGGTCTATTGCTATCACTTAGc                                                                                                                                        >  2:357990/1‑250 (MQ=255)
                                                                                                         cACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGGTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCaa                                                                                                                                      <  1:333682/251‑1 (MQ=255)
                                                                                                                   accacTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATac                                                                                                                             >  1:244405/1‑250 (MQ=255)
                                                                                                                        tGCGCGCAGTTGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGGTTATTGCTTTCAAtttacaaactttataaaaaaa                                                                                                                        >  2:186998/1‑231 (MQ=255)
                                                                                                                         gcgcgcAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTTGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGTCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACaacaac                                                                                                                       <  2:162694/250‑1 (MQ=255)
                                                                                                                          cgcgcAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGGGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCg                                                                                                                     >  2:251473/1‑251 (MQ=255)
                                                                                                                             gcAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGACGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGGGTGCCGATGGTGGCAGATTATCCTTCGATCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCgacg                                                                                                                  >  1:54510/1‑248 (MQ=255)
                                                                                                                                   gggCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGt                                                                                                            >  1:235414/1‑251 (MQ=255)
                                                                                                                                    ggCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTg                                                                                                           <  1:360832/251‑1 (MQ=255)
                                                                                                                                      agAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTAATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTAAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAAAAACAGACGAAAAgcgac                                                                                                         >  1:371669/2‑247 (MQ=255)
                                                                                                                                                 tactacGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGc                                                                                                 >  1:149233/1‑248 (MQ=255)
                                                                                                                                                           cgccaGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGCTGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCa                                                                                     >  2:80523/1‑250 (MQ=255)
                                                                                                                                                              caGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAAAAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTg                                                                                 <  1:447569/251‑1 (MQ=255)
                                                                                                                                                                    aaCTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCt                                                                           <  1:65662/251‑1 (MQ=255)
                                                                                                                                                                                   tAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACgcgtc                                                            >  1:338181/1‑247 (MQ=255)
                                                                                                                                                                                    agagGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAgtc                                                            <  1:31352/250‑3 (MQ=255)
                                                                                                                                                                                         acaaaTCTTCAAGGATGAGTCATGCCGAGCGGGTGCTTGGTCGCGAAGATTCCGGGATGACTCACGAAGAGTTAGCGTTGGCTGTCGTTCTTACTGGTGTGGCGATGGGGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCAGAGATGGGCTATTGCTATCAATTAGCAAAATTAATGCAAAAACCGCAGAAAGGTGACACAACGCCAGACCACCATCAACTGCAAGCtttcagcaagcgagtcatggc                                                       >  2:21230/3‑232 (MQ=255)
                                                                                                                                                                                                   aTGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGAGAACGAGTCATGGCATTGCTgtcg                                             >  2:347005/1‑247 (MQ=255)
                                                                                                                                                                                                           gCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGActc                                         <  2:497378/246‑1 (MQ=255)
                                                                                                                                                                                                           gCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGActc                                         >  1:497378/1‑246 (MQ=255)
                                                                                                                                                                                                                           gtTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAAACGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATa                     >  1:20907/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAAATGCAAGCTTTACGAGAACGAGGCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGAATAAAAACTGGTCGACTGGTTa  >  2:308384/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
CACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTA  >  NC_000913/4641056‑4641541

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