Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,437	0	T	C	41.5%	99.6 / 17.0	65	I38T (ATC→ACC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (28/10);  new base C (0/27);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 3.52e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATAT  >  NC_000913/2421291‑2421581                                                                                                                                                     |                                                                                                                                                 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)                                                                                                        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)                                                                                                                                                     |                                                                                                                                                 ATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTG‑TTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATAT  >  NC_000913/2421291‑2421581

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