Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,105,499	0	A	C	40.9%	121.0 / 12.0	66	Y356S (TAC→TCC)	mltC	membrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (31/8);  new base C (0/27);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAAT  >  NC_000913/3105350‑3105644                                                                                                                                                         |                                                                                                                                                    gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAc                                                                                                                                                    <  1:131822/151‑1 (MQ=255)  ctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAcc                                                                                                                                                   >  2:902496/1‑151 (MQ=255)  ctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAc                                                                                                                                                    <  2:129451/150‑1 (MQ=255)   tgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAccg                                                                                                                                                  >  1:748108/1‑151 (MQ=255)   tgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAccg                                                                                                                                                  >  1:257880/1‑151 (MQ=255)        gTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTAccgccgcc                                                                                                                                             <  1:3707/151‑1 (MQ=255)             ttCGAATGAAAAGATTCAGGCTGCCAATATTATAAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATCCCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATa                                                                                                                                       <  2:995897/150‑1 (MQ=255)              tCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAAt                                                                                                                                     >  1:394953/1‑151 (MQ=255)              tCGAATGATAAGATTCAGGCTGCAAATATTATTAACCCCATGCCGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATaa                                                                                                                                      <  2:326830/150‑1 (MQ=255)               ctaaTGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAAtt                                                                                                                                    >  1:379433/3‑151 (MQ=255)                     aTAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCAtt                                                                                                                               >  1:1003271/1‑150 (MQ=255)                     aTAAGATTCAGGCTCCCAATATTATTACCACCAGTACGCCGGGTGTTGTTTAGCAGACGCTAACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAATTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCa                                                                                                                                 <  2:626301/148‑1 (MQ=255)                               ggCTCCAAATATTATAACCCCCATGCGGCCGGGCGATGTTTACCAGCCGCTAACGCCCCGCCACCCCTCTGCGAATT‑CT‑CGCCGTTATCTTTATAAAGTAAATCCCGCGA‑AAAAATCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGcccc                                                                                                                      <  2:835883/149‑1 (MQ=255)                                       atatTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGc                                                                                                             >  2:278254/1‑150 (MQ=255)                                         attattAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                           >  1:1010421/1‑150 (MQ=255)                                           tattaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                                                                        >  1:333725/1‑151 (MQ=255)                                               aaCCCCATGCCGCGGGGGGATGTTTATCAGCCGCTGCCGCCCCGCCACCCCTCTGCGGAAT‑CT‑CGCCGTTATTTTTATAAAGTAAATCCCGCGC‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGc                                                                                                     <  1:988891/150‑1 (MQ=255)                                                acacCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCa                                                                                                    >  2:587405/1‑150 (MQ=255)                                                  acCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                                  >  1:32641/1‑150 (MQ=255)                                                  acCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                                  >  1:333642/1‑150 (MQ=255)                                                  acCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                                  >  1:374276/1‑150 (MQ=255)                                                             gggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                      >  2:88773/1‑151 (MQ=255)                                                                  aTGTTTATCAGACGCTGAAGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                                                                 >  2:950997/1‑151 (MQ=255)                                                                          cAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                         >  2:12786/1‑151 (MQ=255)                                                                           agccgtgaccgcccccccACCCCTCTGCGAAAT‑CC‑CGCCGTTATTTTTAAAAAGTAAATCCCGCGC‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                          <  1:721707/137‑1 (MQ=255)                                                                           aGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                                        >  1:484235/1‑151 (MQ=255)                                                                              cgttaccccccgcccACCCCTCTGCGAAAT‑CC‑CGCCTTTTTTTTAAAAAAGTAAATCCCGCGA‑AAAAACCCTCCCCCCGCG‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                          <  1:45840/133‑1 (MQ=255)                                                                                cggccgcccccccACCCCTC‑GGGGAATCCT‑CGCCTTTATTTTAAAAAAGGAAATCCCGCAA‑AAAACCCCTCCCCCCGCA‑G‑ATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                                                                   <  2:6847/144‑1 (MQ=255)                                                                                 tgccgccccGCCACCCCTCGGCGAAAT‑CT‑CGCCTTTATTTTTATAAGGTAAATCCCGCGA‑AAAAATCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                         <  1:909487/138‑1 (MQ=255)                                                                                 tgaccccccGCCCCCCCTCGGCGAAAT‑CC‑CCCCTTTTTTTTTAAAAAGTGAACCCCGCGA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                            <  2:885213/135‑1 (MQ=255)                                                                                 tGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                                                                   >  1:304300/1‑150 (MQ=255)                                                                                     gacCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                              >  2:979879/1‑151 (MQ=255)                                                                                      cccccccACCCCTC‑GGGGAAT‑CCCCGCTGTTTTTTTTAAAAAGTAAATCCCGCAA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                      <  2:620791/141‑1 (MQ=255)                                                                                        cccccACCCCTCTGCGGAAC‑CT‑CGCCGTTATCTTTATAAAGTGAATCCCGCGC‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                                           <  1:118100/148‑1 (MQ=255)                                                                                        ccGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                            >  2:837211/1‑150 (MQ=255)                                                                                           ccACCCCTCTGCGGAAT‑CT‑CGCCTTTATTTTTAAAAAGTGAATCCCGCGA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                         <  2:248994/150‑1 (MQ=255)                                                                                           cacccccTCTGCGAATT‑CT‑CGCCGTTTTTTTAAAAAAGTAAATCCCCCGA‑AAAAACCCACCCCCCGCC‑G‑ATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATg                                                        <  2:892325/147‑1 (MQ=255)                                                                                                ccTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                   >  2:568774/1‑151 (MQ=255)                                                                                                ccTCCGGGGAATCCT‑CGCCTTTTTTTTTAAAAAGGGAATCCCGCGA‑AAAAACCCTCCCGCCGCC‑G‑ATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                    <  2:458422/151‑1 (MQ=255)                                                                                                      cGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                             >  1:329563/1‑151 (MQ=255)                                                                                                         gaatccccGCCGTTATTTTTAAAAAGTAAATCCCGCAA‑AAAATCCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa                                                  <  2:1176006/138‑1 (MQ=255)                                                                                                         aaT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                          >  2:570005/1‑151 (MQ=255)                                                                                                           t‑ct‑cGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGcc                                        >  1:1149645/1‑151 (MQ=255)                                                                                                              ctcccgtttATTTTTAAAAAGTAAATCCCGCAA‑AAAAACCCTCCCGCCCCC‑T‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                          <  2:649813/141‑1 (MQ=255)                                                                                                                  cctttATCTTTAAAAAGTGAATCCCGCGA‑AAAATCCCTCCCCCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                      <  2:782522/145‑1 (MQ=255)                                                                                                                  ccGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                                   >  2:264881/1‑151 (MQ=255)                                                                                                                   ctgtttttttAAAAAG‑GAAAACCCCGCGAAAAACCCCACCCCCCCC‑GTATAATTCCCTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                                         <  2:1114029/138‑1 (MQ=255)                                                                                                                     ttATTTTTAAAAATTAAATACGGCAA‑AAAAACCCTCCCCCCGCC‑G‑ATAATTCAATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                          <  1:944886/141‑1 (MQ=255)                                                                                                                     ttATTTTTAAAAAGTAAATCCCGCCA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTa                                 <  1:419058/150‑1 (MQ=255)                                                                                                                     ttATCTTTAAAAAGTAAATCCCGCGA‑AAAAACCCTACCGCCGCC‑A‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTa                                 <  1:61726/150‑1 (MQ=255)                                                                                                                         ttttAAAAATTAAATCCCGCGA‑AAAACCCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTaa                                <  2:257875/146‑1 (MQ=255)                                                                                                                         cTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaa                             >  1:862119/1‑150 (MQ=255)                                                                                                                           ttAAAAAGTAAATCCCGCGA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                           <  1:595537/150‑1 (MQ=255)                                                                                                                           ttAAAAAGAGAACCCCGCAA‑AAAAATCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                           <  2:466935/150‑1 (MQ=255)                                                                                                                             aaaaaGTAAATCCCGCGA‑AAAAACCCCCCCCCCGCC‑A‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                           <  1:848948/146‑1 (MQ=255)                                                                                                                               aaaGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                       >  2:428999/1‑150 (MQ=255)                                                                                                                                aaGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                       >  1:1174436/1‑149 (MQ=255)                                                                                                                                     aaTCCCGCGA‑AAAAACCCTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGt                    <  1:563877/147‑1 (MQ=255)                                                                                                                                        cccGCGG‑AAAAAACCCCCCCCCCCCAG‑ATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCCCAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg               <  2:849812/149‑1 (MQ=255)                                                                                                                                         cccgcaaaaaaaCCCTCCCCCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa              <  1:1171877/143‑1 (MQ=255)                                                                                                                                         ccGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa              >  1:1085785/1‑149 (MQ=255)                                                                                                                                         ccGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa              >  2:1034793/1‑149 (MQ=255)                                                                                                                                              c‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt       >  1:1005306/1‑151 (MQ=255)                                                                                                                                               aaaaaccccTCCCCCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGa        <  2:1010370/143‑1 (MQ=255)                                                                                                                                                aaaaaTCCTCCCCCCGCC‑G‑ATATTCCCTTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt       <  1:1177289/150‑1 (MQ=255)                                                                                                                                                      ccTCCCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAAt  <  1:989380/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAAT‑CT‑CGCCGTTATCTTTATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCC‑G‑ATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAAT  >  NC_000913/3105350‑3105644

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