Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,824,731 T→G 22.1% V73G (GTG→GGG)  trmH → tRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,824,7310TG22.1% 142.7 / 2.7 70V73G (GTG→GGG) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (11/39);  new base G (16/3);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCC  >  NC_000913/3824507‑3824966
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
aTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGAAGGAAAAGCTGGGTACAG‑‑GGGAAAAAACACCGc                                                                                                                                                                                                                                      >  1:435017/1‑238 (MQ=255)
     ttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGTGTAAAAGCTGGGGACCA‑‑GGGAAAAACCAACGCAcc                                                                                                                                                                                                                                   >  1:300706/1‑236 (MQ=255)
      tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTCGTAAAAGCCGGGTACAG‑‑GGGGAAAAACACCCCACCATTg                                                                                                                                                                                                                               >  1:505385/1‑239 (MQ=255)
               aaaccGTTTTGCCCGTATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTTTGCATGGAGCAGTTCCCAAAACCTCATAACGTTTCTGCGATGATTCGTACCGCAGATGCCGTTGGCGTACACGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACCCACCGCACCATTGGCGATGCCg                                                                                                                                                                                                                      <  2:897345/235‑1 (MQ=255)
                    gtttggCACGCATCTGGGAAATGCCCGCCGGGCGTCAGCCTGATCTGCCCGTCTGATTGGCGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATCCGTCCCGCAGATGCCGTTGGCGTACAGGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCt                                                                                                                                                                                                                 <  2:680441/234‑1 (MQ=255)
                      tacgcgcgagtgtgggAAATTCTGGCAAGGCGGCCGCCTGATCTTGCCGTCTGCATGGGAGAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGGCGTTGGCGTACATGAAGTTCACGCCGTCCGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGTGAAAACACACCGCACCATTGGCGATGCCGTCGCt                                                                                                                                                                                                                 <  2:622313/224‑1 (MQ=255)
                      tATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACGTCCGGGGCACGG‑‑GGAAAAACACCCCGCACAATGGGGGAtccgggcgtca                                                                                                                                                                                                                >  1:558235/1‑229 (MQ=255)
                           gcGCCTCTGCGAAAGGCTCCCCGGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTTTTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCa                                                                                                                                                                                                          <  2:435017/238‑1 (MQ=255)
                           ccGCACCGGCAAAATTCCCGCCGGGGGGAAGCCTGATGTGCCCGTTTGGATGGGGCAGGTCCACAACCCTCATAACGTTTGTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCa                                                                                                                                                                                                          <  2:834777/238‑1 (MQ=255)
                              cATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGGACAG‑‑GGGAAAAACAACCGCCCCATTGGCGAGGCCGGCGCTCATCCCAAAg                                                                                                                                                                                                       >  1:523301/1‑239 (MQ=255)
                               aTCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAACTGCGGGTGGGACAAGCTGGGTACAG‑‑GGGAAAAAACCCCCACCCACTGGCGAAGACGTCGCTCacaccaaaag                                                                                                                                                                                                      >  2:19619/1‑230 (MQ=255)
                                 ttGCGAAATTCCCGCCAGGCGGCAGCCTTATCTGACCTTCTGGATGGAGCAGGTCCACAAACCTCATTACGTTTCTGCGATTATGCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                                                                                                                                                                                                    <  2:588905/238‑1 (MQ=255)
                                 ccgagaatggCTCGCCGGGCGGCGGGTTGATTCGACCGTCGGCATGGAGCAGGGCCACAAACCGAATACCGTTTCTGCGGTTATTCGTACCGCAGATGCCGTTTGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                                                                                                                                                                                                     <  2:183429/229‑1 (MQ=255)
                                 cTGGGAAATGGCCGCCAGGGGGAAGCCGTATCTGACCGTTTGCATGGAGCAGGCCCACAACCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                                                                                                                                                                                                    <  2:523301/239‑1 (MQ=255)
                                 cTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTAAAG‑‑GGGTAAAAAACCCGCACCATTGGGGATGACGTCGCCCATCTCCAAAGcc                                                                                                                                                                                                    >  1:177266/1‑239 (MQ=255)
                                         tGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGGTGCCGTTCGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGGGTACAAGCTCGGTCCCA‑‑GCGAAAAACAACCGCCCCCTTGGGGCTGCCGCGCCCCATCTCaaaaagcaaggcgtg                                                                                                                                                                                            >  2:340221/1‑227 (MQ=255)
                                              gCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GGGGAAACACACCGCACCATTGGCGCTTCCGTCGCTCATCTCAAAGGCCCGGGCaatcaaaa                                                                                                                                                                                       >  1:817889/1‑232 (MQ=255)
                                                          ccTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGAACCATGGCTTCGGCAGCGGCGGGTCGTAACAGCGGGGTTCAG‑‑GGGGAAAAACACCCCCCCCATTGTGATGCCCTCGCCCCCCTCAACAGCCAGGGGCCTCCGATTCTGTCAAcccc                                                                                                                                                                         >  2:529694/1‑238 (MQ=255)
                                                              aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCtt                                                                                                                                                                     <  1:909815/239‑1 (MQ=255)
                                                                    ccttctgctggggagaGGTCCACAACCTCCATAACTTTTCTGGGATTATTGGTCCCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGa                                                                                                                                                                <  2:89323/223‑1 (MQ=255)
                                                                        cttcctGGGGCCGGTCCCCAAACCCTATAACGTTTCGGGGATTTTTCGTCCCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACg                                                                                                                                                           <  2:817889/234‑1 (MQ=255)
                                                                        ctcccttggtctggTTCACAAGCCTTATAACGTTTCTGCGATTATTCGGACCGCGGATGCCTTTGGCGTCCATGAAGTTCAGGCCGTCTGGGCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGGACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACg                                                                                                                                                           <  2:847786/227‑1 (MQ=255)
                                                                        cTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGGGTAACACCGGGGGACCG‑‑GGGAAAACACACCGCCCCCCAGGCGAGGGCGTGTCTCCACTCAAAGGCCAGGGACCGCAGAATTCTTCAACCCCCCTTTCTggttaa                                                                                                                                                            >  2:154430/1‑233 (MQ=255)
                                                                        cTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGTCAG‑‑GTGAAAAAACACCGCACCCTTGGCGATGCCGTGCCTCCACTCAAAAGCCAGGGCAAGCAGATTCTGGCAACACACATTTCTGTTaaa                                                                                                                                                            >  1:14805/1‑237 (MQ=255)
                                                                                  cAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACCAACCCCACCATTGGCGAGTCCGCGGCTCATCTCAAAGGCCAGGGCATGCCAATTCTGGCAACCCATCTTTATGGTAACGCTGTCGAtt                                                                                                                                                >  1:42770/1‑238 (MQ=255)
                                                                                          caacccctaaaaCGTTTCGCCGATTATTCGTACCGCGGTAGCCGGTTGGGTACAGGAAGTTCACGCGGTCGGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGaaa                                                                                                                                       <  2:794056/229‑1 (MQ=255)
                                                                                          caAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCGGGGAACAG‑‑GGGAAAACACCCCCCACCCTAGGCGTTGCCGCCGCGCATTTCAAAGGCCCGGGCGTACAGAATTAGGGACACCATCTTTTTGCTAACGCCGTCGCTGTcccgcaaa                                                                                                                                       >  2:884071/1‑233 (MQ=255)
                                                                                               cTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGCAAG‑‑GGGAAAACACACCCCACCATTGGCGATGCCGTCGCTCATCTCCAAAGCCAGGGCATGCAGATTCTCGCAAACCACCCTTCTTATAACCCCGTGGAGTTTCGCGGAAATGAt                                                                                                                                  >  2:481419/1‑239 (MQ=255)
                                                                                                   tAACTTTTTTGCCATTTTTCGGTCCGCGGATGCCGTTGGCGTACATGAAGTTCCCGCCGTCTGGCCCGGGAGCCGGATGCGGACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTaca                                                                                                                              <  2:119051/239‑1 (MQ=255)
                                                                                                          tccgCGATTACTCGTACCGCATATGCCGTTGGGGTACATGAAGTTCACGCCGGCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGccc                                                                                                                       <  1:612198/236‑1 (MQ=255)
                                                                                                            tGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACCGCTGGGTACAG‑‑GGGAAAACACACCGCACCACTAGCGATGCCGCCGCTCATCTCCAAAGCCAGGGCGCGCCGCTTCTGGCCACCCACCTTTCTGTTTACACTGTTCGATTACCGCAAAAtatataaaccccccccc                                                                                                                      >  2:384830/1‑223 (MQ=255)
                                                                                                                 ttattCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCCGGGTGCAG‑‑GGGAAAACACACCCCCCCCCTGGGGATGCCCTCGCTCCACTCAAAGGCCAGCGCCTGCAAATTCGTGACACCCCACTTTTTCATAACGCTTTCCGATTCTCCCAAATTGTTTaccaccccccccccac                                                                                                                 >  2:487326/1‑224 (MQ=255)
                                                                                                                   attCGTAACGCAGATGCCTTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCAt                                                                                                              <  1:529694/239‑1 (MQ=255)
                                                                                                                           cGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGGGATGCCGTCGCTCATCTCACAGGCCAGGGGCTGCGAATTCTGGCAACCCATCTTTCTGTAAAACCGTGCGATTTCCCCGCAATTTATTACACACcccccaaccgacttttggggg                                                                                                      >  2:328856/1‑222 (MQ=255)
                                                                                                                                        ggCGTACATGAAGTTCCCGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAAc                                                                                         <  1:481419/239‑1 (MQ=255)
                                                                                                                                        ggCGTACATGAAGTTCCCGCCGTCTGGCCTGGTAGCCGCATGCGCACCAGGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAAc                                                                                         <  1:791276/239‑1 (MQ=255)
                                                                                                                                          cGTACATGAAGGTCACGGCGGCTGGCCTGGTAGCCGCATGCGGACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACCCCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACg                                                                                        <  1:411666/238‑1 (MQ=255)
                                                                                                                                                tGAAGTTCCCGGCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCa                                                                                <  1:19619/239‑1 (MQ=255)
                                                                                                                                                   atttCCCGCCGTCTGGCCGGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGc                                                                             <  1:566624/237‑1 (MQ=255)
                                                                                                                                                    gtccACGCCGGCTGGCCTGTTAGCCGCATGCGAACCAGGGCTTCGGCGGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGa                                                                                                         <  2:851264/208‑1 (MQ=255)
                                                                                                                                                      tccccccctttGGCCGGGTAGCCGGAGGCGGACCAGGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAgg                                                                          <  2:847275/229‑1 (MQ=255)
                                                                                                                                                         cGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCGGGGTACAG‑‑GGGGAAACACACCCCCCCATTGGCGAATCCGTCGCTCCTCTTAAAAGCCAGGGCCAGCAGATTTTGGCAACCCCTCTTTTTGTTAACCCTGACGATTTCCGCGAAATTGATTTTACTCTCCCCACCCGCACTTTTTTGGGGCGGTGGAAAAAggggaaaaagccgagga                                                                        >  2:692898/1‑225 (MQ=255)
                                                                                                                                                               cTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GGGAAAACACACCCCACCCTTGGCGATGCCGTCGCTCATCTCAAAAGCCCGGGCATCCCGATTCTGGCAACCCCCCTTTTTCTTAACGCTGGTGCGTTCCCCGCAACTGGTTACACTCCCCCGACCTTCATTTTTCTGGGGCACGAGAAAAACGGCATACCGCCGGGAGCAATg                                                                 >  2:195580/1‑238 (MQ=255)
                                                                                                                                                                  gCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGaa                                                                                                                                        <  1:422672/166‑1 (MQ=255)
                                                                                                                                                                     tgtgcgccccagcgcgaCCCTGGCTTCGGCAGCGGCGGGTAGTAAGTGCGGGGTACCG‑‑GTTAAAACCCACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTgg                                                          <  2:888807/223‑1 (MQ=255)
                                                                                                                                                                     tGGTAGCCGCATGCGCCCCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTgg                                                          <  1:556948/239‑1 (MQ=255)
                                                                                                                                                                     tGGTAGCCGCATGCGACCAATGGCTTCGGCAGCGGCGGGTAGTAACAGGTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTgg                                                          <  1:884071/239‑1 (MQ=255)
                                                                                                                                                                     ggttggCCGCGTGCGCACCATGGCTTTGGGCGCGGGGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACCCCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCAGCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTgg                                                          <  2:602548/234‑1 (MQ=255)
                                                                                                                                                                           cgcgatacgcgCCTTGGCTTCGGCAGGGGCGGGTAGTAACCGCGGGGTACAG‑‑GGGAAAACACGCCGCCACATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGAt                                                     <  2:177266/231‑1 (MQ=255)
                                                                                                                                                                           ccccgggtgcgcaaTGGCTTCGGGAGGGGCGGGTAGTAGCAGCTGGGTACAG‑‑GGGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTTCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATc                                                    <  2:558235/226‑1 (MQ=255)
                                                                                                                                                                           ccGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCCCGGTTCAC‑‑GTGAAAACACACCGCACCCTTTGCGGTGCCCTCCCTCATCCCAAAGGCCAAGGCCTACCGATTCTTGCTACCCATCTTTCTCATAACGCTGTCGATTTACCCGAAATTTACTACACCCGCCCCACCTGCATCTTTCTGGGACACTGGAAAACAGCGATCCACCCAGAAACAATTGGCCCTgcgggag                                                     >  2:412554/1‑233 (MQ=255)
                                                                                                                                                                            cGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATc                                                    <  1:381917/238‑1 (MQ=255)
                                                                                                                                                                                  gccccccaggCTCCGGCAGGGGCGGGTAGTAACAGCCGGGGACAG‑‑GTTAAAACCCACCGCACCATTGGCGATGCCGTCGCCCAACTAAAAGGCCAGGGCATGCAGATTCTGGCATCCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGa                                                <  1:340221/228‑1 (MQ=255)
                                                                                                                                                                                   gccccATGGCTTCGGCAGGGGCGGGTTGTAACAGGTGGGTACAG‑‑GGGAAAACACACCGCCCCATTGGCGATGCCGGCCCTCATCCCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCATGGCGGATCAGGAcat                                             <  2:902502/235‑1 (MQ=255)
                                                                                                                                                                                     aCCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACCCACCGCCCCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGAcatcat                                          <  1:66221/239‑1 (MQ=255)
                                                                                                                                                                                      ccATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGGTAAAGCTGTCGATTTCCGCGAAAATTATTAAACTCGCCCGACCTGCATTTTTGTGGGAACGGGGAAAACAGGCCTAACGCCCGAAGCACTTTGCCCTgcgggacagagaaacaaa                                          >  2:585696/1‑222 (MQ=255)
                                                                                                                                                                                          ggCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCa                                                                                                                                                                                                          <  2:660758/80‑1 (MQ=255)
                                                                                                                                                                                          ggCTTCGGCAGCGGCGGGTAGGAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCa                                                                                                                                                                                                          >  1:660758/1‑80 (MQ=255)
                                                                                                                                                                                            cTTCGGCAGGGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCg                                   <  1:277497/239‑1 (MQ=255)
                                                                                                                                                                                             ttCGGCAGCGGCGGGTAGTAACAGCTGGGTTCAG‑‑GTAAAAACCCACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCg                                   <  1:692898/238‑1 (MQ=255)
                                                                                                                                                                                             ttCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCg                                   <  1:293721/238‑1 (MQ=255)
                                                                                                                                                                                                 gccggggcggGTGGTAACCGCTTGGTACCG‑‑GTGAAAACACACCGCCCCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCgatgat                              <  2:300706/234‑1 (MQ=255)
                                                                                                                                                                                                              tgtacacgCTGGGTCCAG‑‑GTGAAAACCACCCGGCCCTTTGGCGATGCCGTCGCTCAGCGCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTGCGGATAACGCTGTCGAGTGCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGAGGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCa                  <  2:405285/232‑1 (MQ=255)
                                                                                                                                                                                                              gtaACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCa                  <  1:384830/238‑1 (MQ=255)
                                                                                                                                                                                                              gtaACAGCTGGGGACAG‑‑GTGAAAACACACCGCTCCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCa                  <  1:74814/238‑1 (MQ=255)
                                                                                                                                                                                                                      cTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGCCAGGAGAAAACGGGCATAACGCAGGAAGCATTGGCCCTGCCGGATCAGGACATCCTCTTTCCGATGACCGGCCTGGTGCGGTccatga          >  1:548275/1‑234 (MQ=255)
                                                                                                                                                                                                                      cTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCCCCCCGGAAGCATTGGCCCTGGCGGGTCAGGAAAACCGCATTCCGATTATCGGCCTGGGGCCAGCGCTGa          >  2:74138/1‑239 (MQ=255)
                                                                                                                                                                                                                         ggTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGCCAGTAGAAAACGGGCATAACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATg       >  1:851858/1‑239 (MQ=255)
                                                                                                                                                                                                                         ggTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAAAGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAAAGGGCATCACGCAGAAACACTTGGCCCCGGCGGATCAGAACATAACAATTCCGATAACCGGCATgggccatgcggtga          >  2:728082/1‑224 (MQ=255)
                                                                                                                                                                                                                              aG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGGTTcc  >  1:561443/1‑239 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAG‑‑GTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCC  >  NC_000913/3824507‑3824966

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