Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,457,105	0	A	G	46.2%	82.5 / 22.1	65	I441V (ATC→GTC)	mpl	UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (6/29);  new base G (29/1);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 2.52e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456956‑4457254                                                                                                                                                      |                                                                                                                                                      atatTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGAt                                                                                                                                                      <  1:592935/151‑1 (MQ=255)     tCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCa                                                                                                                                                   <  1:853968/150‑1 (MQ=255)            caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGaaa                                                                                                                                            >  1:844799/1‑150 (MQ=255)            caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAc                                                                                                                                           <  1:201461/151‑1 (MQ=255)            caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTccct                                                                                                                                                  >  2:883101/1‑142 (MQ=255)             aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                          <  2:26084/151‑1 (MQ=255)              ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCAACAGaaa                                                                                                                                            >  1:628311/1‑148 (MQ=255)                    agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTgg                                                                                                                                     >  1:451398/1‑149 (MQ=255)                     gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTgg                                                                                                                                     >  2:432663/1‑148 (MQ=255)                      aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGTCATCAGAAAAAGCTGGGTgg                                                                                                                                 >  1:593568/1‑151 (MQ=255)                       agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCtgggtggg                                                                                                                                >  2:324010/1‑151 (MQ=255)                       agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATggg                                                                                                                                >  2:393913/1‑151 (MQ=255)                          ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCt                                                                                                                              <  2:86207/150‑1 (MQ=255)                           gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGggtggggtgg                                                                                                                            >  2:385639/1‑151 (MQ=255)                                  ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                       >  2:921703/1‑149 (MQ=255)                                    tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAAGAACTGCTGGGAGGGGTGGGggagaa                                                                                                                     >  1:593903/1‑149 (MQ=255)                                     gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGAAGAAgg                                                                                                                   >  2:288706/1‑150 (MQ=255)                                         tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGGGGGGTCCATCAGAAACCGAAGGAGGGGGGGGCGAAGAAAGCgg                                                                                                                >  1:177902/1‑149 (MQ=255)                                          cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                              <  1:384623/150‑1 (MQ=255)                                             ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcg                                                                                                          <  1:883072/151‑1 (MQ=255)                                                       ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAtt                                                                                                <  1:44852/151‑1 (MQ=255)                                                             gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTCGGcc                                                                                           >  2:886567/1‑150 (MQ=255)                                                             gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGGCGCGCCGGAAATCGGcc                                                                                           >  1:914504/1‑150 (MQ=255)                                                               gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGAAGAAAGGGGGAGCCGCGCAGGAATTTGGCCtt                                                                                         >  2:42169/1‑150 (MQ=255)                                                               gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTTGGccct                                                                                         >  2:410536/1‑148 (MQ=255)                                                                 tgtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGGGGGAGCCGCGCAGTAATTCGGCCTTAGc                                                                                      >  2:761354/1‑151 (MQ=255)                                                                        aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                                <  2:404634/150‑1 (MQ=255)                                                                         cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGAtata                                                                              <  2:501113/151‑1 (MQ=255)                                                                            tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGAAGAAGGCGGGAGACGCGCAGTAATTCGGCCTTTGCCAGATATAGt                                                                            >  1:18967/1‑150 (MQ=255)                                                                             ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                            <  1:762189/149‑1 (MQ=255)                                                                             ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                          <  2:156725/151‑1 (MQ=255)                                                                             ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                          <  1:484587/151‑1 (MQ=255)                                                                               cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGc                                                                        <  2:397251/151‑1 (MQ=255)                                                                               cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGGGAAGAAGGGGGAAAACGCGCAGTAAATTGGCCCTAGGCAGATATAATATgg                                                                        >  2:494291/1‑150 (MQ=255)                                                                                aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                       >  2:880998/1‑151 (MQ=255)                                                                                 gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGGGAAGAAGGCGGGAGGCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGc                                                                        >  2:251057/1‑149 (MQ=255)                                                                                   tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGGGGAGAAAGGGGGAACCGCGCCGGAATTCGGCCTTAGCCCGATATAGAATg                                                                         >  1:882193/1‑146 (MQ=255)                                                                                    atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAg                                                                    <  2:564532/150‑1 (MQ=255)                                                                                      ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                                  <  2:451411/150‑1 (MQ=255)                                                                                       gtggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCg                                                                <  1:432651/151‑1 (MQ=255)                                                                                         ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGa                                                              <  2:116323/151‑1 (MQ=255)                                                                                            gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATaa                                                           <  2:317722/151‑1 (MQ=255)                                                                                             aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                          <  1:46050/151‑1 (MQ=255)                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                        <  1:186238/151‑1 (MQ=255)                                                                                                  cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                     <  2:727323/151‑1 (MQ=255)                                                                                                   gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGGAAgg                                                         >  2:37501/1‑146 (MQ=255)                                                                                                   gCTCAGCCTGGCGACAATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGGAGGGCTGGCGAAAAAGGCGGAAGCCGCGCCGTAATTCGGCCCTAGCCAGATATAGTAAGCGGAGGCCGGGGAAGGCggt                                                     >  2:630245/1‑148 (MQ=255)                                                                                                     tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGACAGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                     >  2:351306/1‑146 (MQ=255)                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                               <  1:921674/151‑1 (MQ=255)                                                                                                            ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                           <  1:154734/151‑1 (MQ=255)                                                                                                              cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                         <  2:837257/151‑1 (MQ=255)                                                                                                                aCCATATTCTGTTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAGCTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGTTAAGGCGTTCACGCCGCATc                                          <  1:852493/148‑1 (MQ=255)                                                                                                                      ttCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATTAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGAAGCCGCGCCGTAATTTGGCCTTTGCCAGATTTAGTATTCGTAAGCCGGGTAAGGGGTTCACGCCGCATCCGGCAtttt                                 >  2:443295/1‑150 (MQ=255)                                                                                                                        cTGGTGATGAGCAACGGCGGCTTTGGGGGGATCCACCCGAAACTGCCGGATGGGGCGGCGAAGAAGGCGGAAGACGGGCAGTAATTCGGGCCTAGCCAGATATAGTATGCGGAGGCCCGGGAAGGGGGTCACGCCGCATCCGGCAttt                                  >  2:556654/1‑148 (MQ=255)                                                                                                                          ggTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGGGGAGAAGGGGGGAGCCGCGCAGTAATTCCGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTAGc                              >  1:279404/1‑150 (MQ=39)                                                                                                                           gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCAt                            <  2:593590/151‑1 (MQ=255)                                                                                                                              acgggaagcgggggTTTTGGAGGGGACCATCAGAAACTGCTGGATAGACTGGCGAAGATGGCGGAGGCCGCGCAGAAATTCGCCCCTACCCAGATATATTATGCCTAGGCCGGCTAAGGCGTTAACGCCGCATCCGGCATTTGAGCCTa                           <  1:788245/141‑1 (MQ=255)                                                                                                                               tgaGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGGGGAGAAGGGGGGAGCCGCGCCGGAATTTGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCCCCTCCCGCCTTTTAGCATAATg                        >  1:471819/1‑151 (MQ=255)                                                                                                                               tgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGt                         <  1:351791/150‑1 (MQ=255)                                                                                                                                gaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGc                       <  2:222865/151‑1 (MQ=255)                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  1:34984/1‑151 (MQ=255)                                                                                                                                                ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg        >  2:291023/1‑150 (MQ=255)                                                                                                                                                   gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc    >  1:153532/1‑151 (MQ=255)                                                                                                                                                    ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCg   >  1:373360/1‑151 (MQ=255)                                                                                                                                                     gATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt  >  2:106704/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456956‑4457254

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