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

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCC  >  NC_000913/3105350‑3105646                                                                                                                                                        |                                                                                                                                                    gctgcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAc                                                                                                                                                    <  1:237508/151‑1 (MQ=255)     cGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgcc                                                                                                                                                >  2:169991/1‑151 (MQ=255)      gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgccg                                                                                                                                               >  1:822419/1‑151 (MQ=255)           ttttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGAt                                                                                                                                          >  2:822824/1‑151 (MQ=255)               cGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGCCGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCTTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATaa                                                                                                                                        <  2:867774/149‑1 (MQ=255)                  aggATAAAATTCAGGCTGTGATTATAATTAACACCATGGCGCCGGGCGATGTTTATCAGACGCTAAAGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTGCCGCCGGCGATAATTcc                                                                                                                                    <  1:297191/148‑1 (MQ=255)                       aaGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAc                                                                                                                              >  1:736908/1‑151 (MQ=255)                        aGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAcc                                                                                                                             <  2:1024352/151‑1 (MQ=255)                         gATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCg                                                                                                                            <  1:857347/151‑1 (MQ=255)                              aGGCTGCAAATATTATAACCCCCATGCCGCGGGGCGATTTTTATCAGCCGCTGACGCCCCGCCACCCC‑TCGGCGAAATCTCGCCGTTATCTTTATAAAGTG‑AATCCCGCAAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGcccc                                                                                                                        <  1:713099/150‑1 (MQ=255)                               ggCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                                      >  1:630976/1‑151 (MQ=255)                                          ttattaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                             >  2:849490/1‑149 (MQ=255)                                          ttattaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                             >  1:759681/1‑149 (MQ=255)                                              taACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGc                                                                                                       >  2:460063/1‑151 (MQ=255)                                                 cacCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                     >  1:75671/1‑150 (MQ=255)                                                  acCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaa                                                                                                    >  1:905952/1‑150 (MQ=255)                                                        aCGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGtt                                                                                              >  2:354492/1‑150 (MQ=255)                                                          gCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                                           >  2:866046/1‑151 (MQ=255)                                                           cggggCGTTGTTTACAAGCCGCTGCGGCCCCCCCACCCC‑CTGGCGGAACCCCCCCTTTATTTTTAAAAAGTA‑ATTCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                                   <  1:521189/140‑1 (MQ=255)                                                           ccGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                                          >  1:613120/1‑151 (MQ=255)                                                            cGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                                                                          >  1:90364/1‑150 (MQ=255)                                                             ggcggttttttATCGGTCGCTGGCGCCCCGCCACCCC‑TCTGCGGAATCCGGCTGTTATCTTTAAAAGGAG‑AATCCCGCGAAAAAATCTTCCCCCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                                             <  2:174689/138‑1 (MQ=255)                                                              ggCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                        >  2:387225/1‑150 (MQ=255)                                                              ggCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                        >  1:303025/1‑150 (MQ=255)                                                                   tttttACCAGCCGCTGCCGCCCCCCCACCCC‑TCTGCGGAATCTCGCCGTTATTTTTAAAAAGTG‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                                    <  2:635006/147‑1 (MQ=255)                                                                     tttATCAGCCGCTGACCCCCCCCCCCCCC‑TCTGCGAAATCTCGCCGTTATCTTTAAAAAGTG‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                                    <  1:378736/147‑1 (MQ=255)                                                                     tttATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAt                                                                                >  2:508737/1‑151 (MQ=255)                                                                        accAGACGCTACCGCCCCCCCACCCC‑TCTGGGAATTCTCGCCTTTATCTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                              <  1:108580/148‑1 (MQ=255)                                                                        aTCAGCCGCTGACGCCCCGCCACCCC‑TCGGGGAAATCTCGCCGTTATCTTTAAAAAGTG‑AATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                              <  2:648572/150‑1 (MQ=255)                                                                          cAGCCGCTGACCACCCGCCACCCC‑TCGGCGGAATCCCGCCGTTATTTTTATAAAGTA‑AATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                           <  2:539764/151‑1 (MQ=255)                                                                             aCGCTGACGCCCCCCCACCCC‑TCTGCGAATTCCCCCCTTTATTTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                            <  1:502727/147‑1 (MQ=255)                                                                               gCTGACGCCCCGCCACCCC‑TCTGGGAAATCTCCCCTTTATTTTTATAAATTA‑AATCCCGCGCAAAAACCCTCCCCCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                                                                        <  2:514474/149‑1 (MQ=255)                                                                               gCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                      >  1:100395/1‑151 (MQ=255)                                                                                   acgacCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAt                                                                  >  2:425184/1‑151 (MQ=255)                                                                                    cgacCGGCCACCCC‑CCTCCGAAATCCCGCCTTTATTTTTAAAAGGTG‑AATCCCCCCCAAAAACCCTCCCCCCTCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaaa                                                                   <  2:858151/149‑1 (MQ=255)                                                                                     cgcccGCCCTCCCTTCGGCGAAACCTCGCCTTTATTTTTAAAAAGAA‑AATCCCGGGAAAAAACCCTCCGGCCGCTGATATTTCAATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                              <  1:1055295/136‑1 (MQ=255)                                                                                       ccccccACCCC‑TCTGGGAATTCCCGCCGTTATCTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                    <  2:615411/145‑1 (MQ=255)                                                                                       ccccccACCCC‑TCGGGGAATTCCCGCCTTTATTTTAAAAAAGTA‑AACCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaaa                                                                   <  1:987012/146‑1 (MQ=255)                                                                                       cccGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                              >  1:524164/1‑151 (MQ=255)                                                                                       cccGCCACCCC‑TCGGGGAAATCTCGCCTTTATTTTTAAAAAGTG‑AATCCCGCGAAAAATTCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                    <  2:1101297/145‑1 (MQ=255)                                                                                               cccctcgcgggAACCCCGCCTTTATTTTAAAAAAGTA‑ATTGCCCCGAAAAACCCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                                        <  1:15958/141‑1 (MQ=255)                                                                                               ccc‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                      >  2:808768/1‑151 (MQ=255)                                                                                                 cccctGCGAAATCCCCCCTTTATTTTTAAAAAGTA‑AATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                     <  1:971512/148‑1 (MQ=255)                                                                                                      gCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                 >  2:945173/1‑150 (MQ=255)                                                                                                        ggAATCTCGCCGTTATTTTTATAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                                <  2:673435/149‑1 (MQ=255)                                                                                                         gAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt                                             >  1:224976/1‑151 (MQ=255)                                                                                                               ccccGTTATTTTTATAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcaca                                       <  1:888162/149‑1 (MQ=255)                                                                                                                gCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                        >  2:541603/1‑149 (MQ=255)                                                                                                                 ccGTTATCTTTAAAAAGTG‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCca                                         <  2:13988/147‑1 (MQ=255)                                                                                                                  ctttttCTTTATAAAGTG‑AATCCCGCGCAAAAACCCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                                     <  2:500566/145‑1 (MQ=255)                                                                                                                     tATCTTTATAAAGTG‑CATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                 >  2:316545/1‑151 (MQ=255)                                                                                                                      aTCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                 >  2:754616/1‑150 (MQ=255)                                                                                                                          ttATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                             >  2:50480/1‑150 (MQ=255)                                                                                                                           tataAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaaa                            >  2:79225/1‑150 (MQ=255)                                                                                                                            ataAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                          >  2:949578/1‑151 (MQ=255)                                                                                                                             taAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                         >  2:177362/1‑151 (MQ=255)                                                                                                                             aaaaaggaaaTCCCGCCAAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaaaa                           <  2:759718/142‑1 (MQ=255)                                                                                                                               aaGTGAAATCCCGCCAAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                         <  1:508712/150‑1 (MQ=255)                                                                                                                                aggggaacccggcgaAAAAACCCCTCCCCGGCCGATAATTCTGTTATCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTTCAATTATTTGCCACAGGTAACaaaaaa                           <  1:1136114/136‑1 (MQ=255)                                                                                                                                 gTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTc                     >  2:412249/1‑151 (MQ=255)                                                                                                                                         ccccgcAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                 <  2:636197/145‑1 (MQ=255)                                                                                                                                         ccGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa               <  2:359774/150‑1 (MQ=255)                                                                                                                                             gcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTg           >  2:114091/1‑150 (MQ=255)                                                                                                                                               aaaatTCCTGCTGCCGCAGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa        <  2:361488/151‑1 (MQ=255)                                                                                                                                               aaaaaCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt         <  2:75674/150‑1 (MQ=255)                                                                                                                                               aaaaaCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa        <  1:734475/151‑1 (MQ=255)                                                                                                                                               aaaaaCCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGa          <  2:731025/149‑1 (MQ=255)                                                                                                                                                     ccTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATc   <  2:416701/150‑1 (MQ=255)                                                                                                                                                     ccTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATcc  >  1:396511/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCC‑TCTGCGGAATCTCGCCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCC  >  NC_000913/3105350‑3105646

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