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.8%	87.3 / 10.3	62	Y356S (TAC→TCC)	mltC	membrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (23/10);  new base C (0/29);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 2.42e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  NC_000913/3105350‑3105647                                                                                                                                                      |                                                                                                                                                     gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTa                                                                                                                                                      <  1:816062/150‑1 (MQ=255)         tCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgccg                                                                                                                                             >  2:712735/1‑151 (MQ=255)           ttttCGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATTTTTATCAGCCGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTAccgccgc                                                                                                                                               <  1:35120/147‑1 (MQ=255)            tttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTACCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGAt                                                                                                                                           <  1:995264/150‑1 (MQ=255)                       aaGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAc                                                                                                                               >  2:727163/1‑151 (MQ=255)                         gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                             >  1:239789/1‑151 (MQ=255)                         gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                             >  2:197661/1‑151 (MQ=255)                          aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGc                                                                                                                            >  2:946734/1‑151 (MQ=255)                          aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGc                                                                                                                            <  1:248276/151‑1 (MQ=255)                             cAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGcccc                                                                                                                         >  2:936496/1‑151 (MQ=255)                                     cAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                                                                  >  1:46557/1‑150 (MQ=255)                                     aaaTATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATTTTTATAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATg                                                                                                                 <  1:687109/150‑1 (MQ=255)                                      aaTATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGc                                                                                                                >  1:748983/1‑151 (MQ=255)                                        tattatTACCCCCATGCCGCCGGGCGATTTTTATCAGACGCTGACGACCCGCCATCCCTCGGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCt                                                                                                               <  1:869238/150‑1 (MQ=255)                                               aaCCCCATGCCGCCGGGCGATTTTTATCAGCCGCTGACGCCCCCCCACCCCTCTGCGGAATCTCGCCTTTATCTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                                                                           <  2:765475/147‑1 (MQ=255)                                               aaCACCACGCCGCCGGGCAATGTTTACCAGCCGCTGCCGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgggg                                                                                                         <  2:88190/149‑1 (MQ=255)                                                 ccccATGACGCCGGGCGATGTTTATCGGACGCTGACGCCCCCCCACCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                      <  1:556528/148‑1 (MQ=255)                                                 ccccAGGACGCGGGGGGATGTTTATCAGCCGCTGACGCCCCGCCACCCCTCGGCGAAATCTCGCCTTTATCTTTAAAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                      <  2:48046/148‑1 (MQ=255)                                                     aTGGCGCCGGGCGATTTTTATCAGACGCTGACGCCCCGCCTCCCCTCTGCGGATTCCCGCCTTTATTTTTAAAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                                                  <  2:332763/150‑1 (MQ=255)                                                     aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                 <  1:676413/151‑1 (MQ=255)                                                       gccGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCTAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                                                  <  2:102919/146‑1 (MQ=255)                                                         cGCCGGGCGATGTTTATCGGACGCTGCCGCCCCCCCATCCCTCTGCGGATTCTCGCCTTTATTTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGt                                                                                                <  2:693857/148‑1 (MQ=255)                                                             gggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                         >  2:849710/1‑151 (MQ=255)                                                              ggCGATGTTTATCAGACGCTGACGCCCCGCCACCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcaca                                                                                        <  1:709687/151‑1 (MQ=255)                                                                cGATTTTTATCAGACGCTGACGACCCGCCATCCCTCGGCGGAATCTCGCCGTTATCTTTAAAAAGTGAATACCGCGAAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                                                                       <  1:16855/150‑1 (MQ=255)                                                                   tGTTTATCAGACGCTGACGCCCCCCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                                                                    <  2:255120/150‑1 (MQ=255)                                                                   tGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                                     <  1:308530/149‑1 (MQ=255)                                                                     tttATCAGCCGCTAACGACCCGCCATCCCTCGGCGGAATCCCCCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                                   <  1:472671/149‑1 (MQ=255)                                                                        aTCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                                <  1:276276/149‑1 (MQ=255)                                                                               gCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                       >  2:828763/1‑151 (MQ=255)                                                                                 tGACGCCCCGCCATCCCTCTGCGAAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                                              <  1:52887/142‑1 (MQ=255)                                                                                 tGACGCCCCCCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGttt                                                                          <  1:240858/146‑1 (MQ=255)                                                                                 tGACGACCCCCCACCCCTCTGCGGAATCCCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                       <  2:468301/149‑1 (MQ=255)                                                                                  gacgccCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTAAATCCCGCGAAAAAAACCTCCCGCCCCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                                                                        <  1:644337/147‑1 (MQ=255)                                                                                   acgccCCGCCACCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                                <  2:714041/138‑1 (MQ=255)                                                                                      cccccccATCCCTCTGCGGAATCTCGCCTTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                                                                <  2:79432/150‑1 (MQ=255)                                                                                       cccGCCCCCCCTCTGCGGAATCTCGCCGTTATTTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                 <  2:874941/149‑1 (MQ=255)                                                                                       cccGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTAAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                 <  1:485690/149‑1 (MQ=255)                                                                                       cccGCCACCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAtt                                                                  <  1:394765/148‑1 (MQ=255)                                                                                        ccGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                               >  2:786102/1‑150 (MQ=255)                                                                                           ccATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATg                                                           >  1:343555/1‑151 (MQ=255)                                                                                            cATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                                          >  2:165093/1‑151 (MQ=255)                                                                                               cccTCGGCGGAATCTCCCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                       <  1:841339/151‑1 (MQ=255)                                                                                                 ctctGCGGAATCTCGCCGTTCTCTTTCTCCCGTGCCTCCCTCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                      >  1:428968/1‑150 (MQ=255)                                                                                                  tctGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa                                                     >  1:78558/1‑150 (MQ=255)                                                                                                  tctGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa                                                     >  2:32473/1‑150 (MQ=255)                                                                                                   ctGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAat                                                   >  2:37144/1‑151 (MQ=255)                                                                                                   cgccgaattctcCCCGTTATCTTTAAAAGGGAAATCCCGCGCAAAAGTCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCaa                                                    <  2:761279/142‑1 (MQ=255)                                                                                                       ggAATCTCGCCGTTATCTTTATAAAGTAAATACCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCaaa                                                   <  2:491115/147‑2 (MQ=255)                                                                                                              cGCCGTTATTTTTAAAAAGGAAATCCCGCGCAAAAACCCACCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAACCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                                            <  2:239766/147‑1 (MQ=255)                                                                                                               gCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                         >  2:457781/1‑149 (MQ=255)                                                                                                                cctttATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                                      <  1:879750/148‑1 (MQ=255)                                                                                                                    tATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACGTCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                  >  2:303695/1‑151 (MQ=255)                                                                                                                     aTCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                  >  2:709015/1‑150 (MQ=255)                                                                                                                         ttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                              >  2:347211/1‑150 (MQ=255)                                                                                                                         ttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                              >  2:578517/1‑150 (MQ=255)                                                                                                                                   aaTCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc                     <  2:770570/149‑1 (MQ=255)                                                                                                                                       ccGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa                <  1:847970/150‑1 (MQ=255)                                                                                                                                       ccGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGAGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAg               <  2:416406/151‑1 (MQ=255)                                                                                                                                       ccGCGAAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa                 <  2:224714/149‑1 (MQ=255)                                                                                                                                              aaaaCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAg        <  2:477242/151‑1 (MQ=255)                                                                                                                                                   ccTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc   <  2:309836/151‑1 (MQ=255)                                                                                                                                                     tACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATccc  >  2:476638/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  NC_000913/3105350‑3105647

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