Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	50.0%	58.7 / 31.6	60	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/23);  new base G (30/0);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 1.74e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213195‑213489                                                                                                                                                   |                                                                                                                                                      gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAggg                                                                                                                                                   >  1:992461/1‑149 (MQ=255)     cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTa                                                                                                                                                 <  1:80523/149‑1 (MQ=255)       cccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGc                                                                                                                                             <  2:966499/151‑1 (MQ=255)           gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGAAAGCTTa                                                                                                                                          >  2:1328012/1‑150 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCg                                                                                                                                     >  2:881840/1‑151 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGc                                                                                                                                      <  2:1119565/150‑1 (MQ=255)                 gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCggg                                                                                                                                   <  2:165710/151‑1 (MQ=255)                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCgggggg                                                                                                                                >  1:1257553/1‑151 (MQ=255)                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                >  2:1424145/1‑151 (MQ=255)                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                                 <  1:365220/150‑1 (MQ=255)                     gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTTAAGCTTAATGCggggggg                                                                                                                               >  1:935243/1‑150 (MQ=255)                     gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGTACAGCTTAATGCGGGAggg                                                                                                                               >  2:614367/1‑150 (MQ=255)                      ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTTATGCgggggggg                                                                                                                              >  1:651135/1‑149 (MQ=255)                      ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCg                                                                                                                              <  1:1303340/151‑1 (MQ=255)                             ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTc                                                                                                                        <  1:631875/150‑1 (MQ=255)                               gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAAAGCGGGGGGGGATATTTGcc                                                                                                                     >  1:1017296/1‑151 (MQ=255)                                     cGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTTGCCCTcccc                                                                                                               >  1:510810/1‑149 (MQ=255)                                      gACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACGGTTTAATGCGGGGGGGGGTATTTCCCCCCcgcg                                                                                                              >  1:120120/1‑151 (MQ=255)                                         gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGc                                                                                                           <  1:606886/151‑1 (MQ=255)                                           tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTTATGCGGGGGGCGATATTCGCCCCCCCCGTGcgg                                                                                                         >  1:192860/1‑149 (MQ=255)                                            aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGa                                                                                                        <  1:43632/151‑1 (MQ=255)                                              cAGTCGCAACTGTGGTGGATTAAAGCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGAc                                                                                                       <  1:868622/150‑1 (MQ=255)                                               aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGGGGCGATATTCGCCCCCCCCGGGCAGACg                                                                                                      >  2:905050/1‑150 (MQ=255)                                                gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                    <  1:432458/151‑1 (MQ=255)                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCATCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                            <  2:211582/150‑1 (MQ=255)                                                              ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCCTGCAGACGAAGCGGGCAgcgtgcg                                                                                      >  2:540829/1‑151 (MQ=255)                                                                 ttAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGCGGTTGGGAGTGGACAACTTTATTCGGGGGGCGAAATTCGCCCTCCGCGTGGAGACCAAACGGTCAgcccgcggtt                                                                                   >  1:1403590/1‑143 (MQ=255)                                                                   aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATATTCGCCCCCCGCGTGCAGACGAAGCGGTCAGCGGGCGgttt                                                                                  >  2:630769/1‑149 (MQ=255)                                                                    aaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCaa                                                                                <  2:1434482/151‑1 (MQ=255)                                                                      tCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGGGCAGACGAAGCGGTCAGCGGGGGGTTCaa                                                                                >  2:462794/1‑149 (MQ=255)                                                                      tCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGGTATTTTCCCTCCCCGTGCAGACGAAGCGGTCAGGGGGCGGTTCaaa                                                                               >  1:442458/1‑150 (MQ=255)                                                                         gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGAGGACAGCATAATGGGGGGGGGGATATTTTCCCCCCGCCGGCAGACGAAGCGGGGAGCGTGCGTGTCAAAacgc                                                                           >  1:1209362/1‑151 (MQ=255)                                                                           cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGATATTTGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGGTTCAAAGCGcc                                                                          >  2:1173752/1‑150 (MQ=255)                                                                           cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCCCCGCGTGCAGACGAAACGGGCAGCGGGCGGTTTAAAGCGccc                                                                         >  2:5322/1‑150 (MQ=255)                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAATGCGGGAGGCGGAATTTCCCCTCCGCCGGGCGACCAAGCGGGCCGCGGGGGTTTTaaaaagcccgggt                                                                     >  1:205823/1‑141 (MQ=255)                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGAt                                                                    <  2:1258484/151‑1 (MQ=255)                                                                                  cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATtgc                                                                 <  2:989264/151‑1 (MQ=255)                                                                                   aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGTGTACAGCTTTATTCGGGGGGCGATATTTGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTCAAAGCGCCA‑GGGTtgct                                                                >  1:632191/1‑151 (MQ=255)                                                                                     aGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATtgctgc                                                              <  2:388042/151‑1 (MQ=255)                                                                                      gCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGGTATTTGCCCCCCCCGTGCAGACGAAACGGGCAGCGGGCGGTTCAAAACGCCC‑GGAGtgctgc                                                              >  1:627321/1‑150 (MQ=255)                                                                                        gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGTGCAGGCTAATGCCGGAGGGGAAATTCTGCCTCCCCCGGCGGACGGCGCGGGGAGCGGGGGGTTCAAAAGGCCA‑gggattctgcgt                                                            >  1:1449183/1‑140 (MQ=255)                                                                                          aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGAGGAAAGCTTAATGCGGGGGGAGAGATTCGGCCTCCGCGGGGAGAGGAAACGGGCAGGGGGGGGTTCAAAACGGCC‑GGATTGGTGCataat                                                         >  2:1378219/1‑149 (MQ=255)                                                                                              aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCCTAATGCGGGAGGGGATATTTGCCCCCCCCGTGCCGACGAAGCGGTCAGCGGGCGGTTCAAAACGCCC‑GGATTGCTGCATATTTTCg                                                     >  1:1363721/1‑151 (MQ=255)                                                                                              aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAAAGCGGGGGGCGGTATTTCCCCCCCCCGGGCAGACGAAACGGGCAGCGGGGGGTTTAAAAAGCCAGGGATTG‑TGCATATTTTCg                                                     >  1:661331/1‑151 (MQ=255)                                                                                               ttGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATATTCTCCCCCCCCGTGCAGACGAAACGGGCAGCGGGCGTTTTAAAGCGGCA‑GGGTTGCTGCCTATTTTCgg                                                    >  1:334978/1‑151 (MQ=255)                                                                                                   gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGATATTTTCCCCCCGCCTGCAGACGAAGCGGGCAGCGTGCGGTTTAAAACGCCC‑GGGTTGCTGCCTAATGTTGGGCgg                                                >  1:1003383/1‑150 (MQ=255)                                                                                                   gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTCAAAGCGCCA‑GGGTTGCTGCATATTGTCGGGCgg                                                >  2:58285/1‑150 (MQ=255)                                                                                                      gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAc                                             <  2:1111118/151‑1 (MQ=255)                                                                                                          cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAcggcg                                         <  2:209869/151‑1 (MQ=255)                                                                                                          cAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGACGGGGAGGGTACAGCTTAATACGGGGGGCGAAATTCGCCCTCCGCGGTGAGACGAAGCGGGCAGCGTGCGTTTAAAAGCGCCA‑GGGTTGCTGCATATTGGCGGGCGTAAcggcg                                         >  1:1377725/1‑151 (MQ=255)                                                                                                                        aCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAg                           <  2:21049/151‑1 (MQ=255)                                                                                                                          cGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGATTAATCTTGGAGGCGATATTCGCCCACCGCGTGAAGACGATGCGGTCAGCGTGCGTTTCAACTCGACA‑CAATTGCTCCATATTGACGGATGTAACGGCGGCCCTAAGCTAAAGaa                         <  1:1160113/151‑1 (MQ=255)                                                                                                                              ggAATTACCGGCGGGGCTGGGGAGTGTAAAGCCTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGAAGAAGCGGGCAGCGTGCGGTTCAAAGCGCCA‑GGGTTGCTGCATATTGTCGGGGGTAACGGCGGGCGTAAGCTAAAGaaaaac                     >  1:1045106/1‑149 (MQ=255)                                                                                                                                     ccGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag              <  2:1072194/151‑1 (MQ=255)                                                                                                                                     ccGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag              <  1:1319901/151‑1 (MQ=255)                                                                                                                                       ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag              >  1:524039/1‑149 (MQ=255)                                                                                                                                        gcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc           <  1:1424150/151‑1 (MQ=255)                                                                                                                                         cggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt          >  2:677557/1‑151 (MQ=255)                                                                                                                                                gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg   >  1:290274/1‑151 (MQ=255)                                                                                                                                                 ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgc  >  2:63748/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213195‑213489

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