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

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCC  >  NC_000913/3105350‑3105646                                                                                                                                                      |                                                                                                                                                    gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTa                                                                                                                                                     <  2:625309/150‑1 (MQ=255)   tgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccg                                                                                                                                                  >  1:917397/1‑151 (MQ=255)     cGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgcc                                                                                                                                                >  2:1200451/1‑151 (MQ=255)      gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccg                                                                                                                                               <  1:953354/151‑1 (MQ=255)      gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccg                                                                                                                                               >  2:450018/1‑151 (MQ=255)         tCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgccg                                                                                                                                            >  2:1024005/1‑151 (MQ=255)               cGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAAtt                                                                                                                                      >  1:384542/1‑151 (MQ=255)                gAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGTAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTc                                                                                                                                     <  1:312063/151‑1 (MQ=255)                 aaTGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTcc                                                                                                                                    >  1:1095052/1‑151 (MQ=255)                 aaTGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTcc                                                                                                                                    >  1:11317/1‑151 (MQ=255)                 aaTGATAAGATTCAGGCTGCAAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTc                                                                                                                                     <  1:1082752/150‑1 (MQ=255)                    gATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAtt                                                                                                                                 >  1:952425/1‑151 (MQ=255)                     aTAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAtt                                                                                                                                 <  1:800768/150‑1 (MQ=255)                         gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                            >  1:618196/1‑151 (MQ=255)                         gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                            <  2:302518/151‑1 (MQ=255)                          aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGc                                                                                                                           >  1:1169883/1‑151 (MQ=255)                          aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGc                                                                                                                           >  2:575713/1‑151 (MQ=255)                              aGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGcccc                                                                                                                        >  1:634963/1‑150 (MQ=255)                              aGGCTGCAAATATTATTAACCCCATGCCGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGcccc                                                                                                                        <  2:12168/150‑1 (MQ=255)                                 cTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGAc                                                                                                                    >  2:181782/1‑151 (MQ=255)                                      aaTATTATAACCCCCAGGACGCGGGCCGATTTTTACCAGACGCGGCCGCCCCGCCACCCCTCTGCGAAATCTCGCCGTTATTTTTATAAAGTAAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGAc                                                                                                                    <  2:14024/146‑1 (MQ=255)                                       atatTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCt                                                                                                              >  1:53693/1‑151 (MQ=255)                                        tattatTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                             >  2:803880/1‑151 (MQ=255)                                          ttattaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                             >  2:311694/1‑149 (MQ=255)                                           tattaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                                                                          >  2:441754/1‑151 (MQ=255)                                             ttaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgggg                                                                                                        >  1:237459/1‑151 (MQ=255)                                                 ccccaggccccggggCGATTTTTACCAGCCGCTGCCGCCCCCCCACCCCTCTGCGAAATCTCGCCTTTATCTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                     <  2:518376/138‑1 (MQ=255)                                                   ccATGACGCGGGGCGATGTTTATCAGACGCTGACGCCCCCCCACCCCTCTGCGAAATCCCGCCGTTATCTTTAAAAAGTGAATCCCGCGCAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaaa                                                                                                  <  2:1157582/151‑1 (MQ=255)                                                    cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                                                 >  2:66020/1‑151 (MQ=255)                                                      tGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGt                                                                                               >  1:1036321/1‑151 (MQ=255)                                                          gCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                                           >  2:753981/1‑151 (MQ=255)                                                           ccGGGCGATGTTTAACAGCCGCTGACCCCCGCCCACCCCTCGGCGAAATCTCGCCGTTTTTTTAAAAAAGTGAAACCCGCCAAAAACCCCTCCCGCCGCGGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                                             <  2:806421/148‑1 (MQ=255)                                                            cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                                          >  2:897392/1‑150 (MQ=255)                                                            cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                                          >  1:217745/1‑150 (MQ=255)                                                                   tttttATCAGACGCTGACGCCCCGCCACCCCTCGGCGGAATCTCGCCTTTATTTTTATAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                                     <  2:566431/146‑1 (MQ=255)                                                                     tttATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAt                                                                                >  1:521041/1‑151 (MQ=255)                                                                        ataagccGCTGGCGACCCACCAACCCGCCGCGGAACCGCGCCTTTATTTTTATAAAGTGAACCCCGCGCAAACATCCTACCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                              <  2:925963/144‑1 (MQ=255)                                                                        accAGACGCTGACGCCCCCCCTCCCCTCGGGGGATTCCCCCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAATTCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                               <  2:1036352/147‑1 (MQ=255)                                                                           agccgctgccccccccccACCCCTCGGCGAAATCCCGCCTTTATCTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                              <  2:29106/135‑1 (MQ=255)                                                                           agccGCTGCCGCCCCCCCACCCCTCTGCGAAATCTCGCCTTTATCTTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                                             <  2:716567/145‑1 (MQ=255)                                                                           aGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                                          >  2:1158491/1‑151 (MQ=255)                                                                              cGCTGCCGCCCCGCCCCCCCTCTGCGAAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                        <  1:42857/150‑1 (MQ=255)                                                                              cGCTGCCGCCCCCCCACCCCTCGGGGAAATCTCGCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                                                                       <  2:284296/151‑1 (MQ=255)                                                                                 tGACGCCCGGCCACCCCTCTGGGAATTCTCGCTGTTTTCTTTATAAAGTAAATCCGGCCCAAAATTCCTCCCGCCGCCGAAATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                        <  1:514275/147‑1 (MQ=255)                                                                                  gacgacCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                    >  2:743286/1‑150 (MQ=255)                                                                                  gacgacCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                    >  1:748547/1‑150 (MQ=255)                                                                                       ccccccACCCCTCGGCGAAATCTCCCCTTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATCCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                                                                       <  2:3134/142‑1 (MQ=255)                                                                                       cccGCCACCCCTCTGCGAAATCTCGCCGTTATTTTTATAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaaa                                                                   <  2:317761/146‑1 (MQ=255)                                                                                       cccGCCACCCCTCGGCGAATTCTCGCCTTTTTTTTTAAAAAGTAAATCCCGCCAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                              <  2:178317/151‑1 (MQ=255)                                                                                         ccccACCCCCCGGCGAAATCCCCCCTTTATCTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGaa                                                            <  2:545773/149‑1 (MQ=255)                                                                                            cATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                                         >  1:675751/1‑151 (MQ=255)                                                                                             aTCCCTCTGCGGACTCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                                        >  1:338144/1‑151 (MQ=255)                                                                                               cccTCTGCGGAATCTCGCCTTTATCTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                                         <  2:962696/148‑1 (MQ=255)                                                                                                   cggggaatccccGCCTTTTTTTTTAAAAAGTAAATCCCGCAAAAAATCCCTCCCCCCCCCAATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAat                                                  <  1:1164461/140‑1 (MQ=255)                                                                                                       ggAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                                              >  1:874725/1‑151 (MQ=255)                                                                                                       ggAATCCCGCCGTTATTTTAAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCaa                                                   <  2:598362/146‑1 (MQ=255)                                                                                                       gaaaTCTCGCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                                <  1:244579/147‑1 (MQ=255)                                                                                                       gaaaTCCCGCCTTTATTTTTAAAAAGTGAATCCCCCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                               <  2:748571/148‑1 (MQ=255)                                                                                                         aacctcGCCGTTATCTTAAAAAAGTAAATCCCGCGCAAAATCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                               <  2:806821/145‑1 (MQ=255)                                                                                                         aaTCTCGCCTTTATCTTTATAAAGTAAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                            <  2:86409/151‑1 (MQ=255)                                                                                                         aaTCTCGCCGTTATTTTTATAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                            <  2:601239/151‑1 (MQ=255)                                                                                                             gcGCTGTTATCTTTATAAAGTGAATACCGCCAAAAAAGCCTCCCGCGGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGcc                                          <  2:763858/148‑1 (MQ=255)                                                                                                              cGCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcaca                                       <  2:653370/151‑1 (MQ=255)                                                                                                              cGCCGTTATCTTTATAAAGTAAATCCCGCGAAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcaca                                       <  2:366844/151‑1 (MQ=255)                                                                                                                cctttATCTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAg                                      <  2:11418/147‑1 (MQ=255)                                                                                                                   ttttCTTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                               <  1:264766/135‑1 (MQ=255)                                                                                                                   ttATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTa                                   >  2:38648/1‑150 (MQ=255)                                                                                                                      tttttATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                 <  1:385256/147‑1 (MQ=255)                                                                                                                        tttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaa                              >  2:218197/1‑150 (MQ=255)                                                                                                                        tttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaa                              >  1:587878/1‑150 (MQ=255)                                                                                                                        tttAAAAAGTAAATCCCGCCAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaa                               <  1:755791/149‑1 (MQ=255)                                                                                                                         ttAAAAAGTGAACCCCGCGAAAAAACCCTCCCGCCGCGGATAATTCCTTTATCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACa                                <  1:321613/147‑1 (MQ=255)                                                                                                                           ataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                          >  1:779607/1‑151 (MQ=255)                                                                                                                           aaaaaGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                          <  2:458616/149‑1 (MQ=255)                                                                                                                               aGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGt                      >  2:1035546/1‑151 (MQ=255)                                                                                                                                  gAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTc                     >  2:161048/1‑149 (MQ=255)                                                                                                                                   aacaCCGCAAAAAAACCCTCCCGCCCCCGATAATTCAATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcg                   <  1:1122150/147‑1 (MQ=255)                                                                                                                                   aaTCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGt                      <  2:353143/147‑1 (MQ=255)                                                                                                                                   aaTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                  >  2:291080/1‑151 (MQ=255)                                                                                                                                   aaTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                  >  2:286089/1‑151 (MQ=255)                                                                                                                                      accacaaaaaaaTCTTCCCGCCGGATATAATTCTATTAACCGCCCCTGACGATGCTCAGATGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcg                   <  1:608408/140‑1 (MQ=255)                                                                                                                                       ccGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa               <  2:679882/150‑1 (MQ=255)                                                                                                                                            aaaaaaCCCCCCCCCCCCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTg           <  1:35469/148‑1 (MQ=255)                                                                                                                                             aaaaaCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGtt            <  2:417945/147‑1 (MQ=255)                                                                                                                                                   ccTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc  >  1:993980/1‑151 (MQ=255)                                                                                                                                                   ccTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc  >  2:829876/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCC  >  NC_000913/3105350‑3105646

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