Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,852	G→A	S51F (TCT→TTT)	yaeF ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,852	0	G	A	98.6%	233.8 / ‑4.4	69	S51F (TCT→TTT)	yaeF	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (32/36);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTT  >  NC_000913/216708‑216996                                                                                                                                                 |                                                                                                                                                 cccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCt                                                                                                                                              >  8:236532/1‑149 (MQ=255)  ccAGAAAGATTGCAACGTGAGTCCCAGAGGAAGTGCTGATGACGCGGATTCCTAATGAGGTTACGCCAAGGCTTGAGGAGTTCAGCAAATCGCCGGGTGTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCtt                                                                                                                                             <  1:54447/149‑1 (MQ=255)  ccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCtt                                                                                                                                             >  1:160568/1‑149 (MQ=255)  ccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCtt                                                                                                                                             >  7:96625/1‑149 (MQ=255)    aGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCttt                                                                                                                                            >  1:226357/1‑148 (MQ=255)    aGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGATCGGTGAAAAAGCttt                                                                                                                                            >  8:11729/1‑148 (MQ=255)         gATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCt                                                                                                                                      <  8:273602/149‑1 (MQ=255)          aTTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTg                                                                                                                                     <  8:336048/149‑1 (MQ=255)             gCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGaa                                                                                                                                  >  8:32196/1‑149 (MQ=255)              cAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCGGGaa                                                                                                                                  <  5:183886/148‑1 (MQ=255)                 cGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAAttt                                                                                                                               >  7:37124/1‑148 (MQ=255)                  gTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGa                                                                                                                             <  7:236532/149‑1 (MQ=255)                   tGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGAc                                                                                                                            <  8:144663/149‑1 (MQ=255)                         acaGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCa                                                                                                                      >  3:61420/1‑149 (MQ=255)                           agagGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTCCCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTATATTGATTGTTCGGTGAAAATGATTTGATGCTGGAATTTGACAGCCCATg                                                                                                                    >  7:288669/1‑149 (MQ=255)                                aaGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGttttt                                                                                                               <  2:4512/149‑1 (MQ=255)                                aaGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGttttt                                                                                                               >  8:37273/1‑149 (MQ=255)                                  gTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGc                                                                                                             <  2:226357/149‑1 (MQ=255)                                     cTGAAGACGCGGATTCCAAATGAGGTTACCCCTAGGCTTGATGAGAACAGCAAATCGCCGGGGTTGAGATCAGGTGACGTTATTTCTTTTATTGATTGTTCGGTTAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTc                                                                                                          <  8:288669/149‑1 (MQ=255)                                     cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTc                                                                                                          >  1:173918/1‑149 (MQ=255)                                     cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTc                                                                                                          >  3:343791/1‑149 (MQ=255)                                      tGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCg                                                                                                         >  2:257925/1‑149 (MQ=255)                                        aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTATAATTTATTGTTCGGTGAAAGAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCg                                                                                                         <  8:118307/147‑1 (MQ=255)                                          gACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCAt                                                                                                     <  5:300233/149‑1 (MQ=255)                                             gcgGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTa                                                                                                  >  6:347606/1‑149 (MQ=255)                                                     cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCtgt                                                                                           <  8:96625/148‑1 (MQ=255)                                                     cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCtgt                                                                                           <  8:37124/148‑1 (MQ=255)                                                     cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCAGTTATTTCTTTTATTCATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCAATGTTTTTGCATCGGCATCTACTGCtgtg                                                                                          >  8:256975/1‑149 (MQ=255)                                                      aaaTGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTgg                                                                                         <  4:115/149‑1 (MQ=255)                                                       aaTGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGc                                                                                        <  8:23082/149‑1 (MQ=255)                                                       aaTGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGc                                                                                        <  4:61420/149‑1 (MQ=255)                                                       aaTGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTAGGTGAAAAAGCTTTGCTGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGc                                                                                        >  5:265785/1‑149 (MQ=255)                                                       aaTGAGGTTACCCCAAGGCTGTAGGAGAACTCCATCTCGCCGTGCTTGAGAACAGGTGCCGTTATTTCTTTTATTGATTGCTCGGAGAAAAAGCTTTGATGCTGTACTTTGACAGCCCATGTTTTTGCCTCTGCATCTACTGCTGTGGc                                                                                        <  4:6218/149‑1 (MQ=255)                                                                 ccccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTct                                                                              <  2:267388/149‑1 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATc                                                                                                    >  1:73214/1‑126 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATc                                                                                                    <  2:73214/126‑1 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGAt                                                                                    <  3:334689/142‑1 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGAt                                                                                    >  4:334689/1‑142 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTctg                                                                             >  7:60641/1‑149 (MQ=255)                                                                  cccAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTctg                                                                             >  3:258277/1‑149 (MQ=255)                                                                     aaGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTctggct                                                                          <  4:258277/149‑1 (MQ=255)                                                                        gCTTGAGGAGAACAGCAACTCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTAATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCtggctggc                                                                       <  6:290363/149‑1 (MQ=255)                                                                        gCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCtggctggc                                                                       <  1:257925/149‑1 (MQ=255)                                                                          ttGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTggctggct                                                                      >  2:133443/1‑148 (MQ=255)                                                                            gaggagAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGa                                                                    <  5:347606/148‑1 (MQ=255)                                                                            gaggagAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGCTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGa                                                                    <  7:11729/148‑1 (MQ=255)                                                                               gagaACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATc                                                                <  7:32196/149‑1 (MQ=255)                                                                                  aaCAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGa                                                                                                                             >  7:31181/1‑85 (MQ=255)                                                                                  aaCAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGa                                                                                                                             <  8:31181/85‑1 (MQ=255)                                                                                           tCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGc                                                    <  7:37273/149‑1 (MQ=255)                                                                                              ccGGGTTTGAGATCAGGTGCCGTTAGTTCTTATATTGATAGTTCGGTGAAAAAGCTTTGATGCTGGAATTAGACAGCCAATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGTTACAATCCg                                                  >  6:363809/1‑148 (MQ=255)                                                                                                gggTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGaaa                                               <  7:53548/149‑1 (MQ=255)                                                                                                          tCAGGTGCCGTTAATTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAaaagaaa                                     <  6:301179/149‑1 (MQ=255)                                                                                                              gTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGata                                                                  >  4:94852/1‑116 (MQ=255)                                                                                                              gTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGata                                                                  <  3:94852/116‑1 (MQ=255)                                                                                                              gTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTg                                 <  8:60641/149‑1 (MQ=255)                                                                                                              gTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTg                                 <  1:133443/149‑1 (MQ=255)                                                                                                                gCCGTTATTTCTTTTATTTATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTggc                               >  6:301752/1‑148 (MQ=255)                                                                                                                 ccGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTggg                               >  8:35828/1‑148 (MQ=255)                                                                                                                  cGTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTggg                               >  4:178146/1‑147 (MQ=255)                                                                                                                   gTTATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTggg                               >  3:165554/1‑146 (MQ=255)                                                                                                                    ttATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGAt                                                                                    <  1:8854/92‑1 (MQ=255)                                                                                                                    ttATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGAt                                                                                    >  2:8854/1‑92 (MQ=255)                                                                                                                    ttATTTCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGgag                            >  6:274029/1‑148 (MQ=255)                                                                                                                         tCTTTTATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACa                       <  3:178146/148‑1 (MQ=255)                                                                                                                           ttttATTGATTGTTCGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAggc                    >  3:358269/1‑149 (MQ=255)                                                                                                                                         cGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCttt      >  6:211827/1‑149 (MQ=255)                                                                                                                                         cGGTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCtt       <  2:160568/148‑1 (MQ=255)                                                                                                                                          ggTGAAAAAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTg     <  4:343791/149‑1 (MQ=255)                                                                                                                                              aaaaAGATTTGATGCTGGAATTTGACAGCCCATGTATTTGACTCTGCATCTACTGCTGAGGCAGATGGGTCTGGCACGCTCATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAAAAGGCGGCAGTACGCCTATGCtt  >  5:263656/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGCTTTGATGCTGGAATTTGACAGCCCATGTTTTTGCCTCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTT  >  NC_000913/216708‑216996

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