Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,101,507	0	T	G	44.4%	119.3 / 11.8	72	intergenic (+77/‑183)	sodA/kdgT	superoxide dismutase, Mn/2‑keto‑3‑deoxy‑D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base T (9/31);  new base G (32/0);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCAT  >  NC_000913/4101360‑4101656                                                                                                                                                     |                                                                                                                                                        aCTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGtttt                                                                                                                                                      <  2:411717/151‑1 (MQ=255) aCTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgttt                                                                                                                                                      >  1:411422/1‑151 (MQ=255)    aCATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGttttt                                                                                                                                                     <  2:5834/149‑1 (MQ=255)    aCATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  1:727157/1‑151 (MQ=255)    aCATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  2:390583/1‑151 (MQ=255)    aCATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  1:39333/1‑151 (MQ=255)     cATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  1:82409/1‑150 (MQ=255)      aTTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  2:921390/1‑149 (MQ=255)        tAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGgtttttt                                                                                                                                                   >  2:317903/1‑147 (MQ=255)        tAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTtgtg                                                                                                                                                <  2:730865/150‑1 (MQ=255)          aaGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGAt                                                                                                                                              >  1:681086/1‑150 (MQ=255)               ttCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATTAAttt                                                                                                                                        >  1:941202/1‑151 (MQ=255)                   ggAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTc                                                                                                                                       <  2:484330/148‑1 (MQ=255)                   ggAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCaa                                                                                                                                     <  1:589961/150‑1 (MQ=255)                     aaCGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGGTCAATTTTaaaa                                                                                                                                   >  2:687122/1‑150 (MQ=255)                     aaCGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTaaaa                                                                                                                                   >  1:765887/1‑150 (MQ=255)                        gtggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCG‑AAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTCCAaaaaaa                                                                                                                               >  2:347531/1‑150 (MQ=255)                               aCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATg                                                                                                                        <  2:277118/151‑1 (MQ=255)                               aCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATaaaaaaaa                                                                                                                          >  1:91892/1‑149 (MQ=255)                                  gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGAt                                                                                                                      <  1:446559/150‑1 (MQ=255)                                  gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGAt                                                                                                                      <  1:166261/150‑1 (MQ=255)                                  gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTCAAAAAAAAAAAAATGATc                                                                                                                     >  1:791144/1‑151 (MQ=255)                                   ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATcc                                                                                                                    <  2:108745/151‑1 (MQ=255)                                       gAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGat                                                                                                                 >  1:362879/1‑150 (MQ=255)                                       gAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑AGGGTTTTTTTTGGTTAATTTCTaaaaaaaaaaaaataaacgat                                                                                                                 >  1:1004436/1‑141 (MQ=255)                                            agcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTTTGATCAATTTCTAAATAAAAAAAAAAAAACGataaaa                                                                                                             >  2:57540/1‑149 (MQ=255)                                                gCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataaaaataaaa                                                                                                       <  1:773291/151‑1 (MQ=255)                                                  aCGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACa                                                                                                     <  1:444935/151‑1 (MQ=255)                                                    gTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGGTTAATTTCAAAATAAAAAAAATGATCCGataaaaataaaa                                                                                                       >  1:213600/1‑147 (MQ=255)                                                      tttGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑AGGGTTTTTTTTGATCAATTTTAAAATAAAAACAAAGATCCGAAAAAAATAAAACa                                                                                                     >  1:298646/1‑147 (MQ=255)                                                           ggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCaa                                                                                          <  1:269440/151‑1 (MQ=255)                                                            gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGGTCAATTTCAAAAAAAAAACAAAGATCCGATAAAAATAAAACAGC‑G‑TTTCaa                                                                                          >  1:425735/1‑150 (MQ=255)                                                             cgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑CGGTTTTTTTGTTATTAATTTTTAAATAAAAAAAAAAATCCGATAAAAATAAAAAAAC‑G‑CTTCAAtt                                                                                        >  1:488454/1‑151 (MQ=255)                                                                      aatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATTAAATTCCAAAAAAAAACAAACATCCGATAAAAAAAAAAAAACGG‑TTTCAATT‑ATGGGGtt                                                                               >  2:676202/1‑151 (MQ=255)                                                                       atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGttt                                                                              <  2:968207/151‑1 (MQ=255)                                                                         cATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTg                                                                            <  1:954593/151‑1 (MQ=255)                                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAATGAGTCGTAAAAAATTAAAAAGC‑G‑ATTCCATTTATGTGGgttt                                                                             >  2:68035/1‑147 (MQ=255)                                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGAt                                                                          <  2:558916/150‑1 (MQ=255)                                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGC‑G‑TTTTAATTGATGTGGGTTTTAt                                                                          >  1:852917/1‑150 (MQ=255)                                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGGTCAATTTTTAAATAAAAAAAAAAGTCCGAAAAAAAAAAAACAGT‑G‑TTTCCATTTAGGAGGGTTTgtt                                                                          >  1:59935/1‑148 (MQ=255)                                                                             tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATc                                                                         <  1:26615/150‑1 (MQ=255)                                                                               ccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAAC‑G‑TTTTAATTGATGTGGGTTTTATTAc                                                                       >  2:643255/1‑150 (MQ=255)                                                                                cgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACtt                                                                     <  1:285913/151‑1 (MQ=255)                                                                                 gccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGGTCAATTTTAAAAAAAAAAAAAATAACCGATAAAAAAAAAAAAAC‑G‑TTTTAATTTATGTGGTTTTGATCACttt                                                                    >  2:39661/1‑151 (MQ=255)                                                                                 gccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGC‑G‑TTTTAATTGATGTGGTTTTTATCACttt                                                                    >  2:512702/1‑151 (MQ=255)                                                                                   ctgctgCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTCAAAATAAAAACAAAGATCCGATAAAAATAAAACAGC‑G‑TTTTAATTGATGTGGTTTTTATCACttttt                                                                  >  1:909854/1‑150 (MQ=255)                                                                                    tgctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGGTCAATTTTAAAAAAAAAACAATAATTCGATAAAAATAAAACAGC‑G‑TTTCCATTTATTTGGTTTTGATCACtttt                                                                   >  1:505360/1‑149 (MQ=255)                                                                                    tgctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGC‑G‑GTTCAATTGATGTGGTTTTTATCACTTTTAt                                                                 >  1:149999/1‑151 (MQ=255)                                                                                     gctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTa                                                                  <  2:150096/149‑1 (MQ=255)                                                                                     gctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTa                                                                  <  2:556581/149‑1 (MQ=255)                                                                                         cAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTtatttatt                                                            >  2:904156/1‑151 (MQ=255)                                                                                          aaTGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTa                                                           <  2:779800/151‑1 (MQ=255)                                                                                             gAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTaaa                                                         <  1:627554/150‑1 (MQ=255)                                                                                               ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                                      <  1:184992/151‑1 (MQ=255)                                                                                               ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATTAATTTTAAAAAAAAAACAAAGATCCGATAAAAATAAAAAAGC‑G‑TTTTAATTTATGTGGGTTTTATTACTTTTTTTTATTaaa                                                         >  2:290810/1‑148 (MQ=255)                                                                                               ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGA‑G‑TTTTAATTTATGTGGTTTTGATCACTTTTATTTATTaaa                                                         >  2:885762/1‑148 (MQ=255)                                                                                                 cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaat                                                    <  1:921403/151‑1 (MQ=255)                                                                                                 cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaat                                                    <  1:981568/151‑1 (MQ=255)                                                                                                  gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaatg                                                   <  2:621300/151‑1 (MQ=255)                                                                                                       ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTAt                                              <  2:629575/151‑1 (MQ=255)                                                                                                         ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑CGGGTTTTTTTTGGACAAATTTAAAAAAAAAAAAAAACACCGGTAAAAATAAAACAAC‑GATTTTAATT‑AAGTGGGTTTTATTACTTTTATTGATTAAAAGAATGTCTATCt                                            >  2:483419/1‑151 (MQ=255)                                                                                                           gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAGC‑G‑TTTTAATTTATGTGGTTTTTATTACTTTTATTTATTAAATGAAAGTCTATCtt                                           >  2:529276/1‑150 (MQ=255)                                                                                                           gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGC‑G‑GTTTAATTGATGTGGTTTTTATCACTTTTTTTGATTAAATGAATGTCTATCtt                                           >  1:596148/1‑150 (MQ=255)                                                                                                             atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGCTGGTTTTTTT‑TGATTTAATTTAAAAAAAAAACAAACAACCGGTAAAAATAAAAAAGC‑G‑GTTTAATTGATGTGGGTTTTGTTACTTTTTTTGATTAAAAGAATGGTTATCtt                                           >  2:211582/1‑148 (MQ=255)                                                                                                                 cAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAAAAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTTAATTTATGTGGGTTTTATCACTTTTATTTATTAAATAAATGTCTATCTTCGttt                                      >  2:718447/1‑149 (MQ=255)                                                                                                                   gCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                                  <  2:604241/151‑1 (MQ=255)                                                                                                                      tAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATc                                 <  1:531652/149‑1 (MQ=255)                                                                                                                      tAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCa                                <  2:200812/150‑1 (MQ=255)                                                                                                                        aaaaaTGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAaca                             <  2:586044/151‑1 (MQ=255)                                                                                                                        aaaaaTGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGC‑G‑TTTTAATTTATTTGGTTTTTTTAACTTTTATTTATTAAATGAATGTCTATCTTCGTTTCCATTAac                              >  2:768580/1‑150 (MQ=255)                                                                                                                         aaaaTGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGC‑G‑TTTCAATTGATGTGGTTTTTATCACTTTTTTTTATTAAATGAATGTCTATCTTTGTTTCCATCAaca                             >  1:178738/1‑150 (MQ=255)                                                                                                                         aaaaTGAACCATCGCCAACGGCGG‑TGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAAAGATCCGGTAAAAATAAAAAAGC‑G‑TTTTAATTTATGTGGTTTTTATCACTTTTTTTGATTAAAAGAATGTCTATCTTTGTTTCCATCAaca                             >  2:485191/1‑150 (MQ=255)                                                                                                                          aaaTGAACCATCGCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTTAATTTATGTGGGTTTGATCACTTTTTTTGATTAAATGAATGTCTATCTTCGTTTTCATCAacac                            >  1:477017/1‑150 (MQ=255)                                                                                                                            aTGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAACAGC‑G‑GTTTAATTTATGTGGTTTTTATTACTTTTTTTGATTAAATGAATGTCTATCTTTGTTTTCCTCAaca                             >  2:653500/1‑147 (MQ=255)                                                                                                                                 ccATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc                    <  2:727147/151‑1 (MQ=255)                                                                                                                                 ccATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc                    <  2:834617/151‑1 (MQ=255)                                                                                                                                      gCCAACGGCGG‑TGGGTTTTTTTTGATCAATTTCAAAATAAAAACAAAGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGGTTTTATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAAAACTGATGCTCCAtt                >  2:410129/1‑150 (MQ=255)                                                                                                                                        cAACGGCGG‑TGGGTTTTTTGTGGTCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAGC‑G‑TTTCAATTGATGTGGTTTTTATCACTTTTATTTATTAAATGAATGTCTATTTTCGTTTCCATCAAAACTGATGCTCCAtt                >  2:135977/1‑148 (MQ=255)                                                                                                                                             gcgg‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAAtt        <  1:898471/151‑1 (MQ=255)                                                                                                                                              cgg‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAAtt        >  2:445951/1‑150 (MQ=255)                                                                                                                                              cgg‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAAtt        >  2:201936/1‑150 (MQ=255)                                                                                                                                               gg‑tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTAc      <  2:681076/151‑1 (MQ=255)                                                                                                                                               gg‑tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTAc      <  1:337425/151‑1 (MQ=255)                                                                                                                                                g‑tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACg     <  2:449784/151‑1 (MQ=255)                                                                                                                                                    gtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCAt  >  2:27185/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                        ACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGG‑TGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGC‑G‑TTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCAT  >  NC_000913/4101360‑4101656

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