Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	38.3%	99.1 / 18.0	61	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (4/33);  new base A (23/0);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 7.50e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508514‑3508809                                                                                                                                                      |                                                                                                                                                    aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTc                                                                                                                                                   <  1:971737/151‑1 (MQ=255)  cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTca                                                                                                                                                  <  1:225849/151‑1 (MQ=255)  cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTca                                                                                                                                                  <  2:1009551/151‑1 (MQ=255)     tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑AGCGGGCGGGAAAAttaa                                                                                                                                                  >  2:149979/1‑146 (MQ=255)        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATCacacac                                                                                                                                             >  1:989059/1‑150 (MQ=255)          ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATCACACACggg                                                                                                                                          >  2:821628/1‑150 (MQ=255)            aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                        <  1:680335/151‑1 (MQ=255)                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGATCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGccccc                                                                                                                                   <  2:653697/151‑1 (MQ=255)                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATCACACACAGTAGccccc                                                                                                                                   >  1:447293/1‑151 (MQ=255)                  ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAAGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGAAAAAATTACACACGGTAGccccc                                                                                                                                   >  1:506100/1‑150 (MQ=255)                    agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATTAAACACGGTAGccccccc                                                                                                                                 >  1:782833/1‑148 (MQ=255)                      agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGTATGTAGGAAAGAAAAGGGCGGCCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGccccccgagc                                                                                                                              <  2:143723/151‑6 (MQ=255)                      agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGCGGGGAAAAA‑ACACAAACAGGAGcccc                                                                                                                                    >  2:508789/1‑145 (MQ=255)                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATc                                                                                                                              <  2:820400/150‑1 (MQ=255)                         tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATTACACACCGGAGCCCCCCCATc                                                                                                                              >  1:752585/1‑148 (MQ=255)                         tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGG‑GGAGCGCGGAAAAAATCACACACGGGAGccccc                                                                                                                                   >  1:330927/1‑143 (MQ=255)                           aagatGCACATGCTCATGGGCGAGGGTGTAACTCGTCAGATCCAAGATGATATTGTGTCTGAGGCTAACCGTAGGGAGCAAAGAAAAGAACGACCTTTTTGGGGAGGT‑GGCGCGCGAAAAAATTAACACACTTTAGCCCCCACATCAGa                                                                                                                           <  2:680141/146‑1 (MQ=255)                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                         <  1:696803/151‑1 (MQ=255)                                gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                    <  1:226525/151‑1 (MQ=255)                                  aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGGGCGGGAAAATTCAAACACGGGAGCCCCCACACCAGAAAACgggg                                                                                                                   >  1:799143/1‑150 (MQ=255)                                              gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCAACAGAAAGCGGGGGGGGCACAcc                                                                                                         >  1:949283/1‑148 (MQ=255)                                                  gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGGGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAGAAAACGGGGGGGGCCCAcc                                                                                                         >  1:716252/1‑144 (MQ=255)                                                          cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCCCACCGAGCTGCGTAGCGat                                                                                          >  1:154254/1‑151 (MQ=255)                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGGAAAATCACACACGCTGGCGCCCACATCAAAAAGCGTGGTGCGCGGACCGAGCCGCGTaaggaaaaacc                                                                                     >  1:941058/1‑141 (MQ=255)                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGC‑TGCGCGCGGGAAAAATCAAACAACGGAGGCCCCACAACAAAAAACGGGGGGCGGGCACCGAGCTCCGTAGCGaacat                                                                                       >  1:921617/1‑149 (MQ=255)                                                                ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAAATCACACACGGTAGCCCCCACAACAGAAAACGGGGGGGGCACACCGAGCTGCGTAGCGATCATccc                                                                                    >  1:2938/1‑151 (MQ=255)                                                                    cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                 <  1:235863/150‑1 (MQ=255)                                                                        aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                            <  1:979806/151‑1 (MQ=255)                                                                         cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCACCAGAAAAAGGGGGGCGCCCACCGCGCTGCGTAGCGATCATCCCTTGCttgt                                                                            >  1:796398/1‑147 (MQ=255)                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                        <  1:316382/151‑1 (MQ=255)                                                                                ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                                    <  1:781063/151‑1 (MQ=255)                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                   <  1:861513/151‑1 (MQ=255)                                                                                     ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAAATCACACACGGGAGCCCCCACCTCAGAAAGCGGGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTgg                                                               >  1:763328/1‑151 (MQ=255)                                                                                      tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGGGCGGGAAAAATCACCAACGGTAGGCCCCCCATCAGAAAGCGTGGGGGGCGCACCCAGCTGCGGAGCGATCATCCCTTGCTGGTAAACCAGCACCgg                                                               >  1:300412/1‑150 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                               <  1:169204/149‑1 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                             <  1:970678/151‑1 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGC‑GGCGCGCGGGAAAAATCACAAACCGTAACCCCCCCATCAGAAAACGGGGGGGGCCCCCCCCGCTGCGTAGCGAACATCCCTTGCTTATAAAACAACACttg                                                               >  2:1012127/1‑147 (MQ=255)                                                                                            aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAAATCACACACCGGAGCCCCCACATCAGAAAACGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGcacc                                                           >  1:42879/1‑148 (MQ=255)                                                                                              aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                      <  1:150248/151‑1 (MQ=255)                                                                                               ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAAATAACACACGGTAGCCCCCACATCCGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCGATAATCCCTTGCTGATAAACCAGCACTGGCCCCACtt                                                       >  2:700736/1‑149 (MQ=255)                                                                                                   agtgaggaAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAAATTAAACAAGGGCACCCCCCCAACAGAAAGCGTGGTGGGGACACCGCACTGCTTAGCCAGAATCCCTCGGTGATAAATCACCACCGGCGCCCCCTCGAgg                                                  >  2:116204/4‑150 (MQ=255)                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                 >  1:166660/1‑151 (MQ=255)                                                                                                     ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGGAAAAATCACACACCGGAGCCCCCACATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTGGCACCACCTCCCggg                                                 >  1:930976/1‑149 (MQ=255)                                                                                                           aagaaaagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtgg                                         <  2:184056/151‑1 (MQ=255)                                                                                                           aagaaaagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                          <  2:658454/150‑1 (MQ=255)                                                                                                             gaaaagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                       <  2:235649/151‑1 (MQ=255)                                                                                                              aaaagaGCGACCATTTTGTGGCGGT‑GGCGCGGGGGAAAAATTACACACCGGAGCCCCCCCCCCAGAAAGCGAGGGGCGCACCCCCCGCCGCCGACCGATCATCCCTTTCTGGTAAAACAACCCTGGCCCCCCCTCCCGGGGAtgggggg                                        >  1:388624/1‑148 (MQ=255)                                                                                                              aaaagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCTGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                      <  2:260190/151‑1 (MQ=255)                                                                                                               aaagaGCGACCATTTTGTGGCGGT‑GGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACCGGCACCACCTCGCGGGGAtggggg                                         >  1:235764/1‑147 (MQ=255)                                                                                                               aaagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgt                                     <  1:867824/151‑1 (MQ=255)                                                                                                                aagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                    <  2:505615/151‑1 (MQ=255)                                                                                                                aagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                    <  1:469087/151‑1 (MQ=255)                                                                                                                aagaGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                    <  2:658716/151‑1 (MQ=255)                                                                                                                            tttGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                         <  1:720171/150‑1 (MQ=255)                                                                                                                             ttGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                        <  2:394437/150‑1 (MQ=255)                                                                                                                              tgtgGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                       <  1:950226/150‑1 (MQ=255)                                                                                                                                  gcggt‑ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                     <  2:77288/148‑1 (MQ=255)                                                                                                                                  gcggt‑ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                     >  1:77353/1‑148 (MQ=255)                                                                                                                                               ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc      <  1:406462/151‑1 (MQ=255)                                                                                                                                                  aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt   <  1:822201/151‑1 (MQ=255)                                                                                                                                                   aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  >  2:591729/1‑151 (MQ=255)                                                                                                                                                   aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  <  1:986431/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508514‑3508809

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