Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	44.9%	100.7 / 28.0	69	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (32/6);  new base C (0/31);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 5.89e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421295‑2421592                                                                                                                                                       |                                                                                                                                                      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)                                                                                                    aaaaaaaaaaaacaacCGGCAGATTTTTT‑TTATAAATGGAACGACCCCCCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  1:640954/140‑1 (MQ=255)                                                                                                        gaaATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               >  1:82721/1‑151 (MQ=255)                                                                                                         aaaTTAACACCCGTATGATTTTTT‑TTATCAATGTGACCATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                              <  1:1082744/151‑1 (MQ=255)                                                                                                            acaacaccCTTCAGGTTATTG‑TATTAAATGTACCCACCACCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                            <  1:1063620/148‑1 (MQ=255)                                                                                                                  acCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      >  2:538877/1‑150 (MQ=255)                                                                                                                   ccGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     >  2:1019791/1‑150 (MQ=255)                                                                                                                      tCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  >  1:651329/1‑150 (MQ=255)                                                                                                                        aGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                               >  1:434685/1‑151 (MQ=255)                                                                                                                           atatTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  2:723815/1‑150 (MQ=255)                                                                                                                            tatTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  1:789094/1‑149 (MQ=255)                                                                                                                             tttttttATAAATGTAACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                                <  1:265362/141‑1 (MQ=255)                                                                                                                                   tATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     >  2:150224/1‑151 (MQ=255)                                                                                                                                       aaTGTGACGCCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 <  1:46421/151‑1 (MQ=255)                                                                                                                                           tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat             >  1:1065661/1‑151 (MQ=255)                                                                                                                                           tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat             >  1:775137/1‑151 (MQ=255)                                                                                                                                                  ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          <  1:798711/147‑1 (MQ=255)                                                                                                                                                   cACTACCCCGCCGATAAAGCGCGCACTTGCGAAGATATCACAGATGCGCTGAATCATCGCACCGTAACGAACAACATTCTTCAGCATGTCGGGAATAACCGTTTCTCTTGGCTGGAATAATTAACTCAGGATGCGCCCGACACCGCCCGt      >  1:471801/1‑150 (MQ=255)                                                                                                                                                   cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:460755/1‑151 (MQ=255)                                                                                                                                                    ccTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  1:366651/2‑150 (MQ=255)                                                                                                                                                      tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:731879/1‑151 (MQ=255)                                                                                                                                                      tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:212403/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421295‑2421592

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: ‑