Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	52.2%	62.0 / 100.1	67	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (27/5);  new base C (0/35);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 4.91e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-10

GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTG  >  NC_000913/907454‑907746                                                                                                                                                       |                                                                                                                                                 gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGAAATCACATTGAGAAACCACACCGAAATCCAGCGGGCGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTAcc                                                                                                                                                >  2:375082/1‑151 (MQ=255)             catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACGTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCa                                                                                                                                     <  2:901689/150‑1 (MQ=255)                 cGCCAGCAGATCGTTTTTCAGCCCCTAAACCACCCCCCGGGGGGGGGGGGGCGGTAGGAAAATAATCACATTGAGAACCCACCCCAAAACCCAGGGGGGGGGGGGGGGGATGGGGATC‑AACCA‑GCGGCGTTTCCCCCCCACCGCCATaaaa                                                                                                                                <  1:880618/151‑1 (MQ=255)                  gCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAAt                                                                                                                               >  2:684610/1‑151 (MQ=255)                   ccAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCGGGGGGGGGTGGGCGGTAGGGAAGTAATCACATTGAAAACCCACCCCAAAACCCAGGGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTCCCCCTCACCGCCATaa                                                                                                                                  <  1:411132/147‑1 (MQ=255)                         gATGGTTTTTCAGCCCCTAAATCACCGCCCGGGGGGGGGGGGGGGGTACGAAAGTAAACACATAGAAAACCCACCCCAAAACCCAGGGGGGGGGGA‑GGGGAAGGGGAAA‑AACCAGGCGGCCCTCCCCGCTCACCGCCATAAAATGTTCaaa                                                                                                                        <  1:126902/151‑1 (MQ=255)                                tttCACCCCCTAAACCAACCCCCGGGGGGGGGGGGGCGGTACGGAAGTAATCACTTTGAAAACCCACACCGAAACCCAGGGGGGGGGGGGGGGGAGGGGGATC‑ACCGG‑GCGGCCCTTCCCGCCCACCGCCATAAAATGTTCAAACTGtt                                                                                                                   <  1:618658/149‑1 (MQ=255)                                   cacccccTAAACCACCGCCCGGGGGGGGGGGGGCGGTAGGAAAGAAACCAAATGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATA‑ACCCA‑GCGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                               <  2:1135578/147‑1 (MQ=255)                                   cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                               >  1:549933/1‑150 (MQ=255)                                      ccccTAAACCACCGCCCGGGGGGTGGGGGGCGGTAGGAAAGTAACCACATGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGACC‑AACCG‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                                             <  2:959095/149‑1 (MQ=255)                                            aaccaccGCCCGGGGGGTGGGGGGCGGTAGGGAAGAAATCACATGGAAAAACCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATC‑ACCCA‑GCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAg                                                                                                      <  2:961129/144‑1 (MQ=255)                                                 aCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                >  1:763585/1‑151 (MQ=255)                                                    cgcggggggtgggGGGGGGGACGGAAAAAACCACATGAAAAACCCCCACGAAAACCCAGCGGGGGGGGAGGGGGGTGGGGATA‑ACCCG‑GGGGCTTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                               <  1:727274/142‑1 (MQ=255)                                                    cccTGGGGGGGGGGGGGGGGCGGGAAGTAATCATATTGAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATC‑ACCCG‑GGGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTgg                                                                                             <  2:301567/151‑1 (MQ=255)                                                    cccTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTgg                                                                                             >  1:741679/1‑151 (MQ=255)                                                     ccggggggtgggGGGGGGTAGGAAAGAAACCACTTTGAAAACCCCCACAAAAACCCGGGGGGGGGGGGGGGGGAGGGGGAAA‑AACCG‑GGGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGt                                                                                            <  1:396359/145‑1 (MQ=255)                                                      cTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTc                                                                                           >  1:330913/1‑151 (MQ=255)                                                       tGGCGGTGGGGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                          >  2:1030433/1‑151 (MQ=255)                                                        ggggggggtggGGGTTAGGAAAGTAATCACATTGAGAACCCACCCCAAAATCCAGGGGGTGGGGAGGGGGATGGGGATC‑ACCGG‑GGGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                         <  2:8724/145‑1 (MQ=255)                                                           ggggggggggggTACGGAAGCAATCACATGAAAAACCCAACCCAAAATCCGGCGGGGGGGGGGGGGGAGGGGGAAC‑ACCCA‑GCGGCCCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                       <  2:137325/144‑1 (MQ=255)                                                                ggggCGGTACGGAAGTAATCACATTGAAAACCCACCCCGAAACCCAGGGGGTGGGGAGGGGGATGGGGATC‑ACCCA‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                 <  2:1109001/150‑1 (MQ=255)                                                                 gggggttaggaaataaaTCAAATGGAGAAACCACACCAAAACCCAGGGGGGGGGGGGGGGGGGGGGGATC‑AACGC‑GGGGCTCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                                <  1:1123410/136‑1 (MQ=255)                                                                 ggggggTAGGGAAGTAATCACATTGAAAACCCCCACCAAAACCGGGCGGGTGGGGGGGGGGAGGGGGATC‑ACCCA‑GCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                 <  1:30157/150‑1 (MQ=255)                                                                      gtaggaaaGTAACCACATGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAATC‑ACCCA‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                           <  2:179614/145‑1 (MQ=255)                                                                      gTAGGGAAATAATCAAATGAAAAACCCCCCCCAAAACCCGGGGGGGGGGGAGGGGGAGGGGGACC‑ACCCA‑GCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                           <  1:1032736/151‑1 (MQ=255)                                                                        agggAATTAATCACATGAAAAAACCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATC‑AACCA‑GCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                         <  2:936970/149‑1 (MQ=255)                                                                          ggAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                       >  1:867183/1‑151 (MQ=255)                                                                            aaGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                     >  1:892772/1‑151 (MQ=255)                                                                            aaGTAACCACATGGAAAAACCACACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATC‑ACCCG‑GGGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                     <  1:1092832/151‑1 (MQ=255)                                                                             aGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                    >  2:455166/1‑151 (MQ=255)                                                                                aaTCACATAGAGAAACCACCCCGAAACCCAGGGGGGGGGGGGGGGGAGGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                   <  2:20894/149‑1 (MQ=255)                                                                                 aTCACATTGAAAAACCACACCAAAACCCAGGGGGTGGGGGGGGGGAGGGGGAAC‑AACCA‑GCGGCGTTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                  <  2:218699/149‑1 (MQ=255)                                                                                   cacaTTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAACC‑ACCCG‑GCGGCGTTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGc                                                                <  2:442690/149‑1 (MQ=255)                                                                                    acaTAGAAAACCCACCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGAATA‑ACCCG‑GCGGTTCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                              <  1:80313/150‑1 (MQ=255)                                                                                     caTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGc                                                            >  1:213841/1‑151 (MQ=255)                                                                                        gaaaaacccACCCCAAAACCCGGCGGGGGGGGGGGGGGAGGGGAATC‑ACCCG‑GCGGCGCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                              <  1:485644/139‑1 (MQ=255)                                                                                          agaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                       >  1:1019841/1‑151 (MQ=255)                                                                                          aaaacccccacCGAAACCCGGGGGGGGGGGGGGGGGATGGGGATC‑AACCG‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGga                                                        <  2:205902/142‑1 (MQ=255)                                                                                          aaaacccACACCGAAACCGAGGGGGGGGGGGGGGGGAGGGGGATC‑AACCG‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                         <  2:529020/144‑1 (MQ=255)                                                                                           gaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                      >  1:667566/1‑151 (MQ=255)                                                                                             aaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                    >  1:13240/1‑151 (MQ=255)                                                                                             aaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                    >  2:697939/1‑151 (MQ=255)                                                                                               ccccacaaaaaccgggCGGGGGGGGAGGGGGATGGGGATC‑AACCG‑GGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                  <  2:435185/136‑1 (MQ=255)                                                                                               ccACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑ACCCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCCGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATc                                                   >  1:432337/1‑150 (MQ=255)                                                                                                cacaCCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                  >  1:90340/1‑150 (MQ=255)                                                                                                   acCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaac                                               >  1:711857/1‑150 (MQ=255)                                                                                                     caaaaTCGAGGGGGTGGGGGGGGGGAGGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                                             <  2:878938/148‑1 (MQ=255)                                                                                                      gAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                           >  1:236542/1‑151 (MQ=255)                                                                                                       aaacccgggggggggggggggggGTGGGGATA‑ACCCA‑GCGGCCCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                          <  1:684578/142‑1 (MQ=255)                                                                                                           ccAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                      >  1:1014159/1‑151 (MQ=255)                                                                                                             aGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                      >  2:30587/1‑149 (MQ=255)                                                                                                               cGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                      >  1:864072/1‑147 (MQ=255)                                                                                                                ggggggggggggggggggGAACC‑ACCCC‑GCGGCGTTCCCCGCCAGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                 <  1:519255/151‑1 (MQ=255)                                                                                                                 gggggggggggggAGGGGGATC‑ACCCA‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg                                <  2:612944/148‑1 (MQ=255)                                                                                                                  ggggggggggggggggggATC‑ACCCA‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGt                               <  2:835044/149‑1 (MQ=255)                                                                                                                    ggggggggggggggggAACAAACCG‑GGGGCTTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTg                              <  1:220961/151‑1 (MQ=255)                                                                                                                       gggggggagggggaTC‑AACCG‑GCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                          <  2:587745/149‑1 (MQ=255)                                                                                                                       gAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                          >  1:1127517/1‑151 (MQ=255)                                                                                                                            ggggggggacC‑ACCCA‑GCGGCGCTCCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                     <  1:247829/147‑1 (MQ=255)                                                                                                                             gatggggatC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                    >  1:24531/1‑151 (MQ=255)                                                                                                                                ggggctC‑AGCAA‑GCGCCACTAACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                 <  2:1126424/151‑1 (MQ=255)                                                                                                                                ggggatC‑AACCA‑GCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                    <  2:227191/148‑1 (MQ=255)                                                                                                                                    atC‑AACCA‑GCGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAg               <  2:892813/149‑1 (MQ=255)                                                                                                                                        aaCCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa           >  1:944080/1‑150 (MQ=255)                                                                                                                                        aaCCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa           >  1:706215/1‑150 (MQ=255)                                                                                                                                            a‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg      >  2:1142675/1‑151 (MQ=255)                                                                                                                                                  cgcTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg  <  1:803475/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATC‑AACCA‑GCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTG  >  NC_000913/907454‑907746

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