Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,077	0	C	A	38.8%	114.6 / 30.0	67	D515Y (GAT→TAT)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (33/8);  new base A (0/26);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 8.30e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCAT  >  NC_000913/1575928‑1576221                                                                                                                                                       |                                                                                                                                                 aCGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATc                                                                                                                                                  >  2:940903/1‑150 (MQ=255) aCGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCa                                                                                                                                                 <  1:268884/151‑1 (MQ=255)    cgcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGa                                                                                                                                              <  2:154959/151‑1 (MQ=255)       aaTAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCATAGTTTTATTTCCTGAAACCGATCCCATCCAGCTATTTCGGATATCACGTAATCCCAGATCAAGAata                                                                                                                                           <  2:911320/151‑1 (MQ=255)          aaCGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTc                                                                                                                                        <  2:112412/151‑1 (MQ=255)                 ttACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCAt                                                                                                                                 >  1:496248/1‑151 (MQ=255)                    ccGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGt                                                                                                                              >  1:718996/1‑151 (MQ=255)                      gATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTa                                                                                                                             >  2:732141/1‑150 (MQ=255)                        ttttCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACg                                                                                                                          >  2:492718/1‑151 (MQ=255)                         tttCTGCTGCAATCCCATTGCCAGTTCCTCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACgg                                                                                                                         >  2:835988/1‑151 (MQ=255)                             tgctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTtat                                                                                                                     >  1:283104/1‑151 (MQ=255)                             tgctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGACATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTtat                                                                                                                     >  2:1085916/1‑151 (MQ=255)                                 gcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAAcc                                                                                                                 >  2:782132/1‑151 (MQ=255)                                      cccATTGCCAGTTCATCGTTATAAGGCTTTTTCAAA‑TCTTGAAAACGCGTCAGATTTTTTTTCCCTGATACCGATTCCGTCCAGTTATTGCGGATATCACGTAATCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACctgc                                                                                                              <  1:207123/149‑1 (MQ=255)                                       ccATTGCCAGTTCATCGTTATAGGGCTTTTTAAAA‑TCCGGAAACCGCGTAAGATTTTTATTCCCTGATACCGTTTCCGCCCAGTTATTGCGGATATCCCTTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                           <  1:791522/151‑1 (MQ=255)                                        cATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                           >  2:888847/1‑150 (MQ=255)                                           tGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                        >  2:801123/1‑150 (MQ=255)                                             ccttttCATTTTTATAGGGCTTTTTCAAA‑TCTTGATACCGCGTCAGATTTTTTTTCCCTGATCCCGTTTCCGCCCAGTTATGGCGAATATCACGTAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTaa                                                                                                         <  1:306689/143‑1 (MQ=255)                                             ccattgcATCTTTATAAGGTTTTTTAAAA‑TCCTGATAACGCTTCAGATTTTTATTCCCTAATCCCGATCCCGTCCAGCTATGGCGGATATCACGTAACCCCATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                        <  1:848956/142‑1 (MQ=255)                                                 ttctttgtaataGGGCTTTTAAAAA‑TCCTTATAAGGCGAAAGATTTTTATTCCCTAATACGAATCCCCTCCACTTATGGGGGATATCACGAAATCCAAAATAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGTTGTAATCATTg                                                                                                   <  2:11300/140‑1 (MQ=255)                                                   cATCGTTA‑AAGGCCTTTTTAAATTCTTAAAAACGCGTCGGTTTTTTATTCCCGGATCCCGTTCCCGCCCAGCTATGGCAAATATCACTAAACCCAATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAgg                                                                                                <  2:1058358/150‑1 (MQ=255)                                                     tttttATAAGTCTTTTTCAAC‑TCCTAATAACGCCTCAGTTTTTTATTCCCTAATCCCGTTCCCTTCCACCTATTGCGGAACTCCCGACATCCAAAATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                  <  1:342631/143‑1 (MQ=255)                                                       gTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                            >  2:825713/1‑150 (MQ=255)                                                          ataAGGCTTTTTAAAA‑TCCTAATAACGCGTCAGATTTTTTTTACCTAATACCGATTCCTCCCAGCTATTGCGGATATCACAAACTCCTATAACAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAt                                                                                          <  2:289852/149‑1 (MQ=255)                                                                 gTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                  >  2:986299/1‑150 (MQ=255)                                                                 gTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                  >  2:167129/1‑150 (MQ=255)                                                                 gTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                  >  1:183559/1‑150 (MQ=255)                                                                  ttttCAAA‑TCCTGATAACGCTTCAGATTTTTTTTCCCTGATACCGTTTCCGTCCAGCTATTGCGGATATCACGTAACCCCATATAAAGATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                  <  1:443920/149‑1 (MQ=255)                                                                    ttCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAatat                                                                              >  1:570869/1‑151 (MQ=255)                                                                       aaa‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATccc                                                                           >  2:338473/1‑151 (MQ=255)                                                                       aaa‑TCCTAATACCGCGTGAGTTTTTTTTTCCCTGATGCAGTTCCCCTCCAGTTATTGCGAATATCCCGTACCCCCATTTAAAGATTATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATccc                                                                           <  1:643552/151‑1 (MQ=255)                                                                        aa‑TCCTGATAACGGGTCAGATTTTTATAACCTGATACCACTTCCGTCCAGTTATGGGGGATATCACGTACCCCAAAATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATccc                                                                           <  2:325614/150‑1 (MQ=255)                                                                             cttaTAACTCGTCAGAGTTTTATTACATGATTCCGTTTCTCCCCCGTTATGGAGGTAATCTCCACCCCCCTTATAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTc                                                                       <  2:438903/147‑1 (MQ=255)                                                                                aTACCGCGTCAGATTTTTATTCCCTAATCCCGATCCCGTCCAGTTATTGCGAATATCACGTAACCCAATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTcc                                                                      <  1:540583/148‑1 (MQ=255)                                                                                     gctttatatttttATTACCTAATACCGCTTCCTTCCAGCTATTGCGAATATAATGTACTCCCATATCAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATaa                                                               <  2:743873/141‑1 (MQ=255)                                                                                      gtacAGATTTTTTTTACCAAATACCGATTCCTCACAGTATTTGCGAATACCAGGTATCCCAAAATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAg                                                             <  2:723263/148‑1 (MQ=255)                                                                                      cttCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCACATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATaa                                                               >  1:402836/3‑149 (MQ=255)                                                                                       gTCAGATTTTTATTCCCTAATCCCTTTCCCGCCCGGCTTTTGCGGTTACCACGTAACCCAATTTAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATaaa                                                              <  1:235119/149‑1 (MQ=255)                                                                                        tCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAgcg                                                           >  1:975681/1‑151 (MQ=255)                                                                                          agatTTTTATTACCTGATCCCTATCCCGTCCAGCTATTGCGGATATCACGTAACCCCATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAgcg                                                           <  1:963630/149‑1 (MQ=255)                                                                                          agagTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGcgcg                                                         >  1:848555/1‑151 (MQ=255)                                                                                           gagTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                        >  1:289489/1‑151 (MQ=255)                                                                                            tgTTTTTTTACCTAATACCGATGCCTCCCGGTTATAGGACAAACTACGAAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                        <  2:469177/149‑1 (MQ=255)                                                                                               tttATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgggg                                                      >  2:447534/1‑149 (MQ=255)                                                                                                ttattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCggggcg                                                    >  2:1054687/1‑148 (MQ=255)                                                                                                     aCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATtg                                              >  2:849796/1‑151 (MQ=255)                                                                                                      ccTAATCCCGATCCCTCCCAGCTATGGCGAATATCACGTAACCCCATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGatat                                                <  1:898964/148‑1 (MQ=255)                                                                                                        tGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTgg                                           >  2:133612/1‑151 (MQ=255)                                                                                                           tcccGCTTCCCTCCAGCTATTGCGAATATCACTAAACCCCATATCAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTgg                                           <  2:79217/146‑1 (MQ=255)                                                                                                             catatcccctCCAGTTTTTGCGAATATCAAATACCCCCATATCAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGt                                          <  2:22940/138‑1 (MQ=255)                                                                                                              cGATCCCGCCCAGTTATGGGGGATATCACGTAACCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                       <  1:347025/149‑1 (MQ=255)                                                                                                               gATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                     >  2:1080668/1‑150 (MQ=255)                                                                                                               gATCCCGCCCAGCTATTGCGGATATCACGTAACCCAATATAAAAATTTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACa                                    <  1:115310/151‑1 (MQ=255)                                                                                                                  tCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                     >  2:135131/1‑147 (MQ=255)                                                                                                                  tCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                     >  2:444885/1‑147 (MQ=255)                                                                                                                   ccGTCCATCTATTGGGAATATCACGTAATCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                       <  2:783563/144‑1 (MQ=255)                                                                                                                   ccGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                                 >  1:615503/1‑150 (MQ=255)                                                                                                                    cGTCCCGCTATGGCTAAAACCACAAAATCCAAAATAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAg                               <  2:435389/151‑1 (MQ=255)                                                                                                                     gTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                              >  1:544493/1‑151 (MQ=255)                                                                                                                       ccAGTTATTGCGGATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                              <  1:941210/149‑1 (MQ=255)                                                                                                                             aTTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcat                        >  1:554676/1‑149 (MQ=255)                                                                                                                               ggCGTATATCACATAACCCCATATAAAGAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTg                            <  2:227209/142‑1 (MQ=255)                                                                                                                                gCGAATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtca                   <  1:986619/151‑1 (MQ=255)                                                                                                                                   tatatCACGTATTCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcata                       <  1:316491/143‑1 (MQ=255)                                                                                                                                         aCGTAATCCCCTATCAAGTATCTTCCCGTCACAGGAACGGTTATAACTTCATGTAATCATAGACGCTTGACTAGTATAATCATCCCATTCCGTCCTAGAGCGCGGGGAGATATTGTGGTTTGATAGATAGTTGTCTTAGCCATAACGCAc           >  2:663733/1‑150 (MQ=255)                                                                                                                                              aTCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCggg     >  1:378514/1‑151 (MQ=255)                                                                                                                                                cccATATAAAGATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCg       <  1:200895/147‑1 (MQ=255)                                                                                                                                                 ccATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCAt  >  2:73559/1‑151 (MQ=255)                                                                                                                                                     aacAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtcata                 <  1:1036530/130‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 ACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAA‑TCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCAT  >  NC_000913/1575928‑1576221

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