Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	42.9%	164.1 / 11.9	85	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/43);  new base G (36/0);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.04e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133158‑3133451                                                                                                                                                 |                                                                                                                                                      aGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGc                                                                                                                                                  <  3:1070570/150‑1 (MQ=255)      ggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATcc                                                                                                                                            >  2:1243618/1‑151 (MQ=255)        cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATcc                                                                                                                                            >  2:1220683/1‑149 (MQ=255)         tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGt                                                                                                                                         >  2:1364114/1‑151 (MQ=255)          gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAg                                                                                                                                          <  2:885514/149‑1 (MQ=255)          gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGtt                                                                                                                                        <  2:431557/151‑1 (MQ=255)           cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTc                                                                                                                                       <  1:423798/151‑1 (MQ=255)            ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCCTccccgt                                                                                                                                        >  2:824074/1‑145 (MQ=255)              cgcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAg                                                                                                                                    <  3:858669/151‑1 (MQ=255)               gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCACCAGTTTaa                                                                                                                                     >  2:559198/1‑149 (MQ=255)                 ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCGGTCAAgg                                                                                                                                   >  3:484134/1‑149 (MQ=255)                   ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCt                                                                                                                               >  2:748286/1‑151 (MQ=255)                    tGTACTGTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCt                                                                                                                               <  2:205330/150‑1 (MQ=255)                    tGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGCTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCCGTCAGTTCGAgg                                                                                                                                   >  2:381420/1‑146 (MQ=255)                     gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCttt                                                                                                                             <  2:807050/151‑1 (MQ=255)                     gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCtt                                                                                                                              <  2:319443/150‑1 (MQ=255)                      tACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAAGATCttt                                                                                                                             >  2:1198860/1‑150 (MQ=255)                         tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGAGTCGCCCATCCAGGTCAAGGATCctt                                                                                                                             >  3:456127/1‑147 (MQ=255)                            tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCa                                                                                                                       <  2:816637/150‑1 (MQ=255)                            tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCCTCAGGTCTAGGATCttt                                                                                                                             >  2:738043/1‑144 (MQ=255)                            tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGATCTTTGCGcc                                                                                                                        >  3:1112085/1‑149 (MQ=255)                             aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTTCGCCAtt                                                                                                                     >  3:188807/1‑151 (MQ=255)                              aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                    <  4:276316/151‑1 (MQ=255)                              aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                    <  2:1289476/151‑1 (MQ=255)                                agtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGaa                                                                                                                  <  3:213610/151‑1 (MQ=255)                                 gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGGTCAAAGATCTTTTCGCCCTTGaa                                                                                                                  >  2:294720/1‑150 (MQ=255)                                  tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCATCCAGTTTAAAGATCctt                                                                                                                             >  2:782757/1‑138 (MQ=255)                                   ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCAACCAGTTCAAGGATCtttttgccctt                                                                                                                     >  2:201859/1‑137 (MQ=255)                                    tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                               <  1:190974/150‑1 (MQ=255)                                     gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTCAAGGATCTTTGCGCCATTGAAt                                                                                                                 >  2:536935/1‑147 (MQ=255)                                      agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGt                                                                                                               >  1:317258/1‑148 (MQ=255)                                       gtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGAGGGTGGCCCATCCATTTCAAAGAGCTTTTCGCCGTTGTATGtgtgt                                                                                                           >  3:584792/1‑147 (MQ=255)                                          tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTg                                                                                                         <  1:381094/150‑1 (MQ=255)                                               tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGGTCAAAGATCTTTGTGCCATTGAATGTCAGTTgggggg                                                                                                    >  3:629628/1‑148 (MQ=255)                                                gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGATCAAGGAAATTTTCGCCCTTTAATGGCAGTTgcggcg                                                                                                    >  2:376234/1‑149 (MQ=255)                                                gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTa                                                                                                  <  1:283861/151‑1 (MQ=255)                                                 cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                                                 <  2:719708/151‑1 (MQ=255)                                                    ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACg                                                                                               <  4:54272/150‑1 (MQ=255)                                                        tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGATCTTTTTGCCCATGGATGTCCGGTTCGGCGGAACGGCg                                                                                            >  4:377631/1‑149 (MQ=255)                                                        tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTCAAGGATCTTTTCGGCATTGAATGTTCGTTGGGGCGTAACgg                                                                                              >  3:100177/1‑147 (MQ=255)                                                        tCGGTAAATCATCCCCTAACCATTTTGGTGATCAGGAAAGCGCGGCGTGAAATGCGCGTTGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGa                                                                                           <  3:96179/150‑1 (MQ=255)                                                                 catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                 <  4:220613/151‑1 (MQ=255)                                                                  atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCCCCCCGGTCAAGGATCTTTGGGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTg                                                                                 >  3:807698/1‑150 (MQ=255)                                                                  atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGCGTAACGGGGAttttttt                                                                                    >  4:360322/1‑146 (MQ=255)                                                                  atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCAGTTCAAGGATCTTTTCGCCCTTTAATTTCAGTTTCGGCGGAACGGCGGTTTTTTTCttt                                                                                >  3:190696/1‑149 (MQ=255)                                                                     cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTCAAGGATCTTTTCGCCATTTGATGTCCGTTGTGGGGGAAAGGCGGGTTTTTTCTGTaa                                                                              >  2:154736/1‑150 (MQ=255)                                                                         aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCt                                                                          <  1:76653/150‑1 (MQ=255)                                                                         aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTg                                                                         <  3:1218102/151‑1 (MQ=255)                                                                          tCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTCAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGCGGAAAGGGGGTTTTTTTCTGTTAAGCTg                                                                         >  1:319073/1‑150 (MQ=255)                                                                          tCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTTAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGGGGAACGGCGGTTTTTTTCTTTTACGc                                                                           >  2:114297/1‑148 (MQ=255)                                                                           ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGc                                                                        <  3:260641/150‑1 (MQ=255)                                                                               tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGt                                                                    <  2:1112032/150‑1 (MQ=255)                                                                               tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGt                                                                    <  3:782792/150‑1 (MQ=255)                                                                                ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCCTTCAGGTTCAGGGTCTTTTCGCCCTTTGATGTCAGGTGCGGGGGAACGGCGGATTTTTTCTGTTACGCTGccc                                                                      >  3:764932/1‑146 (MQ=255)                                                                                   gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAGGTCAAAGATTTTTTTGCCCTTTGATGTCCAATGCCGCGGAACGGCGGTTTTTTTCTGTACAGCTGCCCGtt                                                                   >  3:979834/1‑147 (MQ=255)                                                                                     gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTTTTACGCTGCCCGtttgggg                                                              >  3:1086448/1‑145 (MQ=255)                                                                                      aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTc                                                            <  2:1219550/151‑1 (MQ=255)                                                                                        cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCa                                                           <  4:308954/150‑1 (MQ=255)                                                                                          ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCGGTCAAGGATCTTTTCGCCATTGTATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAACGCTGcccggttggggtccgg                                                         >  2:132764/1‑136 (MQ=255)                                                                                           gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCATCCCGGTCAAAGGTCTTTTCGCCCTTGAATGTCAGGTTGGGGGGAACGGGGGATTTTTTCTGTTAAACTTCCCGTtggtggggggg                                                         >  2:1246130/1‑144 (MQ=255)                                                                                             aaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTTAAGGAGCTTTTTGGCCTTGGATGTCTGTTGCGGCGGAAAGGCGGTTTTTTTCTCTTACGCTGCCCGGtggtggggagggtt                                                      >  3:756709/1‑142 (MQ=255)                                                                                              aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAc                                                    <  2:1164586/151‑1 (MQ=255)                                                                                               gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAcg                                                   <  3:748318/151‑1 (MQ=255)                                                                                                 gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGGTCTTTTCGCCCTTTAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTTACGCTGCCCGTtgggggggagggtt                                                      >  1:23345/1‑138 (MQ=255)                                                                                                  cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACgcgc                                                <  3:1198919/151‑1 (MQ=255)                                                                                                   ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACgcgc                                                <  1:345200/150‑1 (MQ=255)                                                                                                     cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTCTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                              <  3:381456/150‑1 (MQ=255)                                                                                                       tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTCCTGTAACGCTGCCCGTTGGGGGGCAGGTTTACGCGCaaa                                             >  1:335151/1‑149 (MQ=255)                                                                                                               gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGt                                   <  3:1243677/151‑1 (MQ=255)                                                                                                                cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGGTTTTTTCTGTAAAGCTGCCAGTTGGTGGGCAGGGTTTCGCGCAAAAGCGGgttgtt                                  >  3:313887/1‑151 (MQ=255)                                                                                                                    gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcat                               <  3:275894/150‑1 (MQ=255)                                                                                                                     tGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCACCCAGTTCAAGGATCTTTGCGCCATTGGATGTCAGTTGCGGGGTAACGGCGATTTTTTTCTGTAACGCCGCCCGTTTGTGGGCAGGGTTTAGCGCAAAAGCGGGTTGTttat                               >  3:1105288/1‑149 (MQ=255)                                                                                                                      gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                            <  2:188780/151‑1 (MQ=255)                                                                                                                       aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAAt                           <  1:19668/151‑1 (MQ=255)                                                                                                                        cGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCCTCCAGTTCAAGGGTCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAAAGGCGATTTTTTTCTGTAACGCTGCCCGTTGGTGGGCAGGGTTACGCGCAAACGCGGGTTGTTTTTCaa                            >  2:490698/1‑149 (MQ=255)                                                                                                                            tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                       <  1:377635/150‑1 (MQ=255)                                                                                                                            tGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTTTATGTCAGTTGCGGCGGAACGGCGCTTTTTTTCTTTAACGCTGCCAGTTGGGGGGCAGGGTTTCGCGCAAACGCGGGTTGTTTATCAATGcc                        >  3:65491/1‑149 (MQ=255)                                                                                                                             gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCaa                      <  3:559229/150‑1 (MQ=255)                                                                                                                             gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                     <  3:154756/151‑1 (MQ=255)                                                                                                                                     tCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATc              <  3:376268/150‑1 (MQ=255)                                                                                                                                          ggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt        <  1:469651/151‑1 (MQ=255)                                                                                                                                           gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATa       <  1:75716/151‑1 (MQ=255)                                                                                                                                            cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc      >  4:167607/1‑151 (MQ=255)                                                                                                                                              gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg    <  2:764897/151‑1 (MQ=255)                                                                                                                                              gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg    <  1:106504/151‑1 (MQ=255)                                                                                                                                                gtgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGc  <  3:147003/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133158‑3133451

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: ‑