Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	42.7%	65.1 / 53.9	75	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (35/8);  new base A (0/32);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 5.06e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGG  >  NC_000913/3177603‑3177889                                                                                                                                                   |                                                                                                                                              tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaaa                                                                                                                                           >  2:861578/1‑150 (MQ=255) tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaaa                                                                                                                                           >  2:6486/1‑150 (MQ=255)  cACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAt                                                                                                                                         >  2:398692/1‑151 (MQ=255)  cACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAt                                                                                                                                         >  1:382604/1‑151 (MQ=255)   aCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                                                                                        >  2:163794/1‑151 (MQ=255)       tctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc                                                                                                                                    >  2:510105/1‑151 (MQ=255)              tcatcCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCcgccgc                                                                                                                               >  2:805192/1‑149 (MQ=255)                 tcCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTc                                                                                                                          >  2:1233671/1‑151 (MQ=255)                   ccGGCACCCATCCCCAGTAGCCAGGGGG‑TTTCGCTGGTGTCGTCCGCGGAAATCCGAATCTGCTCAACCACCACAACTTCACCACGCACGGGGTAAAAGGGTAGAAAGACTGCGGCGTGACCGCCCTAAAAAATTTCCCGCCGCACCTCAt                                                                                                                        <  2:83114/151‑1 (MQ=255)                   ccGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCa                                                                                                                         >  2:821106/1‑150 (MQ=255)                    cGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg                                                                                                                       >  2:791670/1‑151 (MQ=255)                    cGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg                                                                                                                       >  2:360881/1‑151 (MQ=255)                     gtcacccATCCCCAGTAGCCAAGGGT‑TTTCGCTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCACCACGCATTGGTTCAAGGGGTAGAAAGACTGCGGCGTGCCCGCGCTCAAAAATTCCCCGCCGCACCTCAt                                                                                                                        <  2:229852/144‑1 (MQ=255)                        aaCCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTc                                                                                                                   >  2:978330/1‑151 (MQ=255)                         aCCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGGGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTc                                                                                                                   >  2:817904/1‑150 (MQ=255)                             tctcCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTg                                                                                                              >  1:294680/1‑151 (MQ=255)                              ctcCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCCCCGCACCTCATGACTCATTTg                                                                                                              >  2:470944/1‑150 (MQ=255)                              ccccAGTAGCCAGGGGG‑TTTCGCTGGTGTCGTACGCGGCAACCCGAATCTGCCCAATCAGCCCAACTCCACCACGCACGGGGTAAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCAtt                                                                                                                <  1:1281802/146‑1 (MQ=255)                                  aGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                                                                         >  2:963348/1‑151 (MQ=255)                                    tAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                                                                                        >  2:988898/1‑150 (MQ=255)                                     aGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                       >  1:11348/1‑150 (MQ=255)                                     aGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                       >  2:650092/1‑150 (MQ=255)                                      accAAGGGGTTTTCGCTGGTTACGTACCCGCCAACCCAACTCTGCCAAACAAACACACCTCAAGCCCCCCCGGGATAAAAGGTTACAAAGCCACGGGCGTGACCGCGCCAAAAATTTTCCCGCCGCACCTCATGACTCATTTGccc                                                                                                           <  1:1028148/145‑1 (MQ=255)                                       ccAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACGGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa                                                                                                      <  2:66974/149‑1 (MQ=255)                                        cAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                    >  1:822579/1‑150 (MQ=255)                                        cAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                    >  2:656587/1‑150 (MQ=255)                                           gggg‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                    <  2:2600/147‑1 (MQ=255)                                            ggg‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCACTGGGTAAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                               <  2:608125/151‑1 (MQ=255)                                             gg‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                               >  2:687503/1‑150 (MQ=255)                                                tttGGCTGGTGTTGTGCGAGCAAATCCAATTTTCCTAAACCACCACACCTTCATCACGCGTCGGATCAAAGTGTAAAAAGCTTGAGGGGTCCCCGTGTTAAAAATTTACCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                                                                              <  1:285114/149‑1 (MQ=255)                                                ttccGCTGGTGTCGTACGCGGCAACCCGACTCTGCTCAATCAGCACACCTCCATCACGCCCGGGTTAAAAGGGTAGCAAGCCTGCGGCGTGACCGCGCTCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt                                                                                             <  2:182048/147‑1 (MQ=255)                                                   cGCTGGTTTGGTACGCGGAAATCCGAATCTTCTCAATCAGCACAACTTCATCCCGCCTTGGGCCAAAGGGTAGCAAAGCTGGGGCGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                           <  2:881672/149‑1 (MQ=255)                                                        gtgtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGCCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATCGACGATGACGAAATc                                                                                    >  1:1300395/1‑151 (MQ=255)                                                           tCGGACGCGGAAATCCGAATCTCCAAAACCACCACAACTTCATCACGCACGGGGTCAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                   <  1:990720/149‑1 (MQ=255)                                                            cGTCCGCGGCAATCCGAATCTGCTAAATCAGCACAACTCCATCACGCACGGGGTCAAAGGGTAGCAAGACTGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                <  1:766234/151‑1 (MQ=255)                                                                 gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                               >  1:724455/1‑147 (MQ=255)                                                                 gcgACAACCCGAATCTGCTCAATCAGCCCATCTTCACCCCGCGCTGGGTCAAGGGGTACCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTGCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                            <  1:647378/150‑1 (MQ=255)                                                                  cgGAAACCCAAATCTGCTCAATCAGCACACCTCCACCACGCACTGGGACAAAGGGTAGCAAGACGGCGGCGGGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                           <  1:349387/150‑1 (MQ=255)                                                                    gaaaTCCGAATCTGCACAAACAGCCCACCTCCACCACGCACGGGGACAAAGGGAAGCAAGACTGCGGCGTGCCCGCGCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                               <  1:308548/142‑1 (MQ=255)                                                                      aaTCCGAATCTGCACAACCACCACAACTCCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                         <  1:235747/148‑1 (MQ=255)                                                                         ccGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTTAAAGGGCAGCAAGACTGAGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAACTTTAAAGATCTAAt                                                                      >  1:1177901/1‑148 (MQ=255)                                                                         ccGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                    >  2:885304/1‑150 (MQ=255)                                                                          cGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                    <  1:1181864/149‑1 (MQ=255)                                                                               cTGCTCAATCCGCCCACCTCCATCCCGCAGTGGTTCAAAGGTTAGCAAGACTGCGGCGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                 <  2:841923/147‑1 (MQ=255)                                                                                 gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc                                                           >  1:425937/1‑151 (MQ=255)                                                                                    cAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                        >  1:74584/1‑151 (MQ=255)                                                                                     aaccAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGGGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgc                                                         <  1:1009607/146‑1 (MQ=255)                                                                                     aaTCAGCACAACTTCATCACGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg                                                          <  2:508130/148‑1 (MQ=255)                                                                                             cacctccaccaCGCCCGGGGCAAAAGGTAAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaaa                                                   <  2:722659/138‑1 (MQ=255)                                                                                                   catcaCCCACTGTGTCAAAGGGTAACAAGACTGCGGCGTTCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAACTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                            <  1:835051/148‑1 (MQ=255)                                                                                                      cacgcacgGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                       <  2:565014/150‑1 (MQ=255)                                                                                                        cGCACGGGTTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                       <  1:791640/148‑1 (MQ=255)                                                                                                         gCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                   >  2:710451/1‑151 (MQ=255)                                                                                                                 tCATAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGTCGCGATACAGTGTTTCTCGTGCAATAATTTCt                            >  1:885579/1‑150 (MQ=255)                                                                                                                  cAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                          >  1:731162/1‑151 (MQ=255)                                                                                                                   aaaGGGTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCt                            <  1:510084/148‑1 (MQ=255)                                                                                                                   aaaGGGGAGCAAAACGGCGCCGTGCCCGCGCTCAAAAATTCCCCGCCGCTCCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                               <  2:578726/145‑1 (MQ=255)                                                                                                                    agggggAGCAAGCCTGGGGCGTGACCGCCCTAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACAt                        <  1:830839/149‑1 (MQ=255)                                                                                                                     aGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                       >  2:509212/1‑151 (MQ=255)                                                                                                                      gggTAGCAAGCCTGCGGCGTGCCCGCGCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                       <  1:417624/150‑1 (MQ=255)                                                                                                                      gggTAGAAAGACTGCGGCGTGCCCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                       <  1:978273/150‑1 (MQ=255)                                                                                                                        gTAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtt                    <  1:172782/151‑1 (MQ=255)                                                                                                                         tAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtt                    <  2:952905/150‑1 (MQ=255)                                                                                                                          accAAGACGGCGCCGTGCCCTCTCACAAAAATTACCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                         <  1:992725/142‑1 (MQ=255)                                                                                                                            cAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttttt                 >  1:709831/1‑150 (MQ=255)                                                                                                                             aaagctcgggcgTCCCCACGCTCAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                       <  1:821074/135‑1 (MQ=255)                                                                                                                                acgggggcgTGCCCGCCCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa             <  1:1280920/145‑1 (MQ=255)                                                                                                                                aCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa             >  1:72751/1‑150 (MQ=255)                                                                                                                                   gcggcgTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATg         <  2:20698/151‑1 (MQ=255)                                                                                                                                   gcggcgTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAAt          <  2:980706/150‑1 (MQ=255)                                                                                                                                     ggcgTGACCGCTCTCAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAAt          <  2:210905/148‑1 (MQ=255)                                                                                                                                         tccccGCGCTAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg    <  1:710421/147‑1 (MQ=255)                                                                                                                                          gACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACggg  >  2:1208479/1‑151 (MQ=255)                                                                                                                                          gACCGCGATCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACgg   <  2:540450/150‑1 (MQ=255)                                                                                                                                            ccGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAAc     <  2:53889/146‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGG  >  NC_000913/3177603‑3177889

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