Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,457	0	T	G	36.2%	129.6 / 15.7	69	G232G (GGT→GGG)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base T (9/35);  new base G (25/0);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.26e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGA  >  NC_000913/3995310‑3995602                                                                                                                                                    |                                                                                                                                                   gCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTg                                                                                                                                                 <  1:1035996/151‑1 (MQ=255)   gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCTGTAGGGGTTTAGCAAGGgtttc                                                                                                                                                >  1:876943/1‑148 (MQ=255)   gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCt                                                                                                                                               <  1:1221690/151‑1 (MQ=255)       tCCGGAGCGCGCCAACATCAGTTATATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTGATGAGCGTGGGGCAGGAGAAACCCAGGGCTGCGGCTGAGGCGCGTGTGCGGCGGTTGCAATAGGGATTCAGCAAGGTTTGATg                                                                                                                                              <  2:1123687/148‑1 (MQ=255)         cGGAGCGCGCCAAGAGCGGTGTTATGCCAGGGCTTAAGCGCGGGCATATTCGGGGACGCTTTTAGGCGCGCGGGGCAGGAGAAGCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTGGGGATTCAGCAAGGTTTGCt                                                                                                                                               <  2:397948/145‑1 (MQ=255)               gcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                     >  1:971156/1‑149 (MQ=255)                cgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTa                                                                                                                                  <  1:254503/151‑1 (MQ=255)                cgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTa                                                                                                                                  <  1:652578/151‑1 (MQ=255)                 gcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTTAGCAAGGGTTGCTGGCCgaagaa                                                                                                                                     >  1:106141/1‑147 (MQ=255)                    aaTATCGGTTTTATGCAGGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGTGATTCAGCAAGGTATGCTGGCCGAAGAAGTAcgcg                                                                                                                              <  2:257691/151‑1 (MQ=255)                    aaTATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGTAagcg                                                                                                                              >  1:230786/1‑151 (MQ=255)                         cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGGACGCGTgg                                                                                                                           >  1:502470/1‑149 (MQ=255)                          tgTTTTATGCAAGTGTTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCAATAGGGGATCTGCACGGGTTGAAGGCCGAAGAACAAcgcggcgga                                                                                                                         >  2:869810/2‑145 (MQ=255)                           gTTTTATGCGAGGGGTTAAGCGCGAGCATATTCGTTTGCGCGTTTATGGGCGTGGGGCAGGAGAAACCCATACCTACGGCAGCAGGGCGTGTGAGGCGTTTGCAGTAGGGATTCAGCAAGTTTTGCTGCCCGAAGAAGTACGCGTGGAACt                                                                                                                       <  2:311168/151‑1 (MQ=255)                                  gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcg‑g                                                                                                                <  1:1183947/150‑1 (MQ=255)                                   cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcg‑gcg                                                                                                              >  1:1254436/1‑151 (MQ=255)                                   cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGGAGTAGGGAGTCAGCAAGGGTTGCCGGCCGCAGAAGTACGCGTGGGACACCCcc‑gc                                                                                                               >  2:214252/1‑150 (MQ=255)                                    aaGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCCGGCCGAAGAAGTACGCGGGGAACTCCCcg‑gcgg                                                                                                             >  1:959183/1‑151 (MQ=255)                                                    gCATATTCGTTTACGCGATTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGCGCGGCGGGTGCCGTAGGGCACTAGAAAGGGTTGCCGGCCGAACAAGCAAGCGTGGTAGTCCCCC‑CCGggcgtcctcttaacgcc                                                                                             >  2:922428/1‑140 (MQ=255)                                                     cATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGGAGAAAGAAGCGGGGGAATCCCCC‑GCGGGCGGCTTGGTATCGcc                                                                                             >  1:269007/1‑150 (MQ=255)                                                       tatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTTCTGGGCGAAGAAGTACGCGGGGGACTCCCCG‑GCGGGCGTCTTGATATTGCCCgg                                                                                          >  1:20133/1‑151 (MQ=255)                                                        atTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGTACCCGTGGAACTCCCCG‑GCGGGCGTCTTTATATCGcc                                                                                             >  1:1206786/1‑147 (MQ=255)                                                        atTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTgg                                                                                          <  2:211515/150‑1 (MQ=255)                                                        atTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCCGGCGGGAATCAACACAGGTTGCTTGCCCCAGGAAGACGCGTGGGACTCCCCC‑GCGGTCGTCTTTGGATAGCCTgg                                                                                          >  2:603981/1‑150 (MQ=255)                                                             tttACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGTTTGCCGAAGAAGTACGCGGGGAACTCCCCG‑GCGGGCGTCTTGATATTGCCTGGaaaag                                                                                     >  1:546674/1‑148 (MQ=255)                                                              ttACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTTAGCAAGGTTTGTTGGGCGGAGAAATAGGCGGGGAACTCCCCG‑GCGGTGGTTTTTATATTGCCtcggaacgg                                                                                    >  1:635243/1‑142 (MQ=255)                                                                 cgcgTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCgg                                                                                <  1:1016017/151‑1 (MQ=255)                                                                    gTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGAATTAACAAGGGTTGCTGGGCGGAGAAGTAAGCGTGGAACTCCCCG‑GCGGGCGGCTTTGGATCGCCCGGGAAAGGCCgggg                                                                              >  1:71647/1‑149 (MQ=255)                                                                    gTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAAGGATTCAGCAAGGGTTGCTTGCCGGAGAAGTACGCGTGGAAATCCCCG‑GCGGGCGGCTTGATATTGCCTGGGAAAGTCCcgg                                                                               >  1:779121/1‑149 (MQ=255)                                                                      ttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGGGGGACTCCCCG‑GCGGGCGGCTTGGTATTGCCTGGGAAGGGCCggg                                                                               >  1:89058/1‑147 (MQ=255)                                                                       tATGAGCTTGGGGCCAGAGAAACCCTGTCCTGCTGCATCGGCGCGTTTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTGTGCTCGCCGAAGAGTTACGCGTGGACCGCCCCG‑GCGGTCGTATTGATATCCCCTGGAAAGGTCCGGGTCa                                                                            <  2:870527/149‑1 (MQ=255)                                                                       tATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGCTGCAGGAGGGAATCAGCAAGGGTTGCTGGCCGGAGGAGCACCCGCGGGACACCCCG‑GCGGTCGTCTTTATATAGCCTGGAAAGGGCCcggggatc                                                                          >  2:1187988/1‑146 (MQ=255)                                                                          gAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGt                                                                       <  2:759969/151‑1 (MQ=255)                                                                            gCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGGACGCGTGGGACTCCCCC‑GCGGTCGTCTTTATATTCCCTGGGAAGGTCCGGGGCACCCggt                                                                      >  1:398481/1‑148 (MQ=255)                                                                              gTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTtat                                                                    <  2:1253721/150‑1 (MQ=255)                                                                              gTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTGCTGGCCGGAGCAGTACGGGTGGAACTCCCCC‑GCGGGCGTCTTGATATCGCCTTGAAAGGGCCGGGGCACCCGCtat                                                                    >  1:105565/1‑150 (MQ=255)                                                                                    taGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTCGGAGGAAAGCTCCAGCAAGGGTTTCTGGCCCGAGAAAGAACCGTGGGACTCCCCGCGCGG‑CGCGTTTGTATCGCCTGGGAAAGGCCGGGGGACCCCCtatat                                                                  >  2:832310/2‑146 (MQ=255)                                                                                     aGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTTGGGGAGGGGTTCAAAAAAGGTTGGTGGGCGAAGAAGAACGCCGGGAACTCCCCC‑GCCGGGGGCTTTATATTGCCCGGGAAGGGGCGGGGGACCCGTtatattttac                                                             >  1:117898/1‑145 (MQ=255)                                                                                         gaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGcc                                                        <  1:1023814/151‑1 (MQ=255)                                                                                          aaaCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGccc                                                       <  2:588699/151‑1 (MQ=255)                                                                                           aaCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTCAGCAAGGGTTTCTGGCCGGAGAAGTACGCGGGGGACTCCCCC‑GCGGGCGGCTTTATATCGCCTGGGAAGGGCCGGGGCACCCGTTTTATATTACttggcc                                                       >  1:449592/1‑145 (MQ=255)                                                                                           aaCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcg                                                      <  1:103315/151‑1 (MQ=255)                                                                                             cccAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcgg                                                     <  1:1184398/150‑1 (MQ=255)                                                                                                aGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGTACGCGGGGAACTCCCCG‑GCGGGCGGCTTTATATCGCCTGGAAAAGGCCGGGGCACCCGGTATATATGACTGGccccgggg                                                  >  1:694753/1‑145 (MQ=255)                                                                                                  gCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTGCGCGTGGAACTCCCCG‑GCGTTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTAc                                               <  2:79329/151‑1 (MQ=255)                                                                                                   ccTGCGGCAGCGGCGCGTGTGCTGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAGCTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGGCTGGCCCGTCGGTACa                                              <  2:463056/151‑1 (MQ=255)                                                                                                   ccTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACa                                              <  1:1211786/151‑1 (MQ=255)                                                                                                      gCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATg                                            <  2:230705/150‑1 (MQ=255)                                                                                                          cAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTac                                       <  1:453222/151‑1 (MQ=255)                                                                                                            gcggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacg                                      <  2:1168163/150‑1 (MQ=255)                                                                                                               gcgcGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACgg                                  <  1:1276050/151‑1 (MQ=255)                                                                                                                 gcgTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGGCACCCGTTATATATGACTGGCCCGGCGGGACATGTCTACGACgg                                  >  1:1209159/1‑149 (MQ=255)                                                                                                                   gtgtGCGGCGGTTTCAGGCGGGATTCAACAAGGGTTTCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGGCGTTCTTATATCGCCCGGGAAAGGCCCGGTCACCCCGTTTATATGACAGGGCCGGGGGGACAAGTCttacacggggtt                              >  2:856252/1‑140 (MQ=255)                                                                                                                   gtgtGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGAttt                              <  2:209927/151‑1 (MQ=255)                                                                                                                    tgtgCGGCGGTTGCAGTAGGGATTCACCAAGGGTTGCTTGCCGGAGAAGTACGCGGGGGACTCCCCC‑GCGGGCGTCTTGGTATCCCCTGGAAAGGGCCGGGTCACCCGGTATATAATGCTGGCCCCGCGGGACATGTCCTCGACggggtt                              >  2:714578/1‑146 (MQ=255)                                                                                                                     gtgCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTAt                            >  1:204734/1‑151 (MQ=255)                                                                                                                          gcggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATc                       <  1:311557/151‑1 (MQ=255)                                                                                                                            ggTTGCAGCAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTCGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACACGTCTAAGACGGATTTATTCctct                      >  2:681858/1‑150 (MQ=255)                                                                                                                               tGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATg                   >  2:112535/1‑150 (MQ=255)                                                                                                                                gccgTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGaa                 <  2:1140613/148‑1 (MQ=255)                                                                                                                                gCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGaa                 <  1:958686/151‑1 (MQ=255)                                                                                                                                gCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGa                  <  1:847186/150‑1 (MQ=255)                                                                                                                                gCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGa                  >  1:1278982/1‑150 (MQ=255)                                                                                                                                 cAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGa                  >  1:171121/1‑149 (MQ=255)                                                                                                                                  aGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGaa                 >  1:862226/1‑149 (MQ=255)                                                                                                                                     aGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAc            <  1:189705/151‑1 (MQ=255)                                                                                                                                       ggATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCa          <  1:94488/151‑1 (MQ=255)                                                                                                                                            cAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAgggg      <  1:393211/150‑1 (MQ=255)                                                                                                                                              gCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGgaa    <  2:89037/150‑1 (MQ=255)                                                                                                                                               cAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGgaaga  <  1:577871/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCG‑GCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGA  >  NC_000913/3995310‑3995602

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