Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,105,499	0	A	C	48.3%	86.7 / 10.1	60	Y356S (TAC→TCC)	mltC	membrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/29);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  NC_000913/3105350‑3105647                                                                                                                                                        |                                                                                                                                                     gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTa                                                                                                                                                      <  2:439314/150‑1 (MQ=255)    gcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTAccgc                                                                                                                                                  <  1:541362/151‑1 (MQ=255)       aGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTAccgccgc                                                                                                                                               >  2:390470/1‑151 (MQ=255)        gTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTAccgccgcc                                                                                                                                              >  1:584440/1‑151 (MQ=255)        gTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTAccgccgcc                                                                                                                                              >  1:666274/1‑151 (MQ=255)         tCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTAccgccgccg                                                                                                                                             >  1:182099/1‑151 (MQ=255)              tCGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATaa                                                                                                                                         <  2:2207/150‑1 (MQ=255)               cGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAAt                                                                                                                                        <  2:920200/150‑1 (MQ=255)                gAATGAAAAGATTCAGGCTGCCAATATTATTACCACCATGACGCCGGGCGATGTTTATAAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAAtt                                                                                                                                       <  1:633782/150‑1 (MQ=255)                 aaTGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTcc                                                                                                                                     >  2:125484/1‑151 (MQ=255)                       aaGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTaa                                                                                                                                <  2:770341/150‑1 (MQ=255)                          aTTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGc                                                                                                                            >  2:272546/1‑151 (MQ=255)                                 cTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGAc                                                                                                                     >  2:506548/1‑151 (MQ=255)                                    ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                                                                  >  2:1151198/1‑151 (MQ=255)                                            attaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAggg                                                                                                          >  1:315297/1‑151 (MQ=255)                                             taacccccAGGACGCGGGGCAATTTTTATCAGACGCTGCCGCCCCCCCACCCCTCTGCGGAATCTCGCCT‑TTATCTTTATAAAGGGAATCCCGCGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAg                                                                                                            <  2:1024608/142‑1 (MQ=255)                                                 ccccaggccGCCGGGCGATTTTTATCAGCCGCTGCCGCCCCCCCACCCCTCGGCGGATCCCCGCCG‑TTATCTTTAAAAAGTAAATCCCGCGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                      <  2:795276/142‑1 (MQ=255)                                                   ccATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                                    >  1:862237/1‑150 (MQ=255)                                                     aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                 >  2:577186/1‑151 (MQ=255)                                                      gcccccggcggtttGTTTATAAGCCGTAGACCCCCCCCCTCCCCTCGGCGGATCCCCGCCT‑TTGTCTTTAAAAAGTAAACCCCGCGAAAAATG‑CCACCCCCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                 <  1:365997/137‑1 (MQ=255)                                                             ggggggtGTTTAACAGACGCTGACGCCCCGCCGCCCCTCTGTAAACTCTCGCCG‑TTATCTTTAAAAGGTAAATCCCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                                           <  2:1022503/143‑1 (MQ=255)                                                             gggCGATTTTTTTCGGCCGCTGCCCCCCCCCCACCCCCCGGGGAATTCCCGCCG‑TTTTTTTAAAAAAGTAAATCCCGCCAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                                                             <  2:1080313/147‑1 (MQ=255)                                                             gggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                         >  2:1075565/1‑151 (MQ=255)                                                                   tGTTTATCAGCCGCTAACGCCCCGCCACCCCTCGGGGAATTCTCCCCT‑TTATTTTTAAAAAGTAAATCCCGGGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                                                                       <  2:869063/147‑1 (MQ=255)                                                                    ttttggcAGAGGCAGCCGCCCCGCCA‑CCCGCCGGGGAATGTCCCCCTTTATTTTAAAAAAGTAAATCCCGCGAAAAAAC‑CCACCCGCCGCCGTTATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                                      <  1:845903/141‑1 (MQ=255)                                                                             aCGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                         >  1:666845/1‑151 (MQ=255)                                                                               gCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                       >  2:56419/1‑151 (MQ=255)                                                                                 tgccgcccccccACCCCTCTGCGGAATCTCGCCG‑TTATTTTTAAAAATTGAATCCCGCGAAAAAAT‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGttt                                                                          <  2:1076897/140‑1 (MQ=255)                                                                                  gacgacCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                     >  2:765147/1‑150 (MQ=255)                                                                                       cccccccccccTCTGCGGAATCTCGCCG‑TTTTCTTTAAAAAGTGAATCCCGCGCAAAAAT‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                 <  2:1053110/141‑1 (MQ=255)                                                                                       ccccccaccccctggCGGACTCTCGCCT‑TTATCTTTAAAAAGTGAACCCCCCGAAAAAAC‑CCTCCCGCCCCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                         <  1:342697/127‑1 (MQ=255)                                                                                        ccGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                               >  2:1066349/1‑150 (MQ=255)                                                                                           ccACCCCTCTGCGAATTCTCGCCT‑TTATTTTTAAAAAGTAAATCCCCCAAAAAAAC‑CCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                            <  1:1173785/150‑1 (MQ=255)                                                                                           ccACCCCTCTGCGAAATCCCCCCG‑TTATTTTTAAAAAGTGAATCCCGCGCAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                            <  2:602071/150‑1 (MQ=255)                                                                                              ccccTCGGCGAATTCTCGCCG‑TTATTTTTAAAAAGTAAACCCCCCGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                                          <  2:475749/148‑1 (MQ=255)                                                                                               cccccggCGAAATCCCCCCT‑TTTTCTTTAAAAAGTAAATCCCGCGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt                                                        <  2:279567/144‑1 (MQ=255)                                                                                               cccccgccgaaaTCTCCCCG‑TTATTTTTAAAAAGTGAATCCCGCGAAAAAAT‑CCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                                         <  2:1117741/139‑1 (MQ=255)                                                                                               cccTCGGCGAAATCTCGCCG‑TTATTTTTATAAAGTAAATCCCGCGAAAAAAT‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                                         <  2:577758/149‑1 (MQ=255)                                                                                                ccTCTGCTGAATCCCGCCG‑TAATCTTTAAAAAGAAAATCCCGCGAAAAAAC‑CTTCCCCCCGCCGATAATTCCTTTAGCCGCCCCTGACGATGCTCAGGGCCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                      <  2:459308/151‑1 (MQ=255)                                                                                                         atccccGCCG‑TTTTCTTTAAAAAAAAAACCCCCCGAAAAACT‑CCTCCCCCCGCCAATATTTCCTTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                             <  1:1083540/146‑1 (MQ=255)                                                                                                         aaTCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                             >  1:921179/1‑151 (MQ=255)                                                                                                          aTCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                                            >  2:782681/1‑151 (MQ=255)                                                                                                          aTCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                                            >  2:1075410/1‑151 (MQ=255)                                                                                                            ctcCCCG‑TTATTTTAAAAAAGTAAATACCGCAAAAAAAC‑CCTCCCGCCGCAAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTg                                             <  1:1075511/148‑1 (MQ=255)                                                                                                               gCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                         >  2:909106/1‑149 (MQ=255)                                                                                                                ccG‑TTATTTTTATAAAGGGAACCCCGCGCAAAAAT‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                                      <  1:485671/151‑1 (MQ=255)                                                                                                                ccG‑TTATTTTTAAAAAGGAAATCCCGCCAAAAAAC‑CCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAg                                       <  2:199030/150‑1 (MQ=255)                                                                                                                    ttttttttATAAAGTAAACCCCGCCAAAAAAC‑CCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTa                                    <  1:898639/145‑1 (MQ=255)                                                                                                                         tttAAAAAGTAAACCCCGCAAAAAAAT‑CCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaa                                <  1:752394/149‑1 (MQ=255)                                                                                                                              aaaggaaacaCCGCCAAAAATT‑CCCCCCGCCGCCATTAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCGCATCACAATTTCGTTTTGCAAATTGGGAATGTTTACAATTATTTGCCACAGTTAACaaaa                              <  2:584657/137‑1 (MQ=255)                                                                                                                              aaaGTGAACCCCGCAAAAAAAC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                           <  1:410735/149‑1 (MQ=255)                                                                                                                               aaagtaaagcccgccAAAAATCCCTCCCCCCCCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCa                         <  1:339166/137‑1 (MQ=255)                                                                                                                                       cccGCGAAAAAAC‑CCTCCCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                  <  2:711654/148‑1 (MQ=255)                                                                                                                                       aCCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa                 >  1:150177/1‑150 (MQ=255)                                                                                                                                         ccccaaaaaaccccTCCCCCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                  <  1:364992/136‑1 (MQ=255)                                                                                                                                            gcAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTg            >  1:1157670/1‑150 (MQ=255)                                                                                                                                              aaaaaC‑CCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt          <  1:316506/150‑1 (MQ=255)                                                                                                                                                  aacccTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGaa      <  2:854662/148‑1 (MQ=255)                                                                                                                                                     ccTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc   >  1:382603/1‑151 (MQ=255)                                                                                                                                                      cTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATccc  >  2:557026/1‑151 (MQ=255)                                                                                                                                                       tACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATccc  >  2:85589/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                     GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCG‑TTATCTTTATAAAGTGAATACCGCGCAAAAAT‑CCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑