Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,410,034	G→A	10.1%	intergenic (+28/‑99)	rlmB → / → yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,034	0	G	A	10.1%	229.1 / 14.8	79	intergenic (+28/‑99)	rlmB/yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base A (5/3);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

TGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAC  >  NC_000913/4409888‑4410178                                                                                                                                                   |                                                                                                                                                  tGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAt                                                                                                                                                 <  4:500075/149‑1 (MQ=255)    cTCGCGAACATTGCGATGAGTTGATCAGCATCCCGATGGTTGGAAGCCTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCCGTGCGCCAGCGCAGCTAATTTCTCAGAAGTATGAAAAGCCATCCAGATTTg                                                                                                                                              <  4:576242/149‑1 (MQ=255)     tCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTgg                                                                                                                                             >  1:393389/1‑149 (MQ=255)      cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGa                                                                                                                                            <  3:357494/149‑1 (MQ=255)       gTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                           >  3:590886/1‑149 (MQ=255)       gTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                           >  3:73073/1‑149 (MQ=255)        tGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                                                                          >  2:182114/1‑149 (MQ=255)          aaCATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑t                                                                                                                                       <  3:150145/149‑1 (MQ=255)            cATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttt                                                                                                                                     >  1:666062/1‑149 (MQ=255)                gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttttttt                                                                                                                                 <  4:54209/149‑1 (MQ=255)                gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttttttt                                                                                                                                 <  2:466479/149‑1 (MQ=255)                 cGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tttttttt                                                                                                                                <  3:650023/149‑1 (MQ=255)                  gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGa                                                                                                                                            <  3:20535/137‑1 (MQ=255)                  gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGa                                                                                                                                            >  4:20535/1‑137 (MQ=255)                  gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttttttttt                                                                                                                               <  3:371475/149‑1 (MQ=255)                  gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttttttttt                                                                                                                               <  2:349268/149‑1 (MQ=255)                        ttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tttttttttttttta                                                                                                                         >  4:137555/1‑143 (MQ=255)                           aTCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGc                                                                                                                      >  2:697372/1‑149 (MQ=255)                            tCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCt                                                                                                                     <  1:182114/149‑1 (MQ=255)                               gCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGt                                                                                                                  >  3:331125/1‑149 (MQ=255)                                cATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTa                                                                                                                 <  2:125662/149‑1 (MQ=255)                                 aTCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTaa                                                                                                                <  3:509794/149‑1 (MQ=255)                                    ccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGAt                                                                                                             >  2:746243/1‑149 (MQ=255)                                       aTGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                                          <  1:85697/149‑1 (MQ=255)                                          gttGGAAGCGTTTCTTCCCTGAACGTTTCGTTTGCAACGGGAATTTGTTTATTTGAAGGGGTGCGCCAGCGCAGTTAATTTCTCAAAATAATAAAAAGCCATCCAGATTGGGTTGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                                                                       <  4:283356/147‑1 (MQ=255)                                            tGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                                     >  3:889265/1‑149 (MQ=255)                                            tGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                                     >  2:729382/1‑149 (MQ=255)                                             ggAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                                                    >  2:555869/1‑149 (MQ=255)                                             ggAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                                                    <  4:73073/149‑1 (MQ=255)                                                   gTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaa                                                                                              >  4:137034/1‑149 (MQ=255)                                                    tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                                                             >  4:698098/1‑149 (MQ=255)                                                          cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                                                       <  2:525869/149‑1 (MQ=255)                                                           ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCTATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                      <  4:640698/149‑1 (MQ=255)                                                            cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                     >  1:145725/1‑149 (MQ=255)                                                            cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttttcat                                                                                                                                <  3:786095/106‑4 (MQ=38)                                                            cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttttcat                                                                                                                                >  4:786095/1‑103 (MQ=38)                                                             tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                    <  1:160305/149‑1 (MQ=255)                                                               aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCATCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGAATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                                                                  <  2:461099/149‑1 (MQ=255)                                                               aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                                                                  <  3:491656/149‑1 (MQ=255)                                                               aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                                                                  <  1:746243/149‑1 (MQ=255)                                                                     tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                                            <  2:393389/149‑1 (MQ=255)                                                                      cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                                                           <  1:658844/149‑1 (MQ=255)                                                                       ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                                                           <  4:590886/149‑1 (MQ=255)                                                                         ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                                                        >  2:554131/1‑149 (MQ=255)                                                                          tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                                                                       <  3:137034/149‑1 (MQ=255)                                                                            cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                                                     <  3:698098/149‑1 (MQ=255)                                                                             gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttt                                                                                                                                    <  4:891618/85‑1 (MQ=38)                                                                             gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttt                                                                                                                                    >  3:891618/1‑85 (MQ=38)                                                                               ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                                                                  >  3:17967/1‑149 (MQ=255)                                                                                     tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                                                            >  1:295363/1‑149 (MQ=255)                                                                                       tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                                                          <  1:729382/149‑1 (MQ=255)                                                                                          ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGa                                                       >  3:251144/1‑149 (MQ=255)                                                                                            aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                                                     >  2:177393/1‑149 (MQ=255)                                                                                             tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                    <  3:866215/117‑1 (MQ=255)                                                                                             tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                    >  4:866215/1‑117 (MQ=255)                                                                                               tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                  <  1:555869/149‑1 (MQ=255)                                                                                                 aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                                                          <  2:389048/139‑1 (MQ=255)                                                                                                 aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                                                          >  1:389048/1‑139 (MQ=255)                                                                                                          gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATgg                                       >  2:131203/1‑149 (MQ=255)                                                                                                            tggcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGt                                                                                                                             <  1:630426/59‑1 (MQ=255)                                                                                                            tggcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGt                                                                                                                             >  2:630426/3‑61 (MQ=255)                                                                                                                  aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGtt                               <  2:295363/149‑1 (MQ=255)                                                                                                                  aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGtt                               <  4:331125/149‑1 (MQ=255)                                                                                                                    cTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                             >  2:333353/1‑149 (MQ=255)                                                                                                                        tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCt                         >  3:164362/1‑149 (MQ=255)                                                                                                                        tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGAAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCt                         <  4:889265/149‑1 (MQ=255)                                                                                                                         ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTa                        <  3:721150/149‑1 (MQ=255)                                                                                                                            tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCg                     >  4:654237/1‑149 (MQ=255)                                                                                                                            tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCg                     >  1:122884/1‑149 (MQ=255)                                                                                                                                      gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGc           >  4:926366/1‑149 (MQ=255)                                                                                                                                      gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGc           >  4:438135/1‑149 (MQ=255)                                                                                                                                      aaaaaACCATCCAAATCTGGATGGCTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTg            >  1:710928/2‑149 (MQ=255)                                                                                                                                         aaGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCaaa        <  3:438135/149‑1 (MQ=255)                                                                                                                                          aaccATCCAAATCTGGCTGG‑TTTTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATTTTTAAACTTCATCTATATGGAATAATATATGACATGGAATCGGTTGGCGCTAGCGACGGCGATGAAAAc       >  4:466811/3‑149 (MQ=255)                                                                                                                                             cATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGt    <  4:251144/149‑1 (MQ=255)                                                                                                                                             cATCCAAATCTTGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                        >  3:540468/1‑97 (MQ=255)                                                                                                                                             cATCCAAATCTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                        <  4:540468/97‑1 (MQ=255)                                                                                                                                               tCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAc  >  3:834354/1‑149 (MQ=255)                                                                                                                                               tCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAc  <  1:131203/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAC  >  NC_000913/4409888‑4410178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑