Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,105,489	0	C	A	43.5%	97.9 / 33.1	62	Q353K (CAA→AAA)	mltC	membrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base C (30/5);  new base A (0/27);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGC  >  NC_000913/3105343‑3105630                                                                                                                                                    |                                                                                                                                               gcagCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑aaaa                                                                                                                                           <  2:518898/151‑1 (MQ=255) gcagCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑aaa                                                                                                                                            >  1:119169/1‑150 (MQ=255)     cGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATcc                                                                                                                                       <  1:171653/151‑1 (MQ=255)                 cTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTAccgccgccg                                                                                                                            >  2:410609/1‑150 (MQ=255)                  gtttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGAGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGAt                                                                                                                          >  1:999921/2‑151 (MQ=255)                   tttCGAATGATAAGATTCAGGCTGCCAATATTATTACCACCATGCCGCCGGGCGATTTTTATCAGACGCTGACGACCCGCCACCCCTCGGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGA‑AAAAATCCTACCGCCGCCGAt                                                                                                                          <  2:551785/150‑1 (MQ=255)                           gATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCAtt                                                                                                                 <  2:661432/151‑1 (MQ=255)                              aaGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAAc                                                                                                              >  1:1081570/1‑151 (MQ=255)                                  ttCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGcc                                                                                                          >  1:917149/1‑151 (MQ=255)                                    cAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGcccc                                                                                                        >  2:698572/1‑151 (MQ=255)                                      ggCTGCAAATTTTATAAACACCATGCCCCCGGGCGATGTTTATCAGACCCTGACGCCCCCCCATCCCTCTGCGGAATCTCCCCGTTATCTTT‑ATAAAGTAAATACCGCGA‑AAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                      <  2:307834/151‑1 (MQ=255)                                      ggCTGCAAATATTATTAACCCCAGGACGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCACCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                      <  2:1053505/151‑1 (MQ=255)                                       gCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                      >  1:467424/1‑150 (MQ=255)                                             aattttttaaaCCCCAGGACGCGGGGGGATTTTTACCAACCTCTGCCGCCCCCCCCCCCCTCTGCGAACTCCCGCCTTTATCTTT‑AAAAAGTAAACCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                                                 <  1:577780/140‑1 (MQ=255)                                                   attaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAggg                                                                                         >  2:190164/1‑151 (MQ=255)                                                      aaCCCCATGCCGCGGGGCGATTTTTACCAGACGCTGACGCCCCGCCTCCCCTCTGCGAAATCTCGCCTTTATTTTA‑AAAAATGGAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAggg                                                                                         <  1:193026/148‑1 (MQ=255)                                                          ccAGGACGCGGGGCGATTTTTACCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                    <  2:213347/149‑1 (MQ=255)                                                                  ccGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCGGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATCCCGCGA‑AAAAATCCTACCGCCGCCGACAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                           <  2:658154/150‑1 (MQ=255)                                                                  ccGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                          >  2:193135/1‑151 (MQ=255)                                                                    gggCGTTTTTTATCAGCCGCTGACGCCCCCCCATCCCTCTGCGAAATCTCGCCGTTATTTTT‑ATAAAGTAAATCCCCCGA‑AAAAATCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCca                                                                         <  2:829063/150‑1 (MQ=255)                                                                    gggAGATTTTTATCAGACGCGGACGCCCCCCCACCCCTCTGCGAAATCTCGCCGTTATCTTT‑AAAAAGTGAATCCCGCGA‑AAAAACCCTCCCGCCGCAGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                             <  1:139754/146‑1 (MQ=255)                                                                     ggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                        >  1:133693/1‑150 (MQ=255)                                                                       ctaTGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                                                      >  2:189634/3‑150 (MQ=255)                                                                        gATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                     >  1:237970/1‑150 (MQ=255)                                                                               accAGCCGCTGACGCCCCGCCACCCCTCTGCGGATCCCCGCCTTTATCTTA‑AAAAAGTAAATCCCCCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                             <  2:1047532/149‑1 (MQ=255)                                                                                   gcctctgccccccccccACCCCTCTGCGGATTCTCGCCGTAATCTTA‑ATAAAGTAAATCCCGCGA‑AAAATCCCTCCCGCCGCCGATAATTCAATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGttt                                                         <  2:715776/140‑1 (MQ=255)                                                                                      gcggacgccCCGCCACCCCTCTGCGGAATCTCGCCGTTATTTTT‑ATAAAGTGAATCCCGCGA‑AAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                          <  2:917219/144‑1 (MQ=255)                                                                                      gCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                      >  2:582020/1‑151 (MQ=255)                                                                                        tGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                                                     >  1:1037508/1‑150 (MQ=255)                                                                                        tGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                                                     >  2:783948/1‑150 (MQ=255)                                                                                           cgcccccccACCCCTCGGCGAAATCTCGCCGTTATCTTT‑AAAAAGTAAATCCCGCGA‑AAAAACCCTCCCCCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAt                                                  <  1:1114695/147‑1 (MQ=255)                                                                                              cccGCCACCCCTCTGCGGAATCTCGCCTTTATTTTT‑AAAAAGTAAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAtt                                                 <  2:133727/148‑1 (MQ=255)                                                                                                 cccACCCCTCTGGGAAATCTCCCCTTTTTCTTA‑AAAAAGTGAACCCCGCCA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                           <  1:998575/150‑1 (MQ=255)                                                                                                  cctccccTCGGGGAAATCTCCCCTTTATCTTT‑AAAAAGTGAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                           <  1:489772/146‑1 (MQ=255)                                                                                                  ccATCCCTCTGGGAATTCTCGCCGTTATTTTT‑AAAAAGTGAATACCCCGA‑AAAAACCCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                             <  2:91054/148‑1 (MQ=255)                                                                                                  caaccccTCTGCGTAACCGCGCCGTTTTTTTT‑AAAAAGAAAATACCGCCA‑AAAAACCCCCCCGCCGCAGAAAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                           <  2:751227/146‑1 (MQ=255)                                                                                                   cATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                >  1:399424/1‑112 (MQ=255)                                                                                                      cccccggCGAAATCTCCCCTTTATTTTT‑AAAAAGGAAATCCCGCGA‑AAAAATCCTACCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                        <  1:85695/143‑1 (MQ=255)                                                                                                      cccTCGGCGGAATCTCGCCGTTATTTTT‑AAAAAGTGAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                <  2:399463/109‑1 (MQ=255)                                                                                                      cccTCGGCGAAATCTGGCCGTAATTTTT‑AAAAAGTAAATCCCGCGA‑AAAAACCCTACCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                      <  1:991939/151‑1 (MQ=255)                                                                                                       ccTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                     >  2:452195/1‑151 (MQ=255)                                                                                                        ctctGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                     >  2:462165/1‑150 (MQ=255)                                                                                                            tccgAATCCCGCCTTTATCTTT‑AAAAAGTGAATCCCGCGA‑AAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                      <  1:987751/142‑1 (MQ=255)                                                                                                              ggAATCTCGCCGTTATCTTT‑AAAAAGTAAATCCCGCGA‑AAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                              <  2:856851/151‑1 (MQ=255)                                                                                                              gaactccccctgTTATCTTT‑ATAAAGGATATACC‑CGCAAAAAACCCCGCCCACGCCGAAAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                <  1:567582/138‑1 (MQ=255)                                                                                                                aaTCTCCGCATTATTTTTAATAAGGTGACTACCCCGA‑AAAAACCCTCCCCCCGCCAATAATCCCTTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt                             <  1:140518/151‑1 (MQ=255)                                                                                                                  tctcGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGcc                          >  1:725344/1‑151 (MQ=255)                                                                                                                     cGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                      >  2:987844/1‑136 (MQ=255)                                                                                                                     cGCCGTTATCTTT‑ATAAAGTAAATCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                           <  2:57692/147‑1 (MQ=255)                                                                                                                        cGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                     >  1:1077401/1‑150 (MQ=255)                                                                                                                        cGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGt                    >  2:118825/1‑151 (MQ=255)                                                                                                                              cTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaa               >  2:486325/1‑150 (MQ=255)                                                                                                                              cTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaa               >  1:339865/1‑150 (MQ=255)                                                                                                                                 t‑ataAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaaa            >  1:435416/1‑150 (MQ=255)                                                                                                                                   ataAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc          >  2:263284/1‑151 (MQ=255)                                                                                                                                    taAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc         >  1:452990/1‑151 (MQ=255)                                                                                                                                       aGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAg       >  2:890674/1‑150 (MQ=255)                                                                                                                                           aaTCCCGCGA‑AAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc    <  1:134439/149‑1 (MQ=255)                                                                                                                                           aaTCCCGCAA‑AAAAACCCGCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc    <  1:1108873/149‑1 (MQ=255)                                                                                                                                           aaTCCCCCGA‑AAAAACCCTACCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc    <  1:3603/149‑1 (MQ=255)                                                                                                                                           aaTACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc  >  2:780316/1‑151 (MQ=255)                                                                                                                                           aaTACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGATATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc  >  2:773456/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTT‑ATAAAGTGAATACCGCGC‑AAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGC  >  NC_000913/3105343‑3105630

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