Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,117,707 G→A G301S (GGT→AGT)  trhO → tRNA U34 hydroxylase TrhO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,117,7070GA97.9% 146.4 / ‑4.5 48G301S (GGT→AGT) trhOtRNA U34 hydroxylase TrhO
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/23);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

ACAAAGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAGGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAA  >  NC_000913/1117462‑1117952
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
acaaAGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAgttgt                                                                                                                                                                                                                                                  >  1:54716/1‑251 (MQ=255)
        tGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTg                                                                                                                                                                                                                                           >  1:412280/1‑250 (MQ=255)
                cGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTtgct                                                                                                                                                                                                                                  <  1:318327/251‑1 (MQ=255)
                                   caCGCCCGTAAGGCGCGCGAGCAGGGCTTGTCGGTGCGTTTTATGTGTAAAAATTTTGCTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCCCCAGTGCGGCGCGCCGTGTTATAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGttccc                                                                                                                                                                                                               <  1:496676/250‑4 (MQ=255)
                                     cGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGc                                                                                                                                                                                                             <  1:99175/251‑1 (MQ=255)
                                        ccGTAAGGCGCGCGAGCAGGGCTTGCAGGTGCGGTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGcc                                                                                                                                                                                                            >  1:402078/1‑249 (MQ=255)
                                                           ggCTGGCCGGTGCGTTTTATTGGCAAGAAGTTTGTTTTGGAGGAGCGGCTGGGCGAACGCATAGCGGATGATATTATCGCGCATTGCCACCAGTGCGGCGCGCCTTGGGTCAGCCATACCAACGGTAAAAAGGATGGCGGCCACCTGCGGTTTATTCAGTGTCCAGTATGCGGGGAAAAATACAAAAGTTGTGGTAGTGATACTTGCCGCGAAGAAAGCGCGTGACCGCCagaagatgagctactgcgtcg                                                                                                                                                                                       <  2:436290/251‑16 (MQ=255)
                                                               tGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCggg                                                                                                                                                                                     <  1:368562/249‑1 (MQ=255)
                                                                   ggggCGTTTTATGGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCAGTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGAc                                                                                                                                                                               <  2:512172/248‑1 (MQ=255)
                                                                   ggTGTGTTTTATTGGGAAAAATTTTATTTTTGAAGGGAGGAGGGGAGAACGTATATCGGATGAGATTAGCGCGCAGTGCCACCAGTGCGGCGCGCCGTGCGACAGCCAGACCAACTGTAAAAATGATGGCTGCCACCTGTTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAGGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGAc                                                                                                                                                                               <  2:266078/251‑1 (MQ=255)
                                                                       ggTTTTATTGGCAAAAATTTTGTTGTTGACGAGCGGATGGGCGAACGTATATCGGAGTAGATTATCGCGCTATGCAACCAGGGGGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGGGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGAGTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGa                                                                                                                                                                           <  2:281786/250‑1 (MQ=255)
                                                                           tttttgtcAAAAATTTTTGTTTGGACGTGCGGATGGGCGAACGTAAATCTGAGGAGATTATCGCGCATGTCCACCAGTGCGGCTCGCCGTTGGACAGCCATACCAACTGAAAAAATGATGGCGGCCACCTGCTTTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGaaaa                                                                                                                                                                        <  2:402078/243‑1 (MQ=255)
                                                                                ggCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCa                                                                                                                                                                   <  1:502297/250‑1 (MQ=255)
                                                                                 gCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCCGTGTCCAGTATGCGCGGAAAAATAAAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGACGGGCAGGGCGTGAAAATGGCaa                                                                                                                                                                  >  2:661951/1‑250 (MQ=255)
                                                                                          ttGTTTTTGACGAGCGGTGGGGCGAACGTATATCGGAGGAGATTATCGCGCATTGCAACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCtt                                                                                                                                                         <  2:107080/250‑1 (MQ=255)
                                                                                                     gAGCGGATGGGCGAACGTATCTCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAACATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGCACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATCAGATATTTAATAAGTCTCg                                                                                                                                             >  2:654375/1‑251 (MQ=255)
                                                                                                         ggATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGCGAAAATGGCAATAAGATCTTTAATAAGTCTCGTgg                                                                                                                                          >  2:523803/1‑250 (MQ=255)
                                                                                                          gATGGGCCAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCTTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGAc                                                                                                                                        <  1:152871/251‑1 (MQ=255)
                                                                                                              ggCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCCGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGAACTTTAAAAAGTCTCGTGGACGTCt                                                                                                                                    >  2:597353/1‑251 (MQ=255)
                                                                                                                cgaatgtatttgggATGATATTATCGCTCATTGCCACCAGTCCGGGGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGGGAGATTTGCGGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGa                                                                                                                                  <  2:412280/239‑1 (MQ=255)
                                                                                                                 gAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGa                                                                                                                                  >  1:199410/1‑250 (MQ=255)
                                                                                                                   aCGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGCCGTGAACATGGCAATCAGATCTTTAATAAGTCTCGTGGACGTCTGAAt                                                                                                                                >  1:469990/1‑250 (MQ=255)
                                                                                                                     gTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTATTAAGTCTCGTGGACGTCTGAATac                                                                                                                              >  2:467531/1‑250 (MQ=255)
                                                                                                                      tataTCGGATGAGATTATCGCGCATTGCCAACAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAGAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGACGAACAGCGACGCCGACGGGCAGGACGTGAACATGGCAATAAGATCTTTAATAAGTCTCGTGGCCGGCTGAATaaa                                                                                                                             >  2:652517/1‑248 (MQ=255)
                                                                                                                         atCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATacaaca                                                                                                                          <  1:440827/250‑1 (MQ=255)
                                                                                                                            ggATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACtgt                                                                                                                      <  1:652517/251‑1 (MQ=255)
                                                                                                                              aTGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTta                                                                                                                                                       >  1:267493/1‑216 (MQ=255)
                                                                                                                              aTGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTta                                                                                                                                                       <  2:267493/216‑1 (MQ=255)
                                                                                                                                                 ccaccaGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGaat                                                                                                 >  2:517744/1‑251 (MQ=255)
                                                                                                                                                  ccccaGGGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGGTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGaata                                                                                                <  2:54716/249‑1 (MQ=255)
                                                                                                                                                     caGTGCGGCGCGCAGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGCCGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGCCGTATGAATACAACACTGTGCATTACTGATCCAAAAGaaaaaat                                                                                             >  2:468388/1‑251 (MQ=255)
                                                                                                                                                      aGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCACAGAATAAatat                                                                                           >  1:152277/1‑251 (MQ=255)
                                                                                                                                                             cgcgcCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGc                                                                                     >  1:621542/1‑251 (MQ=255)
                                                                                                                                                                gcCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATACAACAGAATAAATATCATTGCCgg                                                                                  >  2:205674/1‑251 (MQ=255)
                                                                                                                                                                              aTACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCg                                                                    >  1:497267/1‑251 (MQ=255)
                                                                                                                                                                                       gTAAAAATGATGGCTGCCACCTGCTGTGTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACGGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACAt                                                             <  1:661951/249‑1 (MQ=255)
                                                                                                                                                                                            aaTGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCAGGTTGCGTGCCATCCGGCAACATTTCACGc                                                      >  2:122226/1‑251 (MQ=255)
                                                                                                                                                                                                           ctgctgTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTAc                                       >  2:323412/1‑251 (MQ=255)
                                                                                                                                                                                                             gctgTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTAc                                       <  1:69867/249‑1 (MQ=255)
                                                                                                                                                                                                               tgTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACAt                                                             <  1:551163/225‑1 (MQ=255)
                                                                                                                                                                                                               tgTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACAt                                                             >  2:551163/1‑225 (MQ=255)
                                                                                                                                                                                                                     ttCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACt                             >  1:480660/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGTAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCa                   <  1:63473/251‑1 (MQ=255)
                                                                                                                                                                                                                               aGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGcc                                                                                    <  2:568938/186‑1 (MQ=255)
                                                                                                                                                                                                                               aGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGcc                                                                                    >  1:568938/1‑186 (MQ=255)
                                                                                                                                                                                                                                              ctACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAACGCGCGTTACCGACAGAAGAACAGCGACGCCGTCGGGCAGGACGTGCAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATCCAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGAGTGCCATCCGGCAACATTTAACGCTTAATTATGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTg      >  2:68616/2‑249 (MQ=255)
                                                                                                                                                                                                                                                aCAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTgaga   <  1:523803/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aCAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGaa  <  1:205674/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ACAAAGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAGGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAA  >  NC_000913/1117462‑1117952

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