Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,146,736	A→G	6.1%	noncoding (148/169 nt)	psrD →	novel Novel sRNA, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,736	0	A	G	6.1%	169.8 / 3.1	66	noncoding (148/169 nt)	psrD	novel Novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (39/23);  new base G (2/2);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

TTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAG  >  NC_000913/1146588‑1146883                                                                                                                                                     |                                                                                                                                                    ttAGTCATATTTTTTTCCATCATATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCaa                                                                                                                                                     <  1:583086/149‑1 (MQ=255)      cATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                                                                                                                >  2:19773/1‑149 (MQ=255)        tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGCAAATAAAGCAGTAACCCGCAATGTTTGCGAATTATTGGCAAAAGGCAACCACAgg                                                                                                                                              >  1:767780/1‑149 (MQ=255)           ttttttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTTCACGAAAATAAAGCAGTAACACGCAATGTGTGCGAAGTATTGGCAAAAGGCAACCACAGGcgg                                                                                                                                           >  2:378530/1‑147 (MQ=255)                 ccATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCtttt                                                                                                                                     <  1:628294/149‑1 (MQ=255)                   aTCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCttttgc                                                                                                                                   <  1:944687/149‑3 (MQ=255)                   aTCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                                                                                                                   >  2:995238/1‑149 (MQ=255)                     cAGATATAGCGTATTGAAGATAGCGATTTTAAACTATGCGCTGCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCGGCCTTTTTCtt                                                                                                                                 <  1:4050/149‑1 (MQ=255)                      aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                                                                                                                <  1:24356/149‑1 (MQ=255)                       gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                               >  1:2865/1‑149 (MQ=255)                       gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                               >  2:244614/1‑149 (MQ=255)                          ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                                                                                                                            <  2:1102915/149‑1 (MQ=255)                          ataGCGTATTGATGATAGCCATTTTAAACTATGCGATTCGTTTTGCAGGTTGCTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGATTTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                                                                                                                            <  1:349062/149‑1 (MQ=255)                              cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                                                                                                                        <  2:915990/149‑1 (MQ=255)                              cGTATTGATGATAGCCATTTTAAACTAAGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                                                                                                                        <  2:1100055/149‑1 (MQ=255)                                     atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                                                                                                                 <  1:880334/149‑1 (MQ=255)                                        atAGCCATTTTAAACTATGCGCGTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTATCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTGCTTTGACTCTATGACGGTACaa                                                                                                              >  2:901814/1‑149 (MQ=255)                                            ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAgat                                                                                                          >  1:116047/1‑147 (MQ=255)                                            ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt                                                                                                          >  2:143529/1‑149 (MQ=255)                                             cATTTTAAACAATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa                                                                                                         >  2:226770/1‑149 (MQ=255)                                                       tctGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAACTAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgcgc                                                                                               <  2:512812/147‑1 (MQ=255)                                                           cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTa                                                                                           >  1:341261/1‑149 (MQ=255)                                                            ccTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAATCAGTAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGTCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTAt                                                                                          <  1:114239/148‑1 (MQ=255)                                                              ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTATTCTTAGACTCTATGACGTTACAAAGTTAATAAGCGCGCCCTATGc                                                                                        >  1:78199/1‑149 (MQ=255)                                                                          gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGGTGCCTTTTTCTTTGCCTCTATGACGTTACAAAGTTAAAATGCTCGCCCTATGAAAAAGGTAAAAt                                                                            >  1:815124/1‑149 (MQ=255)                                                                          gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTGTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                            >  1:296079/1‑149 (MQ=255)                                                                             gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAc                                                                         <  2:845534/149‑1 (MQ=255)                                                                              aTGTTTGTTATCAGCACTGAACGAAAATACAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAcc                                                                        <  1:244614/149‑1 (MQ=255)                                                                                gtttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAAGTTACAAAGTTAATATGAGCGCCCTATGCAAAAGGTAAAATTAcccc                                                                      >  2:474851/1‑149 (MQ=255)                                                                                 tttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCATAAGGTAAAATTACCCCt                                                                     >  1:578024/1‑149 (MQ=255)                                                                                    gttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGCCTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGGAAATTTACCCCTGaa                                                                  >  1:958472/1‑148 (MQ=255)                                                                                       aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                               <  2:135777/149‑1 (MQ=255)                                                                                       aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                               <  1:55057/149‑1 (MQ=255)                                                                                       aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                               >  1:490293/1‑149 (MQ=255)                                                                                       aTCAGCACTGAACCAAAATACAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGTCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTGACAAAGTTAATATGTGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                               <  2:1063580/149‑1 (MQ=255)                                                                                          aGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGa                                                            <  1:143529/149‑1 (MQ=255)                                                                                             aCTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATcc                                                         >  1:1064116/1‑149 (MQ=255)                                                                                               tGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                       <  1:19773/149‑1 (MQ=255)                                                                                                gAACGAAAAGAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAATGGCAACCACAGGCTGCATTTTTCTTTGACTCTATGACGTTACAAAGATATTATGCGCGCGCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                      >  2:877080/1‑149 (MQ=255)                                                                                                 aaCGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGtt                                                     >  1:639297/1‑149 (MQ=255)                                                                                                 aaCGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGtt                                                     >  1:395135/1‑149 (MQ=255)                                                                                                   cGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCg                                                   >  2:290640/1‑149 (MQ=255)                                                                                                     aaaaTAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTA‑CCCTGACTCTCGATCCGGTTCGTAc                                                >  1:418245/1‑149 (MQ=255)                                                                                                         tAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCATATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCCGATCGTACgga                                             >  2:1116823/1‑148 (MQ=255)                                                                                                           aaGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCGTTGACTCTATGCCGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCGGGTTCGTACGGCTc                                           >  2:727165/1‑149 (MQ=255)                                                                                                                gTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGACTTTTTCTTTGACTCTATGACTTTACAAAGTTAATATGCTCGCCCTATGCAAAAGGTAAAATTACCACTGACTCTCGATACGGTTCGTACAGCTCaaaaa                                      >  1:815871/1‑149 (MQ=255)                                                                                                                 tAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTAAATATGCGCGCCCTATGCAAAACGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                                     <  2:116047/149‑1 (MQ=255)                                                                                                                    cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTATATGACGTTACAAAGTTAAAATGCGCGCCCTATTCAAAAGGAAAAATTACCACTCACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                                  >  2:774490/1‑149 (MQ=255)                                                                                                                     atgCATTGTGTTCTAATTATTGGCAAAAGGCAAGCACAGGGTTTCTTTTTCTTTGACTCTATGACGTTACTAAGTTAATATTCGCGACATATGCATAAGGTAATATTAGGCGTGACTCTCGATGAGGTTCGTTGGGCTCAAAAACGcat                                 <  2:465540/147‑3 (MQ=255)                                                                                                                           tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTa                           <  1:995238/149‑1 (MQ=255)                                                                                                                               tgCGAATTATTGTCTAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAg                       <  2:296079/149‑1 (MQ=255)                                                                                                                                gCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAgg                      <  1:304904/149‑1 (MQ=255)                                                                                                                                 cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGCCTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAggg                     >  1:594643/1‑149 (MQ=255)                                                                                                                                 cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCGTGTTTACCAggg                     >  1:768193/1‑149 (MQ=255)                                                                                                                                    attattGGCAAACGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTaa                  >  2:978206/1‑148 (MQ=255)                                                                                                                                         aagaaaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAATTTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata             >  1:174692/5‑149 (MQ=255)                                                                                                                                         aagaaaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata             >  2:497923/5‑149 (MQ=255)                                                                                                                                           gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc           >  2:434473/1‑149 (MQ=255)                                                                                                                                           gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc           >  1:580761/1‑149 (MQ=255)                                                                                                                                           gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATAACCAGGGTATCTATAcc           >  1:543321/1‑149 (MQ=255)                                                                                                                                           gCAAAAGGCAACCACAGGCTGCCGTTTTCATTGGCTCTATGACGTTACACAGTTAATATGCGCCCCCTATGCAAAAGGTAAAATTACCCGTGACTCACGCTACGGTGCGTACGGCTCAAAAACCCGATTCTTACCAGGGTATCTATAcc           >  2:1060957/1‑149 (MQ=255)                                                                                                                                            cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGCTCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAccc          >  2:1009255/1‑149 (MQ=255)                                                                                                                                            aaaaaGGCAAACACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAATTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAccc          <  2:641378/148‑1 (MQ=255)                                                                                                                                             aaaaGGCAACCACAGTCTTCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCGTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCTAAAACGCCTTGATTACCAGGGTATCTATAcccc         <  2:78199/149‑1 (MQ=255)                                                                                                                                             aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAACCGCCTTGATTACCAGGGTATCTATAcccc         >  2:729649/1‑149 (MQ=255)                                                                                                                                               aaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGAATCTCGATCCGGTTAGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATCCCCCTg       >  1:541635/1‑149 (MQ=255)                                                                                                                                                  gCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATc    >  1:537698/1‑149 (MQ=255)                                                                                                                                                    aaCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAg  >  1:18142/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAG  >  NC_000913/1146588‑1146883

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