Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,640,908 G→T intergenic (+137/‑34) ytjE → / → yjtD protein YtjE/putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,640,9080GT93.0% 126.8 / NA 43intergenic (+137/‑34)ytjE/yjtDprotein YtjE/putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (23/17);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATG  >  NC_000913/4640666‑4641152
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
tttAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTcgcg                                                                                                                                                                                                                                                <  1:300487/250‑1 (MQ=255)
   aataatCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGtt                                                                                                                                                                                                                                             <  1:616995/250‑1 (MQ=255)
      aatCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTAtt                                                                                                                                                                                                                                         <  2:16088/250‑1 (MQ=255)
         ccAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTgg                                                                                                                                                                                                                                      <  1:90986/250‑2 (MQ=255)
             cTAGTTGCATCATACAACTAATACACGGGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTCAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCagg                                                                                                                                                                                                                                 >  2:204285/1‑249 (MQ=255)
                     atcatACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTgg                                                                                                                                                                                                                         <  2:6144/251‑1 (MQ=255)
                                aataaACGAGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACACTTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAACTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTCCATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTAGCGGTTATTTGGTCAAGAAGAGTGGCACTATGCGt                                                                                                                                                                                                               >  1:403553/1‑250 (MQ=255)
                                       gtggtgAATCCAATTGTCGCGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACAGCGTAATAAGTTGACAGTGATCACCCTGTTc                                                                                                                                                                                                                                                   >  2:375916/1‑208 (MQ=255)
                                       gtggtgAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTc                                                                                                                                                                                                                                                   <  1:375916/208‑1 (MQ=255)
                                                             tttatttTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGcc                                                                                                                                                                                    >  1:337615/1‑248 (MQ=255)
                                                                                    aaGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcggc                                                                                                                                                          >  2:390702/1‑251 (MQ=255)
                                                                                              aaGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCGTTTTTAATCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGTTTATTTGGTCAATAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATg                                                                                                                                                <  2:165539/251‑1 (MQ=255)
                                                                                                             tttttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTCCGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATCTTGGGGCCGCGGCGCGGGCAATGAAAACGATggtattt                                                                                                                                 >  1:390440/1‑246 (MQ=255)
                                                                                                                ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGCAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGc                                                                                                                              >  1:530579/1‑251 (MQ=255)
                                                                                                                ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGTGGCAGCGGCTCGGGCaa                                                                                                                                                  <  1:180207/231‑1 (MQ=255)
                                                                                                                ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCaa                                                                                                                                                  >  2:180207/1‑231 (MQ=255)
                                                                                                                 tAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGc                                                                                                                              <  1:303010/250‑1 (MQ=255)
                                                                                                                     acttCAACCGTTAGTACCATCAGGAAATAGTTTCGCCTTTTTTACGCTAAGTACAGGGATTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACAGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCAGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGAGGGGTTTTAGCGaact                                                                                                                         >  2:437795/4‑248 (MQ=255)
                                                                                                                       gtCAACCGTTAGTACATTCAGTAAAGAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGctaactg                                                                                                                        <  2:28716/249‑7 (MQ=255)
                                                                                                                          aaCCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCgg                                                                                                                     >  1:16731/1‑250 (MQ=255)
                                                                                                                                   tACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGAt                                                                                                            >  1:50148/1‑250 (MQ=255)
                                                                                                                                                          ttttttAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGagccagcc                                                                                    >  1:349534/1‑251 (MQ=255)
                                                                                                                                                           tttttAAGCAAAGTAAAGGGATTTTTCTGCGACTTACTTTAAGGATTTGTAAATTCGCGCCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTAGTTGGTCAAGAAGAGTGGCAATATGCTTATAACGATTATTCTGTTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGATAACGATGGGTTTTAGCGAACTGCGGATTGTCGAGTGGCAGGCACACCTGGaaccagc                                                                                     <  2:563152/249‑1 (MQ=255)
                                                                                                                                                               tAAGCTAAGTAAAGGGCGTTTTCTGCGCCTTACGTTAAGGATTTGTAAATTCGAACCGCGTAATAAGTTGACAGTGACCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAAGATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCAccc                                                                                <  2:143946/250‑1 (MQ=255)
                                                                                                                                                                 aGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAACTATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGAGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCt                                                                             >  2:3815/1‑251 (MQ=255)
                                                                                                                                                                   cTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTgg                                                                           >  2:141697/1‑251 (MQ=255)
                                                                                                                                                                     aaGTAAAGGGATTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCATAGCAGAAACTATTGGGGCAGCGGCGCGGGCAATGAACACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCAGGAGCCAGCCACCCGCTgggg                                                                         >  1:110734/1‑250 (MQ=255)
                                                                                                                                                                       gTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                                                                      <  1:375327/124‑1 (MQ=255)
                                                                                                                                                                       gTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                                                                      >  2:375327/1‑124 (MQ=255)
                                                                                                                                                                          aaGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGcac                                                                    >  1:403144/1‑251 (MQ=255)
                                                                                                                                                                             ggCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATg                                                                 <  1:135122/251‑1 (MQ=255)
                                                                                                                                                                             ggCTTTTTATGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATg                                                                 <  2:213796/251‑1 (MQ=255)
                                                                                                                                                                                  tttCGGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGGTCGCGGTTATTGGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACTCCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATc                                                             <  2:530579/250‑1 (MQ=255)
                                                                                                                                                                                   ttCTGCGACTTACGTTAAGGATTTGTAAGTTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTg                                                           >  2:162424/1‑251 (MQ=255)
                                                                                                                                                                                    actggtaCTTACGTTAAGTATTTGTAAATTCGCACCGTGTAATAAGTTGACATTGATCCCCCTTTTCGCGGTTATTTGTTCAAGAAGAGTGGTAATATGGTTAGAACGATTATTCGGGGCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAAGGAAAACGATGGGGTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGTAAGCCACCCGCTGGGTCGcacaatgaactg                                                           <  2:262232/244‑9 (MQ=255)
                                                                                                                                                                                        cGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGCGCCAGCCACCCGCTGGCTCGCACATGGATCTGGTgct                                                      >  2:148228/1‑249 (MQ=255)
                                                                                                                                                                                                 ttAAGGATTTGTAAATTAGCACCGCGTAATACGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGATGGGTCGCACATGGATCTGGTGATATTATTGat                                             >  1:458450/1‑251 (MQ=255)
                                                                                                                                                                                                         gtGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGTTTCGCGGTTATTTTGTCAAGAAGAGTGGAAATATGCGTATAACGATTATTCTGGGCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACACGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTa                                      <  2:148156/249‑1 (MQ=255)
                                                                                                                                                                                                           gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTaa                                     >  1:537279/1‑249 (MQ=255)
                                                                                                                                                                                                              aaTTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGttt                                >  2:130806/1‑251 (MQ=255)
                                                                                                                                                                                                                       cggttaataaGTGGACAGTGATTACCCGGTGCGCGGTTAGTTTGTAAAGAAGAGTGTCAATATGGGTATAACGATTATTCTGGTGGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCTGGCAATGAAAACGATGGGTGTTAGCGAACTGCGGAGTGTCGATAGGCAGGCACACCTTGAGCCAGTCTCCCGCTGGGTCGCACATGGAGCTGGTGATATTATTGATAATATTAAAGTTTGTCCGACac                       <  2:389553/247‑2 (MQ=255)
                                                                                                                                                                                                                             ataaGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTg                 >  1:148497/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cACCCTGTTCGCGGTTATTTGGTCAAAAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATg  <  2:16731/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
TTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATG  >  NC_000913/4640666‑4641152

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