Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,105,483	0	A	C	44.3%	119.3 / 11.5	70	T351P (ACC→CCC)	mltC	membrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (34/5);  new base C (0/31);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 5.36e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCC  >  NC_000913/3105335‑3105628                                                                                                                                                     |                                                                                                                                                  aggcgCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATAcc                                                                                                                                                 >  2:601125/2‑151 (MQ=255)           agCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAAt                                                                                                                                       >  1:346510/1‑151 (MQ=255)              gTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCt                                                                                                                                    >  1:142727/1‑151 (MQ=255)               tgctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCt                                                                                                                                    >  2:620818/1‑150 (MQ=255)                  tgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccg                                                                                                                                >  2:1060617/1‑151 (MQ=255)                     gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccg                                                                                                                             >  1:222826/1‑151 (MQ=255)                     gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCAGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccg                                                                                                                             >  1:693541/1‑151 (MQ=255)                       gTCTTTTCGAATGATAAGATTCAGGCTGCAAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgcc                                                                                                                           <  1:574563/151‑1 (MQ=255)                              cGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAAtt                                                                                                                    >  2:761909/1‑151 (MQ=255)                              cGAATGAAAAGATTAGGGCTGCGAATATTATTACCCCCAGGACGCGGGACGATGTTTATCAAACGTTAACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATaa                                                                                                                      <  1:364937/149‑1 (MQ=255)                                    aTAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAtt                                                                                                               >  2:29739/1‑150 (MQ=255)                                        gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                          >  1:56354/1‑151 (MQ=255)                                          ttCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGcc                                                                                                        >  2:17478/1‑151 (MQ=255)                                             aGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCt                                                                                                     >  1:893849/1‑151 (MQ=255)                                               gCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                    >  1:780963/1‑150 (MQ=255)                                                 tGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACg                                                                                                 >  1:915948/1‑151 (MQ=255)                                                 tGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACg                                                                                                 >  2:965928/1‑151 (MQ=255)                                                     aaTATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGc                                                                                             >  2:776217/1‑151 (MQ=255)                                                        attattAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                           >  1:485933/1‑150 (MQ=255)                                                         ttataacccccATGCCCCCGGGCGATGTTTATCAGACGCTGACGCCCCCCCACCCCTCTGCGGAATCTCCCCTTTATTTTTATAAAGTAAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCt                                                                                            <  1:359464/139‑1 (MQ=255)                                                            ttacccccATGCCCCGGGGCGTTTTTTATCAGACGCTGACGCCCCCCCACCCCTCGGCGAATCCTCGCCTTTATTTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAggg                                                                                       <  2:193862/144‑1 (MQ=255)                                                              aaCACCATGCCCCGGGGGTATTTTAACCAGACCTAGACCCCCCGCCCCCCCTCTGCGAAATCCCGCCTTTTTTTTTAAAAAGTAAATCCCGCGAAAAACCCCTCCCGCCGCCAATAATCCCATTAACCGCCCCTGACGATGCTc                                                                                           <  2:7349/144‑1 (MQ=255)                                                                 cccATGCCGCGGGGCGTTGTTTATCAGCCGCTGACGCCCCCCCACCCCTCGGCGAATTCTCGCCGTTTTTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                   <  2:429972/148‑1 (MQ=255)                                                                  ccATGACGCGGGGCGTTTTTTATCAGCCGCTGACGCCCCCCCACCCCCCTGCGGATTCTCGCCGTTATTTTTAAAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                 <  1:71810/150‑1 (MQ=255)                                                                    aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                              >  2:631060/1‑151 (MQ=255)                                                                          cggggggttGTTTATCGGACGCTGACGCCCCCCCCCCCCTCTGCGAAATCCCGCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAAATCCTCCCCCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                               <  2:1010098/136‑1 (MQ=255)                                                                          cggggCGATTTTTATCAGCCGCTAACGCCCCCCCACCCCTCTGCGAAATCTCGCCGTTTTTTTTATAAAGGGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                                          <  2:150222/147‑1 (MQ=255)                                                                          ccGGGCGATGTTTATCAGCCGCTGACGACCCGCCACCCCTCTGCGGAATCTCGCCGTTATCTTTAAAAAGTGAATACCGCGAAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                           <  1:537640/148‑1 (MQ=255)                                                                           cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                        >  1:618969/1‑150 (MQ=255)                                                                             ggCGTTTTTTATCAGCCGCTGCCGCCCCCCCACCCCCCTGCGAAATCTCCCCTTTTTCTTTAAAAAGGAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcaca                                                                     <  2:361063/151‑1 (MQ=255)                                                                                gATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                   >  1:114603/1‑150 (MQ=255)                                                                                gATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                   >  1:369319/1‑150 (MQ=255)                                                                                  tGTTTATCAGCCCTTGACGCCCCGCCACCCCCCGGCGAAATCTCGCCTTTATTTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                  <  2:368219/149‑1 (MQ=255)                                                                                   gTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                >  1:76129/1‑150 (MQ=255)                                                                                    tttATCAGCCGCTGCCGCCCCCCCTCCCCTCTGCGAATTCCCGCCGTTATTTTTAAAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                <  2:469033/149‑1 (MQ=255)                                                                                    ttaaTCAGAGCCTGACCCCCCCCCACCCCGCTGCGGAATCTGGCGGTTTTCTTAATAAGGTGAACCCCGCGAAAAAATCCTCCCCCGGCCGATATTTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                                                 <  2:436940/145‑1 (MQ=255)                                                                                             cgttgccgccccGCCACCCCCCGGCGAATCCCCGCCGTTATTTTTAAAAAGTAAATCCCGCGAAAAAACCCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                      <  2:612882/141‑1 (MQ=255)                                                                                             cGCTGACGCCCCGCCACCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                           <  1:55403/145‑1 (MQ=255)                                                                                              gCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                    >  2:748732/1‑151 (MQ=255)                                                                                              gCTGACACCCCCCCACCCCTCTGCGAAATCCCGCCGTTATTTTTATAAAGGGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                        <  2:899300/147‑1 (MQ=255)                                                                                                   cgccccGCCACCCCTCGGCGGAATCCCGCCGTTATTTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                           <  1:620795/136‑1 (MQ=255)                                                                                                    gacCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                              >  2:930901/1‑151 (MQ=255)                                                                                                      ccccccACCCCTCTGGGAAATCTCGCCTTTATTTTTAAAAAGGAAATCCCGCGAAAAAATCCTCCCGCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                              <  2:908176/149‑1 (MQ=255)                                                                                                      ccccccACCCCTCTGCGGAATCCCGCCTTTATCTTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                          <  1:1030098/137‑1 (MQ=255)                                                                                                      ccccccACCCCTCTGCGAATTCTCGCCTTTTTTTTTATAAATTAAATCCCGCGAAAAAATCCTCCCGCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                              <  1:962641/149‑1 (MQ=255)                                                                                                      cccGCCACCCCTCGGCGGAATCTCGCCGTTATCTTTAAAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                            <  1:262965/151‑1 (MQ=255)                                                                                                        cGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGa                                           >  2:1074377/1‑150 (MQ=255)                                                                                                        cGCCATCCCTCGGCGAATTCTCGCCGTTATTTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGaa                                          <  1:889981/151‑1 (MQ=255)                                                                                                           cATCCCTCTGGGGACTCTCGCCTTTTTCTTTATAAAGTGAATACCGCGCAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                       <  1:1060575/151‑1 (MQ=255)                                                                                                            aTCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                      >  2:11740/1‑151 (MQ=255)                                                                                                            aTCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                      >  2:240747/1‑151 (MQ=255)                                                                                                              cccTCTGCGAATTCCCGCCTTTATTTTTAAAAAGTAAAACCCGCCAAAAACCCCTCCCGCCGCCAATATTCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt                                     <  2:1057026/150‑1 (MQ=255)                                                                                                                    gCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                               >  2:188654/1‑150 (MQ=255)                                                                                                                      gaaaccccccctttATTTTTAAAAAGTGAATCCCGCGAAAAACTCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                            <  2:915982/139‑1 (MQ=255)                                                                                                                        aaTCTCGCCTTTATTTTTAAAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                          <  1:744817/151‑1 (MQ=255)                                                                                                                        aaTCTCGCCGTTATTTTTAAAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt                           <  1:987008/150‑1 (MQ=255)                                                                                                                        aaTCCCCCCTTTATTTTTAAAAAGTAAATCCCGCAAAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                               <  2:346517/146‑1 (MQ=255)                                                                                                                         aTCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                         >  2:530275/1‑151 (MQ=255)                                                                                                                          tctcGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGcc                        >  2:354443/1‑151 (MQ=255)                                                                                                                            ccccctttATATTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGAAAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                          <  2:533733/141‑1 (MQ=255)                                                                                                                             cGCCTTTATTTTTAAAAAGGGAACCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt                           <  2:142739/145‑1 (MQ=255)                                                                                                                                cGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGt                  >  2:259961/1‑151 (MQ=255)                                                                                                                                 gTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                   >  1:785722/1‑149 (MQ=255)                                                                                                                                    aTCTTTATAAAGTGAATCCCGCGCAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTaa                <  1:17479/149‑1 (MQ=255)                                                                                                                                     tCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACa              >  1:806386/1‑150 (MQ=255)                                                                                                                                       tttAAAAAGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc               <  1:654395/147‑1 (MQ=255)                                                                                                                                          aaaaaggaaaTCCCGCGCAAAACCCCCCCCCCCGCAGTTAATTCCATTACCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc        <  2:135797/143‑1 (MQ=255)                                                                                                                                            aaaGTAAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCa      <  2:386015/151‑1 (MQ=255)                                                                                                                                             aaGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc       >  2:44744/1‑149 (MQ=255)                                                                                                                                                  aaTCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc  <  1:631037/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGGCGCAGGCAGCGTGCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCC  >  NC_000913/3105335‑3105628

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