Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,146,740	C→G	6.1%	noncoding (152/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,740	0	C	G	6.1%	182.7 / 4.9	66	noncoding (152/169 nt)	psrD	novel Novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (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 = 3.56e-01

CATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAG  >  NC_000913/1146593‑1146883                                                                                                                                                    |                                                                                                                                                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                                                               <  1:55057/149‑1 (MQ=255)                                                                                  aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                               <  2:135777/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:395135/1‑149 (MQ=255)                                                                                            aaCGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGtt                                                     >  1:639297/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)                                                                                                                                                    |                                                                                                                                                CATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAG  >  NC_000913/1146593‑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: ‑