Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,056,115	0	A	C	27.1%	142.1 / 11.7	59	G75G (GGT→GGG)	glnL	sensory histidine kinase NtrB
Reads supporting (aligned to +/- strand):  ref base A (39/4);  new base C (0/16);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 4.42e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGG  >  NC_000913/4055979‑4056253                                                                                                                                         |                                                                                                                                            gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTTAAACC‑t                                                                                                                                         >  2:1147296/1‑140 (MQ=255)  ttCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑tt                                                                                                                                        >  1:1339547/1‑140 (MQ=255)   tCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTg                                                                                                                                       >  1:101444/1‑140 (MQ=255)      tGACTTAGGCGGCGCTGGTTATCCATGGGAGCCATCCCCAGCAGGATCATGCCGCCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGccc                                                                                                                                    <  2:1538507/140‑1 (MQ=255)      tGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGccc                                                                                                                                    >  1:877293/1‑140 (MQ=255)      tGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGccc                                                                                                                                    >  2:691456/1‑140 (MQ=255)      ggACTTAGGGGGCGCTGGTTACCCATGGGGGCCATCCCCAGCGGGATCAGGCCGCCCGGCATCCGCTGGGCCGTCACAGAAAGGATATGCGGGCGCCGGTCGATGCCCAGCGTCACTTCGTTATGGGTAAACCC‑TTGccc                                                                                                                                    <  2:1534261/139‑1 (MQ=255)                 gcgcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGAc                                                                                                                         >  1:586274/1‑140 (MQ=255)                 gcgcTGGTTATCCAGCGGGCCCATCCCCAGCAGGTTCTTGCCCCCCGGCAACCGCGGGGCCTTCACAGAAAGGATAGGCGCGCCCCCGTCGTTGACCAGCGTCACTTCGTTATCGTTAAACCC‑TTGCCCCGCCTCCAGAc                                                                                                                         <  2:635563/140‑1 (MQ=255)                  cgcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACt                                                                                                                        >  1:3606/1‑140 (MQ=255)                         tATCCAGCGGGCCCACCGCCAGCAGGGTCGTGCCGCCGGGCATGCGTGGGGCGTTCACAGAAAGGATATGCGGGCGCCCGTGGATGCCGAGCGTCACTTCGTTACCGGAAAATCT‑TTCCCCCGCCTCCAGACTTTCTTGc                                                                                                                 <  2:1370712/140‑1 (MQ=255)                            ccATCGGAGCCATCCCCAGCAGGTTCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAACCC‑TTCCCCCGCCTCCAGACTTTCTTGCATc                                                                                                              <  2:556874/140‑1 (MQ=255)                               tCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGc                                                                                                           >  1:508761/1‑140 (MQ=255)                                 ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                         >  1:1543468/1‑140 (MQ=255)                                  gcgCCACCCCCGGCGGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATGGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCg                                                                                                        <  2:1343442/138‑1 (MQ=255)                                  gAGCCACCTCCAGCAGGACCATGCCGCCCGCCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCCCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCg                                                                                                        <  2:448072/140‑1 (MQ=255)                                           ccAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAAt                                                                                               <  1:1417664/140‑1 (MQ=255)                                               cagGATCATGCCGTCCGGCATACGCTGGGTCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgaga                                                                                           <  1:523773/140‑1 (MQ=255)                                                  gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                                                                                        >  1:147227/1‑140 (MQ=255)                                                  gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                                                                                        >  2:308183/1‑140 (MQ=255)                                                  gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                                                                                        >  1:622232/1‑140 (MQ=255)                                                    tCATACCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                      >  2:908493/1‑140 (MQ=255)                                                    tCAGGCCGCCCGGCATCCGCGGGGCCGCCACAGAAAGGATATGCGGGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGAAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                      <  2:652272/140‑1 (MQ=255)                                                         ccggccggCAAACGCTGGGCCTCCACAGAAAGAATATGCGAGCGCCCGGGGAGGACCAGCTTCTCTTCGTAATCGTAAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCaa                                                                                 <  2:1339547/140‑1 (MQ=255)                                                         ccGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCaa                                                                                 <  2:309813/140‑1 (MQ=255)                                                           gcccGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCTTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                               <  2:594717/138‑1 (MQ=255)                                                           gTCCGGCATCCGTTGGGCCGTCACAAAAAGGTTCTGGGGGCGCCCGTGGAGGGCCAGCGTCCCTTCGTTA‑CGGGAAAACCTTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                               <  2:781021/140‑1 (MQ=255)                                                           gTCCGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGGGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                               <  2:614588/140‑1 (MQ=255)                                                                gCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTcc                                                                          >  2:851254/1‑140 (MQ=255)                                                                     cggggggCCGTCAAAGAAAGGATATGCGGGCGCCCGTCGATGCCCGGCGTCTCTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa                                                                     <  2:508761/136‑1 (MQ=255)                                                                     cGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa                                                                     <  1:207222/140‑1 (MQ=255)                                                                      gCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAc                                                                    <  1:1127518/140‑1 (MQ=255)                                                                      gCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAc                                                                    >  2:909993/1‑140 (MQ=255)                                                                       cTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACg                                                                   <  2:527815/140‑1 (MQ=255)                                                                        ggcgCCGTCACAGAAGGGATAGGCGGCCGCGCTGCGTGGCCCAGCGTTTCTTCTTTTTGGTTAAACCT‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACgg                                                                  <  2:633060/137‑1 (MQ=255)                                                                         gggCCGTCACAGAAAGGATAGGCGGGCCCCCGTCGATGCCCAGCGTCATTTCGTTATCGGAAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgt                                                                 <  2:3606/140‑1 (MQ=255)                                                                          ggCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtg                                                                >  2:989120/1‑140 (MQ=255)                                                                          ggCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtg                                                                >  1:37028/1‑140 (MQ=255)                                                                           gCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtgt                                                               >  1:1248027/1‑140 (MQ=255)                                                                               tCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                           >  1:562627/1‑140 (MQ=255)                                                                                 acaGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCaaa                                                         <  1:1213249/140‑1 (MQ=255)                                                                                      aaGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAtt                                                    >  2:360071/1‑140 (MQ=255)                                                                                         gATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGc                                                 >  1:1428431/1‑140 (MQ=255)                                                                                          atatGCGGGCGCCCGTCGATGCCCGGCGTCCTTTCGTTTGCGGGAACCCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCg                                                <  2:622232/140‑1 (MQ=255)                                                                                          atatGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCg                                                <  2:147227/140‑1 (MQ=255)                                                                                           tatGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCgg                                               >  1:1020399/1‑140 (MQ=255)                                                                                            atGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCggg                                              >  1:876266/1‑140 (MQ=255)                                                                                             tGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGa                                             >  2:395665/1‑140 (MQ=255)                                                                                             tGCGAGCGCCCGTCGATGACCAGCGTCACCTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGa                                             >  2:826554/1‑140 (MQ=255)                                                                                                 cgcgcCCGTCGAGGACCAGCGTCATTTCTTTTCCGTTAAACCT‑TTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                         <  2:586274/139‑1 (MQ=255)                                                                                                 aGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                         >  1:437025/1‑140 (MQ=255)                                                                                                 aGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                         >  1:734818/1‑140 (MQ=255)                                                                                                 aGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                         >  1:871964/1‑140 (MQ=255)                                                                                                   cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTg                                       >  1:668434/1‑140 (MQ=255)                                                                                                   cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTg                                       >  1:1162160/1‑140 (MQ=255)                                                                                                   cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTg                                       >  2:784914/1‑140 (MQ=255)                                                                                                   cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTg                                       >  1:883843/1‑140 (MQ=255)                                                                                                     cccGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTggg                                     <  1:851254/140‑1 (MQ=255)                                                                                                     cccGTCGATGACCAGCCTCACTTCTTTTTCGTAAAACCC‑TTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTggg                                     <  2:1543468/140‑1 (MQ=255)                                                                                                       cGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCg                                   >  1:1576202/1‑140 (MQ=255)                                                                                                                   gCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgc                       >  2:794017/1‑140 (MQ=255)                                                                                                                    cGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgcg                      >  2:934700/1‑140 (MQ=255)                                                                                                                    cGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgcg                      >  2:1349989/1‑140 (MQ=255)                                                                                                                        ttttCGTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGc                  <  2:1428431/138‑1 (MQ=255)                                                                                                                         cTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCa                 >  1:17307/1‑140 (MQ=255)                                                                                                                         cTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCa                 >  1:936022/1‑140 (MQ=255)                                                                                                                             tttATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt             <  2:1223723/139‑1 (MQ=255)                                                                                                                             gTTATCGGTAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt             <  1:165699/140‑1 (MQ=255)                                                                                                                             gTTACCGGAAAACCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt             <  1:1147296/140‑1 (MQ=255)                                                                                                                              ttATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGtt            >  2:922646/1‑140 (MQ=255)                                                                                                                                aTCGGAAAACCC‑TGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTgg          <  2:29120/140‑1 (MQ=255)                                                                                                                                    gTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTa      >  1:1063492/1‑140 (MQ=255)                                                                                                                                      aaaaCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAAt    >  2:1156612/1‑140 (MQ=255)                                                                                                                                       aaaCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATg   <  1:989120/140‑1 (MQ=255)                                                                                                                                        aaCC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATgg  >  2:485024/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                            GTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACC‑TTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGG  >  NC_000913/4055979‑4056253

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