Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,604,237	0	C	T	40.0%	75.9 / 32.2	61	S64L (TCG→TTG)	lsrD	autoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (6/30);  new base T (24/0);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 3.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604091‑1604386                                                                                                                                                   |                                                                                                                                                      cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                                   >  1:1035703/1‑151 (MQ=255) cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGtt                                                                                                                                                    <  2:1232854/150‑1 (MQ=255) cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGt                                                                                                                                                     <  2:963292/149‑1 (MQ=255)    tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTgg                                                                                                                                                 <  2:129817/150‑1 (MQ=255)    tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTgg                                                                                                                                                 <  1:38135/150‑1 (MQ=255)     aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTttgttttgtt                                                                                                                                               >  1:780390/1‑151 (MQ=255)     aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCACTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGt                                                                                                                                                <  2:566093/150‑1 (MQ=255)         agaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCGATATTTGGTTTTGtt                                                                                                                                               >  1:617612/1‑147 (MQ=255)         agaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCgtc                                                                                                                                           <  1:484274/151‑3 (MQ=255)           aTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCa                                                                                                                                         <  2:1094015/151‑1 (MQ=255)              gTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATc                                                                                                                                       <  2:1215622/150‑1 (MQ=255)                 gCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGcc                                                                                                                                   <  1:1008131/151‑1 (MQ=255)                  cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTTTTTTGTTTGTCCCCCGGcc                                                                                                                                   >  1:1121416/1‑150 (MQ=255)                  cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGAAGGTATTTTGTTTTGTTTGAACCTCgg                                                                                                                                     >  2:1256080/1‑148 (MQ=255)                   aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCctc                                                                                                                                 <  1:473253/151‑1 (MQ=255)                         tGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTCGACCATCGGCCTCTGCGcc                                                                                                                           >  1:578326/1‑151 (MQ=255)                          gCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGCCGATATTTTGTTTTGTTCGACCCTCGGCCTCTGCGcc                                                                                                                           >  1:479999/1‑150 (MQ=255)                           cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAACGAAATTTTCTTTGGTTCGACCATCGGCCTCTg                                                                                                                               >  1:786575/1‑145 (MQ=255)                            aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGTGTCGATATTTCGTTTGGGTTGACCCTCGGCCTCCGCGCcctt                                                                                                                        >  1:1077718/1‑151 (MQ=255)                            aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCCGCCTCTGCGCcctt                                                                                                                        >  1:1183101/1‑151 (MQ=255)                            aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTTGACcctccgcccctccgcccct                                                                                                                        >  1:886419/1‑148 (MQ=255)                                  cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcattg                                                                                                                  <  1:866514/151‑1 (MQ=255)                                       gAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTGGTTTGCCCATCCGCCCCCGCGCCCTTGCATTggggg                                                                                                              >  2:1152314/1‑148 (MQ=255)                                       gAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGt                                                                                                             <  1:177021/151‑1 (MQ=255)                                           gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGTTTGATATTTTTGCTTGTTTGACCCTTCGCCTCTTCGCCCATGCCATGGGCGt                                                                                                             >  2:603921/1‑147 (MQ=255)                                            ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGt                                                                                                        >  1:1067697/1‑151 (MQ=255)                                            ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTTGACCATCGGCCTCCGCGCCATTTCCTTTGGCg                                                                                                              >  1:40010/1‑145 (MQ=255)                                             tAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCCGCCTCTGCGGCCTCGCATTTGGCGTGCTGtt                                                                                                       >  2:881377/1‑151 (MQ=255)                                             tAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGGTTTACCCTCGGCCGCTGCGCCCTTTCATTGGGGGGGCTGtt                                                                                                       >  1:725775/1‑151 (MQ=255)                                                    cAATATGTTGCTGTTCAGCGCCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCaaa                                                                                                  <  1:1188881/149‑1 (MQ=255)                                                       tatGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGCGTGCAAATTTTTGTTGGTTTTGCCATCGGCCCCTTCCCCCTTTCATTTGGGGGACTGttttcaatggg                                                                                              >  1:747485/1‑142 (MQ=255)                                                          gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTTCATTTGGCGTACTGTTTCAAAGTGGggtg                                                                                          >  1:30711/1‑151 (MQ=255)                                                          gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGgtgt                                                                                           <  1:1162474/150‑1 (MQ=255)                                                            tGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGACCATATTTTTTTTTGTTTGACCAACGGCCTCTGCGCCCTTTCATTTGGGGGTCTTTTTTTAAGTGGggtg                                                                                          >  2:743196/1‑149 (MQ=255)                                                               tGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATATTTTGTTTGGTTCGGCCATCGGCCCCTGCGCCATTTCATTTGGGGTACTGTTTTTAAgtgggg                                                                                            >  2:546119/1‑142 (MQ=255)                                                                       aCCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATGGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTgg                                                                             <  2:891354/151‑1 (MQ=255)                                                                            tGACTTTATCTGCATTGGCATTGTCGCCC‑ACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGAATTTTGTTTTGTTTCACCCTCGGCCCCTGCGCCCCTGCCTTTGGGGGCATGTTTTAAAGTGGGGTGGCGGGGCCGCCgggg                                                                           >  2:1248561/1‑145 (MQ=255)                                                                              aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTGGTTTGACCATCGGCCTCTTGGCCATTGCATTGGGGGTACTGTTTTTAAAGGGGGGGGCGATGCCCCTGGCGATACt                                                                      >  1:1072405/1‑151 (MQ=255)                                                                               cTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGAAAATTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTTCATTTGGGGGACTGTTTCAAAGAGGTGTGGCCGTGCCGCCGGCGAttc                                                                       >  2:525235/1‑147 (MQ=255)                                                                                tttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGGCCCTCCGCCCCTGCCCCCTTGCATTTGGCGGAGTGGTTTAAAATGGGGTGGCGATTCCCCCGGCGGGACTcc                                                                    >  2:455315/1‑151 (MQ=255)                                                                                       gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGAcc                                                               <  1:866734/149‑1 (MQ=255)                                                                                       gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGAcc                                                               <  2:700751/149‑1 (MQ=255)                                                                                          ttggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGTTCCGTTTTTCGTTTTTTTTGACCCACCGCCTCCGGCCCCTTTCAATTGGGGTACTTTTTCATAATGGTGTGCCGCCGCCGCTGGCGGTACTCCTCAcc                                                               >  1:1051010/1‑146 (MQ=255)                                                                                           tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACt                                                          <  1:303289/150‑1 (MQ=255)                                                                                                ttgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATTGGGCTCTGCGCCCTTTCATTTGGGGTACTGTTTCAAAGTGGTGGGGCGCTGGCGCCGGCGGTACTCCTTACCtt                                                             >  1:42166/1‑142 (MQ=255)                                                                                                  gTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCAGCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGgc                                                   <  2:293257/150‑1 (MQ=255)                                                                                                      cccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTCGTTTGGTTCGCCCATCGGCCCCTGCGCCACTGCCTTGGGCGGACTGTTTTTAAATGGGGTGCCGATGCCGCTGGCGATACTCCTTTCCTTACTGCTCGGCGCAt                                               >  1:764381/1‑150 (MQ=255)                                                                                                       ccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtg                                             <  2:30709/151‑1 (MQ=255)                                                                                                          aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                           <  1:388842/150‑1 (MQ=255)                                                                                                               tAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCATTTGGGGGACTGTTTCAAAGTGGGGTGCCGATGCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGCATTGTGCggg                                       >  1:601436/1‑149 (MQ=255)                                                                                                                      ggTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAAc                              <  2:1103745/151‑1 (MQ=255)                                                                                                                             gTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGa                        <  2:1087697/150‑1 (MQ=255)                                                                                                                               cAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATTCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatt                      <  2:984019/150‑1 (MQ=255)                                                                                                                               cAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatt                      <  2:812151/150‑1 (MQ=255)                                                                                                                                     cGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTAt                <  2:786546/150‑1 (MQ=255)                                                                                                                                     cGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATc               >  1:799973/1‑151 (MQ=255)                                                                                                                                      gggATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                  <  2:646943/147‑1 (MQ=255)                                                                                                                                         aTCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtata           <  1:387505/151‑1 (MQ=255)                                                                                                                                         aTCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtat            >  1:394549/1‑150 (MQ=255)                                                                                                                                            gATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa        <  1:95782/151‑1 (MQ=255)                                                                                                                                            gATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa        <  1:148092/151‑1 (MQ=255)                                                                                                                                                  tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt   >  2:1009677/1‑150 (MQ=255)                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa  >  2:978902/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604091‑1604386

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