Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,397	0	G	A	43.5%	76.3 / 43.0	62	E25K (GAA→AAA)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (0/27);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 1.99e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCGGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTT  >  NC_000913/2421254‑2421547                                                                                                                                                 |                                                                                                                                                        cGTTCGGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTaa                                                                                                                                                  <  1:305418/151‑1 (MQ=255)  gtacGGATAGTTAACAGAAACAGGTTCCCGTCTATTATTTCATCCTAAGTAAAACAACTGAGACCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAAAC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTaa                                                                                                                                                  <  2:677234/147‑1 (MQ=255)   ttCGGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAAtt                                                                                                                                                    <  1:962774/147‑1 (MQ=255)      ggATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAAcc                                                                                                                                             >  2:516112/1‑151 (MQ=255)       gATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCg                                                                                                                                            <  1:201877/151‑1 (MQ=255)          aGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCa                                                                                                                                         <  1:353785/151‑1 (MQ=255)               aCAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGata                                                                                                                                    <  1:64714/151‑1 (MQ=255)                     aCAGGTTCTCGTGTATTATTTCAGCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTa                                                                                                                             <  1:1015675/151‑1 (MQ=255)                          ttCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAAt                                                                                                                        >  2:186838/1‑151 (MQ=255)                            ctcGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAAtgt                                                                                                                      >  2:383990/1‑151 (MQ=255)                              cGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGa                                                                                                                    <  2:168561/151‑1 (MQ=255)                               gtgtATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGAc                                                                                                                   <  1:803063/151‑1 (MQ=255)                                 gtATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACg                                                                                                                  >  2:928216/1‑150 (MQ=255)                                      attTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCAc                                                                                                            >  1:835046/1‑151 (MQ=255)                                        ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACt                                                                                                           <  1:527961/150‑1 (MQ=255)                                          cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTAc                                                                                                         <  1:315873/150‑1 (MQ=255)                                          cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TCATCAATGTGACGATCCACTAc                                                                                                         <  2:544649/150‑1 (MQ=255)                                            tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTAcccc                                                                                                      >  1:867514/1‑151 (MQ=255)                                             ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTAcccc                                                                                                      >  1:1004678/1‑150 (MQ=255)                                             ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCccg                                                                                                     <  2:384376/151‑1 (MQ=255)                                              cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCccgc                                                                                                    >  2:172082/1‑151 (MQ=255)                                                   aaaaaCAACGAACACCCAGAAAGGCCACACCCTGCCTATTTTATTCATAAAAAGACCC‑TTCTTTTGCGTACTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTAAGGATATTG‑TTATCAATGTGACGACCCACTCCCCCGCCGATa                                                                                               <  1:832474/150‑1 (MQ=255)                                                            ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                       >  2:86169/1‑150 (MQ=255)                                                             gagaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAGCCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCAc                                                                                     >  1:700205/1‑151 (MQ=255)                                                                 aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                 >  1:1037779/1‑151 (MQ=255)                                                                   cTGCAATGGCACAACCTCCCGTTATTATTGGTATAAACAACC‑TTCGTTTGCGTACGGTTATCGGAATTAAGGAAGAAGCAGTTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCGGCCGATAAACGGGGCACTTGc                                                                                 >  1:1015422/1‑149 (MQ=255)                                                                       aatgcaacacCCTGCCTTTTTTATTTTTAAAAAAACCC‑TTCTTTTGCTTACTTTTAGCGGAATTAAGAAAAAAAAAATAAACACCCGTAAGGATTTTG‑TTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                           <  2:315012/145‑1 (MQ=255)                                                                       aaTGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                               >  1:576365/1‑147 (MQ=255)                                                                          ggCACACCCGGCCATTTTTATTGTTAAAAAAACCC‑TTGTTTTGCGTATTTTTATCGAAATAAGGAAAGAAAAAATAACCACCCGCAAGGATTTTT‑TTATAAATGGGACATCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                         <  2:854482/150‑1 (MQ=255)                                                                                 cttcccttttttattaGAAAAAAGACCT‑TTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAACCACCCGAAAGTATTTTG‑TTAAAAATGTGACATTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                         <  1:243929/133‑1 (MQ=255)                                                                                  cTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                >  2:374541/1‑151 (MQ=255)                                                                                   tGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                               >  1:169310/1‑151 (MQ=255)                                                                                     ccgttttttttCGTAAAAAAACCC‑TTCTTTGGCTAACTTTTATGGAAATAAAGAAAAAAAAAATTAACAACGGAAAGATTTTTT‑TTATAAATGGGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcgc                                                              <  2:7974/141‑1 (MQ=255)                                                                                       gCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                            >  1:774025/1‑150 (MQ=255)                                                                                           ttattaaaaaaaaGACCC‑TTCGTTTGCGTACTTTTACCAAAATAAAAAAAAAAAAAATAAACACCCGCCAGGATATTT‑TTACAAATGTGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                        <  1:189865/140‑1 (MQ=255)                                                                                              ttctaaaaaaGACCT‑TTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGCCAGAATTTTT‑TTATAAAGGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                      <  2:584729/142‑1 (MQ=255)                                                                                              tgcgaaaaaaGCACC‑GTCTTTGGCGTTCTTTTAGGGGAAAAAAGGAAAAAAAAATAAACAGCCGTCAGGTTTTTG‑TTAACAAGGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                     <  1:1023223/143‑1 (MQ=255)                                                                                               ttcaaaaaaaaaCCTTTCGTTTGCGTATTTTTATGGAAATAAAGAAAAAAAAAAAAACCACCCATAAGAATATTT‑TTAAAAATGTGACCCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                      <  1:115992/143‑1 (MQ=255)                                                                                                cGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                   >  2:781235/1‑150 (MQ=255)                                                                                                   ataAAGACCC‑TTGTTTTGCGTACGTTTATGGGAATAAAGGAAAAAAAAATAAACACCCGTAGGGATTTTT‑TTATAAATGTGACGATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                  <  2:403433/148‑1 (MQ=255)                                                                                                   aaaaaGACCT‑TTCTTTTGCGTATGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTT‑TTACAAATGTCACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                 <  2:423832/147‑1 (MQ=255)                                                                                                      aaGACCC‑TTCTTTTGCTTATTTTTATCAAAATAAAGAAAAAAAAAATAACCACCCGACAGGATTTTT‑TTATAAATGTGCCAACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                            <  1:930832/151‑1 (MQ=255)                                                                                                       aaaacccttttttggCGTACTTTTATCGAAATAAAGGAAAAAAAAATTAAAACCCGCGGGGATTTTT‑TTATAAATGTGCCTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                            <  2:890377/137‑1 (MQ=255)                                                                                                       aGAACC‑TTCTTTTGCGTTCTTTTATCGAAATAAAGGAA‑AAAAAAATAACCACCCTGCGGATTTTGTTTATCATGTTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                           <  1:949917/151‑1 (MQ=255)                                                                                                         accc‑TTCGTTTGCGTACGTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTGAGGATTTTT‑TAATAAATGTGCCAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                          <  2:550097/148‑1 (MQ=255)                                                                                                           cc‑TTCTTTGGCGTACTTTTATCGGAATAAAGAAAAAAAAAATTAACACCCGCAAGTATTTTT‑TTATAAATGGAACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                       <  2:574810/151‑1 (MQ=255)                                                                                                              ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                      >  1:800000/1‑150 (MQ=255)                                                                                                                  tttGCGTCCTTTTATGGGAATAAAGAAAAAAAAAATAAACACCCGTCAGAATTTTT‑TAATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                  <  2:963365/150‑1 (MQ=255)                                                                                                                   ttGGGTACGTTTAACGAAATAAAGAAAGAAAAAATAACCAACCGAAAGTTTTTTT‑TTAACAATTTAACGTCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                   <  1:232331/148‑1 (MQ=255)                                                                                                                         aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                          >  2:449507/1‑151 (MQ=255)                                                                                                                           ttttattaaaaTAAAGAAAAAAAAAATAAACACCCGTGAGAATATTG‑TTATAAAGGAGCAGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                      <  1:809566/128‑1 (MQ=255)                                                                                                                           ttttAGCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGAATATTT‑TTATAAATGTAACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                         <  1:467846/149‑1 (MQ=255)                                                                                                                           tttaaTCGAAATAAAGAAAAAAAAAATAAACAACGGTAAGGTTATGT‑TTACCAAACTGCCTACCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                           <  1:783416/144‑1 (MQ=255)                                                                                                                              tgtCGAAATAAAGGAAAAAAAAAATAAAAACCTTAAGGATTTTT‑TTAAAAATGTGACAACCCAATCGCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                      <  1:646444/148‑1 (MQ=255)                                                                                                                              tATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGGATTTTG‑TTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                      <  2:315609/150‑1 (MQ=255)                                                                                                                               aTCGGAATAAAAGAAAAAAAAATAAACAACCGGCGGGTTTTTT‑TTATAATGGGGACCTCCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                     <  2:1064045/150‑1 (MQ=255)                                                                                                                                  ggAATTAAGGAAGAAGAAATTAACACCCGTAAGGATATTG‑TTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                   <  2:450575/149‑1 (MQ=255)                                                                                                                                  gaaaTTAAGGAAGAAAAAATAAACACCGGTAAGGATTTTG‑TTATCATTGTAACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                   <  1:210234/147‑1 (MQ=255)                                                                                                                                    aaTAAAGGAAAAAAAAATTAACACCCGTCAGAATTTTT‑TTACAAATGTACCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                   <  2:713085/147‑1 (MQ=255)                                                                                                                                     aTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt               >  1:518311/1‑150 (MQ=255)                                                                                                                                        aaGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc           >  1:1032633/1‑151 (MQ=255)                                                                                                                                        aaGAAAAAAAAAATAAACAACCATAAGAATTTTG‑TTATCAATGTGACCCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                  <  1:210594/144‑1 (MQ=255)                                                                                                                                             agaagaaATTAACAACCGTCCGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt       >  1:1083590/1‑150 (MQ=255)                                                                                                                                                 gaaATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt  >  1:82721/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        CGTTCGGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTT  >  NC_000913/2421254‑2421547

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