Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,488,209	0	G	T	39.5%	142.7 / 30.1	81	G298C (GGC→TGC)	lptG	lipopolysaccharide export ABC permease
Reads supporting (aligned to +/- strand):  ref base G (17/32);  new base T (32/0);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 1.90e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCT  >  NC_000913/4488064‑4488359                                                                                                                                                  |                                                                                                                                                         ctctATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGGGGTGATGCTGATGGGGCTGTCGGTTATTTTTTGcc                                                                                                                                                       >  1:671019/1‑149 (MQ=255)   ctATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCAc                                                                                                                                                    >  1:486804/1‑150 (MQ=255)    tATCAGTGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGccccc                                                                                                                                                    >  2:505246/1‑147 (MQ=255)    tATCAGCGGTTTGCACAACTATGTGAAGTATCTAAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTgg                                                                                                                                                         <  2:27638/144‑1 (MQ=255)     aTCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTg                                                                                                                                                  <  1:471676/150‑1 (MQ=255)     aTCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGGGGTGATGCTGATGGGGCCGTCGTTTCTCTTTTGccccc                                                                                                                                                    >  2:1402910/1‑146 (MQ=255)          cGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCCCTgcgtag                                                                                                                                             >  2:1268497/1‑150 (MQ=255)           ggTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTTATGCTGGTGGCGCTGTCGCTTCTTTTTTGCCCACTgcgtagcg                                                                                                                                           >  2:598164/1‑151 (MQ=255)           ggTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGCGATGATGCCGGTGGCGTCGTCGTTCCTCTTTTGGCCCCTgcgt                                                                                                                                               >  2:124113/1‑147 (MQ=255)            gTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATTGGGCTGGCGTTCATCTTTTGCCCCCTgcgtagcgg                                                                                                                                          >  2:1000436/1‑150 (MQ=255)            gTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGGGCTGTCGTTTATCTTTTGCCCCCCgcggag                                                                                                                                             >  2:1297950/1‑148 (MQ=255)              ttGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTgcgtagcgta                                                                                                                                         <  2:868225/150‑1 (MQ=255)              ttGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTgcgtagcgta                                                                                                                                         <  1:1298014/150‑1 (MQ=255)                 cacaACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCCGGCGGTCTTCTTTTGGCCCCTCCGGAGCGTAc                                                                                                                                        >  2:1372119/1‑148 (MQ=255)                   caACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCCCTGCGCCGGGTAcc                                                                                                                                       >  1:233813/1‑147 (MQ=255)                   caACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGGTGGTCATCTTTTGCCCACTGCGTAGCGTAcc                                                                                                                                       >  2:974153/1‑147 (MQ=255)                     aCTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGTTGGCGCCGTCGTTCATCTTTTGCCCCCCGCGGAGCGGACCGaatgg                                                                                                                                 >  2:118431/1‑147 (MQ=255)                       tATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCCCTGCGTAGCGTACCGATGGgcg                                                                                                                               >  1:1138603/1‑151 (MQ=255)                        aTGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGgc                                                                                                                                <  2:347397/149‑1 (MQ=255)                          gtgAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGgcgtgc                                                                                                                            <  1:893302/151‑1 (MQ=255)                              aGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGggcgggcg                                                                                                                           >  1:445669/1‑148 (MQ=255)                               gTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGGCCACCGCGTAGCGTACCGATGggcgggcggg                                                                                                                         >  1:40448/1‑147 (MQ=255)                                 aTCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGCTGGCGCCGTCGTTCATTTTTTGCCCCCTGCGTAGCGTACCGGTTGGCGTGGGTGTgg                                                                                                                      >  2:505418/1‑150 (MQ=255)                                  tCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGGGTGCGTGTgg                                                                                                                      >  2:113458/1‑149 (MQ=255)                                     gAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTc                                                                                                                    <  1:1202562/148‑1 (MQ=255)                                      aaGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTTTGCCCACTGCGTAGCGGACCGATGGGCGTGCGTGTGGTCAcc                                                                                                                 >  1:377513/1‑150 (MQ=255)                                      aaGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCAcc                                                                                                                 <  2:1270725/150‑1 (MQ=255)                                       aGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGCGGAGATGGTGTTGGGGGTGGTGTTTATCTTTTGCCCCCTGCCGTGCGTAACGACGGgggggcg                                                                                                                           >  2:942588/1‑139 (MQ=255)                                          cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGGGCTGTCGTTCATCTTTTGCCCCCTGCGTAGCGTACCGATGGGGGGGGGTGTGGTCACCGGGAt                                                                                                            >  2:1147703/1‑151 (MQ=255)                                          cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTa                                                                                                             <  1:545288/150‑1 (MQ=255)                                           gAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑c                                                                                                          <  1:278978/151‑1 (MQ=255)                                               ggTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGTTGATGGTGAAGGCGCTGTCGTTCTTTCTTGGCCCACCGCGGAGCGTACCGATGGGCGTGCGTGTGGTCACCGGGAT‑CAGt                                                                                                       >  1:1408966/1‑150 (MQ=255)                                               ggTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATTGCGCTGTCGTTCCTCTTTTGCCCCCTGCGTAGCGTACCCATGGGCGGGCGTGTGGTCACCGGGAT‑CCGtt                                                                                                      >  1:714831/1‑151 (MQ=255)                                               ggTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTTTTCTTTGGGCCACTGCCTAGCGAACCGATGGGGGTGCGTGGGGTCAcccgcgtcag                                                                                                         >  1:523250/1‑141 (MQ=255)                                                 tCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCTCTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGttt                                                                                                     <  2:743982/150‑1 (MQ=255)                                                  cAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCg                                                                                                   <  1:707603/151‑1 (MQ=255)                                                   aGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTTTCGTTCCTCCTTTTCCCCCCGCGTTGCGTTCCCCTTGGCGTGCGTGTGGTCCCCGGttt                                                                                                            >  2:662903/1‑140 (MQ=255)                                                   aGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGGCCCCCGCGTAGCGTGCCGATGGGGGGGGGGGTGGGCACCGGGAT‑TAGTTTTgg                                                                                                  >  2:1076204/1‑151 (MQ=255)                                                      aTGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGGTGCTGAATGCGCCGGCGTTCATCCTTTGCCCGATGAGTAGCCTACCGGGCGGTGGGCGTGGGGGCACCGGTAT‑CGGTTTCg                                                                                                   >  1:150065/1‑147 (MQ=255)                                                        gCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGtttt                                                                                              <  2:930315/150‑1 (MQ=255)                                                           ggACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGGAT‑CAGTTTCGGGTTTGtct                                                                                          >  1:344320/1‑151 (MQ=255)                                                               gTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGGTCACCGGTAT‑CAGTTTCGGTTTTGtctt                                                                                         >  1:122574/1‑148 (MQ=255)                                                               gTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGACAAGATGCTGATGGCGGTGGCGTTCCTCTTTGGCCCAATGCGTAGCGTAACGAAGGGCGCGTGTGTGGGCCCCCGAAC‑CAGTTTCTGCTTTCtctt                                                                                         >  1:273171/1‑148 (MQ=255)                                                                  aTCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATTCTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGGTCACCGGTAT‑CAGGTTCGGTTTTTtcttct                                                                                       >  1:1137588/1‑147 (MQ=255)                                                                  aTCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACg                                                                                    <  2:750467/150‑1 (MQ=255)                                                                  aTCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGt                                                                                   <  1:786313/151‑1 (MQ=255)                                                                  aTCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGCCGTGATGCCGATGGCGCTGTTGGTTATCTTTGGCCCACTGCCAAGCGTTCCGATGGGCGTGCGCGGGGCGAACGGTGT‑CCGTTTTGGGTTTGTTTTTTAc                                                                                     >  1:754355/1‑149 (MQ=255)                                                                       cTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTg                                                                               <  2:1001378/150‑1 (MQ=255)                                                                            cATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGTTGATTCTGTAGGCGCCGTAGCTCCTCTTTGGCCCCCCGCGCCGCGTACCGCTCGGCGCTCGGGTGGTCACCGGCATCCAGTTCGGGTTTGGTC‑TCTACGtggggg                                                                              >  2:11665/1‑141 (MQ=255)                                                                             aTGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGGGCTGTCGTTCATCTTTTGGCCACCGCGTAGCGTACCGGTGGGCGTGCGTGTGGTCCCCGGTAT‑CCGTTTCGGTTTTGTCTTCTTCGTACTGGAcc                                                                           >  1:1042848/1‑148 (MQ=255)                                                                              tgtgGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGGCACCGGTAT‑CAGTTTCGGGTTTGTTTTCTACGTACTGGAcc                                                                           >  1:1419279/1‑147 (MQ=255)                                                                               gtgGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAt                                                                       <  2:594103/150‑1 (MQ=255)                                                                                tgGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGGGCACCGGGAT‑CAGTTTCGGTTTTGTCTTCTTCGTACTGGGCCAGATC‑t                                                                    >  1:368024/1‑151 (MQ=255)                                                                                  gAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTc                                                                  <  1:1362091/151‑1 (MQ=255)                                                                                    gCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCg                                                                 <  2:612376/150‑1 (MQ=255)                                                                                       aaaTCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGccc                                                             <  2:1101837/151‑1 (MQ=255)                                                                                       aaaTCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGccc                                                             <  1:878629/151‑1 (MQ=255)                                                                                        aaTCTTCCAGCCGCTATCTGTGGCGGTGGTGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGGACCGATGGGGGGGGGGGTGGTCACCGGTAT‑CAGTTTTGGGTTTGTTTTTTTCGCGCTGGGCCAGATC‑TTCGGccc                                                             >  2:970522/1‑150 (MQ=255)                                                                                        aaTCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTTATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGGGCACCGGGAT‑CAGTTTCGGGTTTTTCTTTTACGTACTGGACCCGATC‑TTCGGcccc                                                            >  1:145938/1‑150 (MQ=255)                                                                                          tCTTCCAGCCGCTATCTGTGGCGGTGATGATGGTGATGGGGCTGTCGTTCATCTTTGGCCCACTTCGGAGCGGAACGATGGGCGGGCGTGTGGGCACCGGGAT‑CAGTTTCGGGTTTTGCTTCTACGTTCTTGGCCAGATC‑TTCGGGCCGc                                                           >  1:596503/1‑150 (MQ=255)                                                                                           cTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGGAT‑CAGTTTCGGTTTTTTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCt                                                          >  1:10111/1‑150 (MQ=255)                                                                                             tCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCACACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCt                                                          <  1:150104/148‑1 (MQ=255)                                                                                              ccAGCCGCTATCTGTGGCGGTGGTGATGCTGATGGGGCCGGTGTTCATCTTTTGCCCACCGCGGTGCGGACCCATGGGCGGGCGTGGGGTCACCGGGAT‑CAGTTTTGGGTTTTTTTTTTACGTAATGGACCAAATCTTTCGGCCC‑CTGAc                                                       >  2:1108873/1‑150 (MQ=255)                                                                                              ccAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGAc                                                       <  2:715617/150‑1 (MQ=255)                                                                                              ccAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGAc                                                       <  1:1435681/150‑1 (MQ=255)                                                                                               cAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGGTCACCGGGAT‑CAGTTTCGGGTTTGTCTTCTACGTACTGGACCAAAAC‑TTCGGCCCGCTGAc                                                       >  1:957820/1‑149 (MQ=255)                                                                                                aGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTTTCTTCTTCGTACTGGACCAGATC‑TTCGGCCCGCTGACggt                                                    >  1:1229121/1‑149 (MQ=255)                                                                                                           gtgGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTAt                                             <  1:949180/147‑1 (MQ=255)                                                                                                             ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCAt                                        <  2:1111822/150‑1 (MQ=255)                                                                                                             ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCAt                                        <  1:1140111/150‑1 (MQ=255)                                                                                                                 gTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCccgc                                   >  1:1445972/1‑151 (MQ=255)                                                                                                                           tGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCgg                         >  1:301625/1‑151 (MQ=255)                                                                                                                             aTGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcg                       >  2:706440/1‑151 (MQ=255)                                                                                                                                gcgcTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCAc                    <  2:1010537/151‑1 (MQ=255)                                                                                                                                  gcTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTg                  >  2:691543/1‑151 (MQ=255)                                                                                                                                     gTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTCCTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTg               <  2:1254081/151‑1 (MQ=255)                                                                                                                                     gTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGtt                <  1:302283/150‑1 (MQ=255)                                                                                                                                          tCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCaa           <  1:1111795/150‑1 (MQ=255)                                                                                                                                           cATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTAACCGGTAT‑CAGCTTCGGTTTTGTCTTATACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTAAGGCATGCCGCCGATCATCGGCCCACTGTTGccca           >  1:810443/1‑147 (MQ=255)                                                                                                                                             tCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAgcgc       >  2:246475/1‑151 (MQ=255)                                                                                                                                                 tGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCCCCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAg    <  2:933841/150‑1 (MQ=255)                                                                                                                                                 tGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAg    <  1:704966/150‑1 (MQ=255)                                                                                                                                                  ggCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGct  >  1:726874/1‑151 (MQ=255)                                                                                                                                                  ggCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGc   <  2:1344308/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                         CTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTAT‑CAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATC‑TTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCT  >  NC_000913/4488064‑4488359

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