Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,324	0	T	G	41.5%	121.1 / 19.4	65	G33G (GGT→GGG)	ynbB	putative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (5/33);  new base G (27/0);  total (32/33)
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

GATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAT  >  NC_000913/1478181‑1478471                                                                                                                                                 |                                                                                                                                                    gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAg                                                                                                                                                  <  1:280310/147‑1 (MQ=255)            gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGGGAAAAGGGGGGGAGGGTACATTACGTAtt                                                                                                                                   >  1:906092/1‑151 (MQ=255)            gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGGGAAAAGGGGGGGAGGGTACATTACGTAtt                                                                                                                                   >  1:22484/1‑151 (MQ=255)            gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTAtt                                                                                                                                   <  1:851990/151‑1 (MQ=255)            gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCCGGGGGAAAAGGGGGGGGGGGTAAATTTCGTAtt                                                                                                                                   >  1:6423/1‑151 (MQ=255)             tCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGGAAAAGGGGGGGGGTTACCTTTCGGAtt                                                                                                                                   >  1:120419/1‑150 (MQ=255)               ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATCGGGGGGAAAAGGGGGGGGGGGTTCATTTTGTAttt                                                                                                                                  >  2:992850/1‑148 (MQ=255)                  gtgtGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTAc                                                                                                                              <  2:494536/150‑1 (MQ=255)                      gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGGAAAAGGGGGGGGGGTACATTACGTATTCGTACCttg                                                                                                                          >  1:27952/1‑148 (MQ=255)                       cttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCtggtg                                                                                                                        >  1:673115/1‑151 (MQ=255)                           ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAAGGGGGGAGGGTATATTTCGTATTTGTACCTGggggggg                                                                                                                    >  2:17438/1‑150 (MQ=255)                             aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATCTGGGAGCAAAGGGGGGGAGGTGAAGTTACGGATTCGTGCCTTGTtggg                                                                                                                     >  2:910220/1‑148 (MQ=255)                             aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAGGGGGGGGGGGTTCATTACGTATTCGTACCTGgggggg                                                                                                                     >  1:679371/1‑148 (MQ=255)                              tGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTACGGGGAGAAAAGGGGGGGAGGGTACATTTCGTTTTTGTTCCTGGCGGGGcgt                                                                                                                  >  2:730389/1‑148 (MQ=255)                               gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                                                                                <  1:52750/151‑1 (MQ=255)                                    cTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGGGAAAAGGGGGGGAGGTTACATTACGTATTCGTACCTGGTGGGGGCTCGTCCtt                                                                                                           >  1:783409/1‑151 (MQ=255)                                      gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAttt                                                                                                          <  1:718419/150‑1 (MQ=255)                                         atatAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGttt                                                                                                      <  2:43187/151‑1 (MQ=255)                                          tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGTGGGGGATCGTCATTTGtttt                                                                                                     >  1:986594/1‑151 (MQ=255)                                           ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAAGGGGGGGGGGTATATTACGTATTTGTTTCCGGTGGGGGCTCGTCGTTTGttttt                                                                                                    >  2:13576/1‑151 (MQ=255)                                           ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGtttt                                                                                                     <  1:386316/150‑1 (MQ=255)                                             aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTc                                                                                                   <  2:135330/150‑1 (MQ=255)                                             aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCa                                                                                                  <  1:420247/151‑1 (MQ=255)                                             aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCa                                                                                                  <  2:67886/151‑1 (MQ=255)                                               tGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGCATCGTCCTtttttttt                                                                                                    >  1:363830/1‑147 (MQ=255)                                                      aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGAtt                                                                                         <  1:424718/151‑1 (MQ=255)                                                           ctctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATCTGAGGGAAAAGGGGGGGAGGGTACAATTAGTATTTGTACCCTGTGGGGGCTCGTCCTTTGTTTTTTTATTGTGAgtttcgg                                                                                    >  1:937686/1‑147 (MQ=255)                                                             ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                                                                      <  2:511543/147‑1 (MQ=255)                                                               ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTc                                                                                  <  2:666104/149‑1 (MQ=255)                                                               ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCa                                                                                <  2:354065/151‑1 (MQ=255)                                                                        tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGa                                                                         <  1:598936/149‑1 (MQ=255)                                                                         ttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTAACGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTATCTGGGGGGGGATCGTCATTTTTTTTTTAATGGGGATTTCTGGGCCACCCTGGAtg                                                                       >  1:813148/1‑150 (MQ=255)                                                                          tCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGAAAAGGGGGGGAGGGTTAATTACGTATTTGGACCTGGGGGGGGATTGGCATTTTTTTTTTTATTGTGATTTTTGGGCCACGCTGGGtggc                                                                     >  2:617657/1‑149 (MQ=255)                                                                           cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGAGGAAAGGGGGGGGGGGTTAATTTCCTATTTGTACCTGGTGGGGGGTCGTCATTTGTTTTTTAATTGGGAGTTTAGGGCCCCGCCGGGttacg                                                                    >  2:941674/1‑151 (MQ=255)                                                                           cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacg                                                                    <  2:490602/151‑1 (MQ=255)                                                                              gttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttg                                                                 <  1:606870/151‑1 (MQ=255)                                                                               ttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTAACGGGGGGAAAAAGGGGGGGGGGTTTATTTCGTATTTGGACCTGGTGGGGGATTGTCATTTGTTTTTTTATGGGGGTTTTTGGTCCCCGCTGGGtgacgtt                                                                  >  2:13495/1‑149 (MQ=255)                                                                                tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGGGAAAAGGGGGGGAGGGTACATTTTGTTTTTCTACCTGGGGGGGGATCGTAATTTTTTTTTTAAAGGGGGTTTTATGTGCACGCTGGAtggag                                                                    >  2:580729/1‑143 (MQ=255)                                                                                 gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgac                                                               <  2:519679/150‑1 (MQ=255)                                                                                   tAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTAGAtgacgttgacgt                                                             <  2:603956/150‑1 (MQ=255)                                                                                    aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCCGTTGAGAAAAAGGGGGGAGGGTTCATTACGTATTTGTTCCTGGTGGGGGGGCGTCCTTTTTTTTTTAATTGGGATTTCAGGGTCACGCTGGGTGACGTTGACGttt                                                           >  2:818910/1‑151 (MQ=255)                                                                                    aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATgacgttgacgtt                                                            <  2:986509/150‑1 (MQ=255)                                                                                       tttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACCGG‑GAGAAAAGGGGGGGAGGGTTCATTTCGTATTCGTTCCCGGGGGGGGGTCGTCCTTTTTTTTTTCATGGGGGTTTTTGGGCCACGCTGGGTGGCGGTTACGtttttt                                                        >  2:953958/1‑151 (MQ=255)                                                                                         tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTg                                                       <  2:736887/150‑1 (MQ=255)                                                                                          aGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgc                                                      <  1:301644/150‑1 (MQ=255)                                                                                             aaCGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTg                                                  <  1:373367/151‑1 (MQ=255)                                                                                                 cTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGGGGGAAAAGGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGTCATTTGTTTTTTAATGGTGGTTTTAGGTCCACGCTGGGTGGCGTTGACGGTTTTTTCGCTggt                                                >  1:990765/1‑147 (MQ=255)                                                                                                 cTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGCGAGAAAAAGGGGGGGGGGTACATTATGGATTACTTCCTGGGGGGGGCTCGTCATTTGTTTTTTATTGGGGATTTTAGGTTCACGCTGGATGACGTTGACGttttttttg                                                     >  1:618795/1‑141 (MQ=255)                                                                                                     tAAACCGTTTTCTGCTATGGCGG‑TTAACCGAAGGAAAAGGGGGGGAGGGTTCATTACGTATTCGTAACAGGGGGGGCGACGTCCTTTTTTTTTTTTTTGGGATTTCCGGTCCACGCTGGGTGGCGCTGACGGTTTTTTTGCTGAttaaatt                                          >  1:227254/1‑145 (MQ=255)                                                                                                      aaaCCGTTTTCTGCTATGGCGG‑TTACCGGAGGGGAAAGGGGGGGGGGGTACATTAAGTATTCGTGCCTGGGGGGGGCTAGTCATTTGTTTTTTTATGGGGATTTTAGGTCCCCCCTGGGTGACGTTTACGTTTTTTTCGCTGATAAGattt                                         >  1:971346/1‑151 (MQ=255)                                                                                                      aaaCCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtttt                                         <  1:735701/151‑1 (MQ=255)                                                                                                                       ggCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAata                        <  1:931145/151‑1 (MQ=255)                                                                                                                        gCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAatat                       <  2:667371/151‑1 (MQ=255)                                                                                                                        gCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAatat                       <  1:818990/151‑1 (MQ=255)                                                                                                                         cGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATAtt                      <  2:844257/151‑1 (MQ=255)                                                                                                                               aCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                 >  2:101590/1‑151 (MQ=255)                                                                                                                               aCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                 >  1:265163/1‑151 (MQ=255)                                                                                                                                ccGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGc                >  1:93620/1‑151 (MQ=255)                                                                                                                                    agagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat            >  2:385042/1‑151 (MQ=255)                                                                                                                                        aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt        <  2:171046/151‑1 (MQ=255)                                                                                                                                        aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt        <  1:140539/151‑1 (MQ=255)                                                                                                                                            gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc    <  1:798983/151‑1 (MQ=255)                                                                                                                                            gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCGGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc    <  1:218056/151‑1 (MQ=255)                                                                                                                                             ggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa   <  2:378241/151‑1 (MQ=255)                                                                                                                                              gggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAt  >  2:443886/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAT  >  NC_000913/1478181‑1478471

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