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

TGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCAT  >  NC_000913/3105352‑3105649                                                                                                                                                     |                                                                                                                                                       tgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccg                                                                                                                                                     >  1:372698/1‑151 (MQ=255)  gcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgc                                                                                                                                                    >  2:277318/1‑151 (MQ=255)  gcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccg                                                                                                                                                     <  1:774487/150‑1 (MQ=255)   cGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgcc                                                                                                                                                   >  1:197978/1‑151 (MQ=255)       tCTTTTCGAATGATAAGATTCAGCCTGCCAATATTATTAACCCCATGCCGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgcc                                                                                                                                                <  1:7944/150‑1 (MQ=255)          tttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgccg                                                                                                                                               <  1:89716/148‑1 (MQ=255)          tttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgccg                                                                                                                                               <  2:157659/148‑1 (MQ=255)          tttCGAATGATAAGATTCAGCCGGCAAATATTATTACCCCCAGGCCGCCGGGCGATGTTTATCAGACGC‑TGACGCCCCGCCACCCCTCTGCGAAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTACCGCCGCCGAt                                                                                                                                             <  1:71136/150‑1 (MQ=255)             cGAATGATAAGTTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAAtt                                                                                                                                         >  2:314130/1‑151 (MQ=255)              gAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAAtt                                                                                                                                         <  2:534869/150‑1 (MQ=255)                aTGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCa                                                                                                                                      <  1:340888/151‑1 (MQ=255)                  gATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGGATCTCGCCGTTATCTTTATAAAGAGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAt                                                                                                                                     >  1:330104/1‑150 (MQ=255)                        aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                               <  2:56674/150‑1 (MQ=255)                           cAGGCTGCCATTATTATTAACCCCATGACGCCGGGCGTTGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGcc                                                                                                                             <  1:482947/149‑1 (MQ=255)                           cAGGCGGCCAATATAATTAACCCCATGGCGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCCCTGCGAAATCTCCCCGTTATCTTTATAAAGTGAATCCCGCGCAAAATTCCTCCCGCCGCCGATAATTCCATTAACCGcc                                                                                                                             <  1:1181428/149‑1 (MQ=255)                                tGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACg                                                                                                                      >  2:206942/1‑151 (MQ=255)                                  ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGa                                                                                                                     >  2:489280/1‑150 (MQ=255)                                          attaACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAggg                                                                                                            >  1:598520/1‑151 (MQ=255)                                             aaCCCCATGCCGCCGGGCGTTGTTTATCAGACGC‑TGCCGCCCCGCCACCCCTCGGCGGAATCTCGCCTTTATTTTTATAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgggg                                                                                                           <  2:831043/149‑1 (MQ=255)                                             aaCCCCATGCCGCCGGGCGATTTTTATCAGCCGC‑TGACGCCCCGCCACCCCCCGGCGGAATCTCCCCTTTATTTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCa                                                                                                         <  2:139048/151‑1 (MQ=255)                                                acCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                                       >  1:151932/1‑150 (MQ=255)                                                 ccATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                                      >  1:853459/1‑150 (MQ=255)                                                   atacggcggggCGATTTTTATCAGCCTC‑AGACGCCCCGCCTCCCCTCGGCGAATCCCCGCTTTTATTTTTATAAGGTGAATACCCCGAAAAATCCCTCCCGCCGCAAATAATCCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                   <  1:902443/143‑1 (MQ=255)                                                   aTGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                   >  2:382602/1‑151 (MQ=255)                                                        gCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                                              >  1:394678/1‑151 (MQ=255)                                                         ccGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                                              >  1:923688/1‑150 (MQ=255)                                                           gggCGTTTTTTACCAGCCGC‑TGACGCCCCCCCATCCCTCTGGGAAATCTCGCCGTTATTTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                                                               <  2:961450/147‑1 (MQ=255)                                                           gggCGATGTTTA‑CAAACCCTTGACAACCCCCCACCCCTCTGCGGACTCCCCCGGTTATCTTTAAAAAGTAAATACCGCGCAAAAATCCTCCCCCCCCCGATAATTCCTTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGtt                                                                                                 <  2:290867/145‑1 (MQ=255)                                                               gATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                                        >  2:1066389/1‑150 (MQ=255)                                                                 tttttATCAGACGC‑TGACGCCCCCCCACCCCTCTGGGAAATCCCGCCTTTATTTTTATAAAGTGAATCCCGCGCAAAATCCCTCCCCCCGCCAATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                                       <  2:522125/147‑1 (MQ=255)                                                                  gTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                                     >  1:1069154/1‑150 (MQ=255)                                                                  gTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                                     >  1:365792/1‑150 (MQ=255)                                                                   tttATCGGCCGT‑TGAGTACCCCCCTTCCCTCTGGGAAATCACGCCGTTCTCTTTAAAAAGTAAATCCCCCGAAAAACTCTTCCCGCCGCCGAAAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACaa                                                                                    <  2:131120/150‑1 (MQ=255)                                                                   tttACCAGCCGC‑TGACGCCCCGCCACCCCCCTGCGGAATCCCGCCGTTATCTTTAAAAAGGAAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                                                                         <  2:753611/145‑1 (MQ=255)                                                                      accggACGC‑TGACGCCCCGCCACCCCTCGGCGAAATCTCGCCGTTTTTTTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                                  <  2:1184171/145‑1 (MQ=255)                                                                      accAGCCGC‑TGACGCCCCCCCACCCCTCGGGGAATTCTCGCCTTTATTTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                                 <  1:313648/148‑1 (MQ=255)                                                                       tCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                               >  1:182297/1‑151 (MQ=255)                                                                                tgccgcccccccACCCCTCTGCGAATCCTCGCCTTTATCTTTAAAAGGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGttt                                                                            <  2:825197/140‑1 (MQ=255)                                                                                 gcctcccccccACCCCTCTGCGGAACCCCGCCGTTATCTTTAAAAAGTAAATCCCGCGAAAAACCCTTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaaa                                                                      <  2:184851/146‑1 (MQ=255)                                                                                   cgcccccccACCCCTCTGCGAAACCTCCCCTTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCCCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                         <  2:1009274/143‑1 (MQ=255)                                                                                   cgacCCGCCACCCCTCTGCGGAATCCCGCCTTTATCTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGttt                                                                            <  2:1054998/143‑1 (MQ=255)                                                                                    gacCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                   >  2:277143/1‑151 (MQ=255)                                                                                      ccccccACCCTTCGGCGAATTCGCGCCGTTATTTTTAAAAGGTGAATCCCGCAAAAAAACCCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                                                                  <  1:178728/150‑1 (MQ=255)                                                                                      cccGCCACCCCTGTGGGGATTCCCGCGTTTATCTTTAAAATGTAAATCCCGGCAAAAAACCCTCCCCCCGCCAATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                   <  2:1194680/149‑1 (MQ=255)                                                                                        cGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                                                >  1:51541/1‑150 (MQ=255)                                                                                        cGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                                                >  1:199843/1‑150 (MQ=255)                                                                                          ccATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATg                                                             >  2:1157134/1‑151 (MQ=255)                                                                                              cccTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                         >  1:759129/1‑151 (MQ=255)                                                                                              cccTCGGCGAAATCTCGCCTTTATTTTTAAAAAGTAAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                         <  2:554025/151‑1 (MQ=255)                                                                                               ccccggCGAAATCTCGCGGTTTTCTTAATAAAGTAAATCCCGCGCAAAATTCCTCCCGCCGCCGAAAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                        <  2:784534/146‑1 (MQ=255)                                                                                                ctctGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                        >  2:664353/1‑150 (MQ=255)                                                                                                   tggggACTCTCGCCTTTATCTTTAAAAGTTGAATCCCCCGAAAAAATCTTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa                                                       <  1:478670/145‑1 (MQ=255)                                                                                                    gCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                                   <  2:538719/151‑1 (MQ=255)                                                                                                      ggAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                                                 >  2:750782/1‑151 (MQ=255)                                                                                                      ggAATCGCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACACTTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                                                 >  1:40746/1‑151 (MQ=255)                                                                                                      gaaaTCCCGCCTTTATTTTTAAAAAGTAAATCCCGCCCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCaa                                                      <  2:242713/144‑1 (MQ=255)                                                                                                        aaTCTCGCCTTTTTTTTTAAAAAGGAAATCCCGCGAAAAATTCCTCCCGCCGCCGATAATTCCTTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                                  <  1:1043072/148‑1 (MQ=255)                                                                                                        aaTCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                               >  2:936914/1‑151 (MQ=255)                                                                                                        aaTCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                               >  1:612940/1‑151 (MQ=255)                                                                                                             cGCCGTTATTTTTAAAAAGTAAACCCCCCGAAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcaca                                          <  1:588040/151‑1 (MQ=255)                                                                                                             cGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                           >  2:220785/1‑150 (MQ=255)                                                                                                               cgtttATTTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAg                                         <  1:540643/147‑1 (MQ=255)                                                                                                                   tATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                    >  2:738208/1‑151 (MQ=255)                                                                                                                          ataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                             >  2:679464/1‑151 (MQ=255)                                                                                                                          ataAAGTAAATCCCGCACAAAAACCCTCCCGCCCCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaaa                               <  2:912094/149‑1 (MQ=255)                                                                                                                           taAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                            >  1:50144/1‑151 (MQ=255)                                                                                                                            aaaTTGAATACGGCGAAAAACTCCTCCCGCCGCCGATAATTCAATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCa                           <  2:330111/151‑1 (MQ=255)                                                                                                                            aaaGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                            <  2:1011102/150‑1 (MQ=255)                                                                                                                             aagtaaatcccccaaaaaaaCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAg                          <  2:823717/136‑1 (MQ=255)                                                                                                                                  aaTCCCGCGAAAAATCCCCCCCGCCGCGGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                     <  1:729300/151‑1 (MQ=255)                                                                                                                                      ccGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAg                 <  2:691135/151‑1 (MQ=255)                                                                                                                                      ccGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa                  <  2:612014/150‑1 (MQ=255)                                                                                                                                      ccGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAg                 <  2:575880/151‑1 (MQ=255)                                                                                                                                                     aCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCcat  >  2:562096/1‑151 (MQ=255)                                                                                                                                                     aCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCcat  >  1:366629/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       TGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGC‑TGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCCAT  >  NC_000913/3105352‑3105649

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