Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,398,552	G→A	5.4%	A24V (GCT→GTT)	uspE ←	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,398,552	0	G	A	5.4%	235.0 / 2.2	74	A24V (GCT→GTT)	uspE	stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base A (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGA  >  NC_000913/1398407‑1398699                                                                                                                                                  |                                                                                                                                                    gATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTTGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTc                                                                                                                                                  >  6:64715/1‑149 (MQ=255)  aTCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcg                                                                                                                                                 >  3:85364/1‑149 (MQ=255)     cAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGTGTGGTCATTTCTTAAGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              <  1:337908/149‑1 (MQ=255)      aGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              >  3:31853/1‑148 (MQ=255)      aGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              >  7:143325/1‑148 (MQ=255)       ggCTGTACGCTGGCTGATGCCGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccgc                                                                                                                                             >  3:262090/1‑148 (MQ=255)       ggCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccgc                                                                                                                                             >  2:314956/1‑148 (MQ=255)              cGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtgg                                                                                                                                     <  2:219018/149‑1 (MQ=255)                 tggctgATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggttg                                                                                                                                  >  5:339189/1‑149 (MQ=255)                   gctgATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggttg                                                                                                                                  <  6:43579/147‑1 (MQ=255)                          aCGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTg                                                                                                                         <  1:282643/149‑1 (MQ=255)                             cccTGACGCATAGCGGTACGTTCGTCCTGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTTCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGtt                                                                                                                      <  6:161527/149‑1 (MQ=255)                               cTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAg                                                                                                                    <  4:31853/149‑1 (MQ=255)                                      aTAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGa                                                                                                              >  7:145222/1‑148 (MQ=255)                                      aTAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGAt                                                                                                             <  1:314956/149‑1 (MQ=255)                                           ggTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATaacaa                                                                                                        >  8:175920/1‑149 (MQ=255)                                            gTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATaacaac                                                                                                       <  6:72787/149‑1 (MQ=255)                                                gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACg                                                                                                      <  6:235321/146‑1 (MQ=255)                                                gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACg                                                                                                      >  5:235321/1‑146 (MQ=255)                                                gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAg                                                                                                    >  2:337183/1‑148 (MQ=255)                                                    gTCCGGGTAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGAAAATAAACAGCTCGCCGCAATGCTCGATGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                                                               >  2:222895/1‑149 (MQ=255)                                                      ccGGGGAGAGCAGGGTGGTCATTTCGTATTAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATTTAAATAAACAGCTCGCCGCAATGCTGGTGGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             <  1:240339/149‑1 (MQ=255)                                                      ccGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             >  5:41265/1‑149 (MQ=255)                                                       cGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             <  4:262090/148‑1 (MQ=255)                                                       cGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCt                                                                                            <  8:143325/149‑1 (MQ=255)                                                         gggAGAGCAGGGTGGTCATTTCGTATTATAAGTCATAGATAGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATTATGGTTGGTCTTCCTGGTTAGGATAGATAACAACGAGAATGTTCTGa                                                                                          >  3:261020/1‑149 (MQ=255)                                                           gagaGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATa                                                                                        <  4:85364/149‑1 (MQ=255)                                                             gagCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACa                                                                                      >  8:348683/1‑149 (MQ=255)                                                               gCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                                                                    <  6:301744/149‑1 (MQ=255)                                                                cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                                                                   >  5:627/1‑133 (MQ=255)                                                                cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                                                                   <  6:627/133‑1 (MQ=255)                                                                cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAg                                                                                   <  3:266568/149‑1 (MQ=255)                                                                cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAg                                                                                   <  6:41265/149‑1 (MQ=255)                                                                     tggtCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATa                                                                              <  7:300558/149‑1 (MQ=255)                                                                        tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAg                                                                           >  1:162024/1‑149 (MQ=255)                                                                         cATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAgg                                                                          <  6:339189/149‑1 (MQ=255)                                                                          aTTTCGTATGAGAAGTCATAGATCGGCAAAAATGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAgg                                                                          <  4:261020/148‑1 (MQ=255)                                                                          aTTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAggg                                                                         >  8:361491/1‑149 (MQ=255)                                                                           tttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGt                                                                        >  7:131376/1‑149 (MQ=255)                                                                           tttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTAGACGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGt                                                                        >  5:176661/1‑149 (MQ=255)                                                                                tATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCt                                                                   <  7:175920/149‑1 (MQ=255)                                                                                tATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCt                                                                   >  4:180619/1‑149 (MQ=255)                                                                                  tGAGAAGTCATAGATCGGCAAAAAGTCCTTAATTTTGCCACCAATCCGTTTATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGCTAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTa                                                                 <  1:222895/149‑1 (MQ=255)                                                                                    agaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTac                                                                >  1:370280/1‑148 (MQ=255)                                                                                     gaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGtt                                                                                              <  5:239922/117‑1 (MQ=255)                                                                                     gaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGtt                                                                                              >  6:239922/1‑117 (MQ=255)                                                                                         tCATAGATCGGCAAAAAGGCTTTAATTTTGCCCCCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACaacaac                                                          >  1:106653/1‑149 (MQ=255)                                                                                          cATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACaacaac                                                          <  1:337183/148‑1 (MQ=255)                                                                                           aTAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTg                                                        >  4:321525/1‑149 (MQ=255)                                                                                               aTCGGCAAAAAGGCTTTAATTTTGCCCCCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCaa                                                    >  6:247213/1‑149 (MQ=255)                                                                                                      aaaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                                             >  4:255869/1‑149 (MQ=255)                                                                                                      aaaaGGCTTTAATTTTGCCACCAAACCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATACCAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                                             >  6:102416/1‑149 (MQ=255)                                                                                                       aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGttt                                            >  8:300593/1‑149 (MQ=255)                                                                                                       aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGttt                                            <  3:255869/149‑1 (MQ=255)                                                                                                          ggCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         <  3:105392/149‑1 (MQ=255)                                                                                                             tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAAt                                       >  2:295197/1‑148 (MQ=255)                                                                                                               tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGATGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACACCTGTCAACGCAGTATGTAATTa                                     >  4:47907/1‑148 (MQ=255)                                                                                                                aaTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaa                                    >  8:348264/1‑148 (MQ=255)                                                                                                                 aTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaaa                                   <  7:161173/148‑1 (MQ=255)                                                                                                                         accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  3:34189/1‑149 (MQ=255)                                                                                                                          ccgaTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         <  3:9868/130‑1 (MQ=255)                                                                                                                          ccgaTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         >  4:9868/4‑133 (MQ=255)                                                                                                                          ccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  8:85014/1‑148 (MQ=255)                                                                                                                          ccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  1:154404/1‑148 (MQ=255)                                                                                                                                 gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  4:34189/149‑1 (MQ=255)                                                                                                                                 gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  >  3:281892/1‑149 (MQ=255)                                                                                                                                 gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  2:154404/149‑1 (MQ=255)                                                                                                                                 gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTACTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  8:145222/149‑1 (MQ=255)                                                                                                                                  ttGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTgg                 <  1:295197/149‑1 (MQ=255)                                                                                                                                         aaataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaatg           >  8:283718/1‑148 (MQ=255)                                                                                                                                         aaataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCAGGtgaatg           >  1:134280/1‑148 (MQ=255)                                                                                                                                           ataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa        <  2:162024/149‑1 (MQ=255)                                                                                                                                                cAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGtg   <  7:300593/149‑1 (MQ=255)                                                                                                                                                 aGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGtga  <  8:340145/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGA  >  NC_000913/1398407‑1398699

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