Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,379,638	0	G	A	39.0%	81.2 / 21.6	60	F45F (TTC→TTT)	frdC	fumarate reductase (anaerobic), membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base G (27/9);  new base A (0/23);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 2.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCG  >  NC_000913/4379498‑4379788                                                                                                                                                |                                                                                                                                                         ggTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAAGCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAggt                                                                                                                                               <  2:251195/151‑3 (MQ=255)  gTGCCAGTTCAAACCAGGTGTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTc                                                                                                                                              <  2:336453/151‑1 (MQ=255)    gCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCa                                                                                                                                             >  1:362611/1‑150 (MQ=255)        gTTCAATCCAGGTTTTGGTGTGCAGCAGATCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAGCCAGG‑TTTTGTAAAAAGTCGACGAAGCCCGCACA‑GGCTTCCGGGCCATTTTGCAGGGTAAA‑CAGCACG‑AAAATCAGTTCAATg                                                                                                                                          <  2:290910/149‑1 (MQ=255)        gTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCt                                                                                                                                        <  2:1126504/151‑1 (MQ=255)           cAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑a                                                                                                                                     <  1:140988/150‑1 (MQ=255)                 aGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACcacac                                                                                                                              <  1:1221584/151‑1 (MQ=255)                      ttGGTGTGCAGCAGAGCTGTCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTC‑AACCACACAGCCg                                                                                                                         >  2:705613/1‑151 (MQ=255)                        ggTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCgg                                                                                                                        >  2:358574/1‑150 (MQ=255)                         gtgtGCAGAAGAGCTGCCGCCAGAGTTATCAGGTTAATGATCACGATAGCCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCT‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGaa                                                                                                                      <  2:1038952/151‑1 (MQ=255)                              cagcagAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAAccgcc                                                                                                                 >  2:250132/1‑151 (MQ=255)                               agcagAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGCCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAAccgccg                                                                                                                >  1:1152719/1‑151 (MQ=255)                                        gccgccAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTcgc                                                                                                       >  1:319965/1‑151 (MQ=255)                                           gccAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCgcgc                                                                                                     >  1:199540/1‑150 (MQ=255)                                               gagTGATCAGGCTAATGATCCCGACAGCCGCG‑TTTTGTAAAAAGTAGACGACTCCAGCCCA‑GGCTTCCGGTCCATTTTTCAGGGCAAA‑CAGCCCG‑ATAATCAGTTCAATGCTG‑ACCTACACAACCGGAAACGCCGTGCCCTCGCGCAGCAt                                                                                                >  2:246106/1‑151 (MQ=255)                                                  tGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGt                                                                                              >  1:550340/1‑150 (MQ=255)                                                  tGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGt                                                                                              >  1:1262350/1‑150 (MQ=255)                                                     tCAGGTTAATGATCACGATAACGGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCTCACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTaaa                                                                                           >  2:350241/1‑150 (MQ=255)                                                     acAGGTAAATAATCACAATACCGGGT‑TTTTAAAAAAAGTCGACAACCCCCCCCCA‑GGCTCCCGGGCCTTTTTCCGGGGAAAC‑CAGCCCA‑AAAATCATTTCAATGCTA‑ACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTaaaa                                                                                          <  1:300733/150‑1 (MQ=255)                                                      cAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTaaa                                                                                           >  1:849651/1‑149 (MQ=255)                                                        gtttAATGATCACGATAACCGGT‑TTTTGAAAAAAGTCGACGAACCCCGCCCA‑GGCTCCCGGGCCTTTTTTCAGGGCAAA‑CAGCCCA‑AAAATCAGTTAAATGCTG‑AACCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGc                                                                                        <  1:1165505/148‑1 (MQ=255)                                                         ttaaaTGAGCAGGAAAACGGGT‑TTTTGAAAAAGCCCGACAAACCCCCCCCG‑GGCTTCCGGCCCTTTTTTGGGGGCAAA‑CAGCCAA‑AAAACCAGTTCAATGTAA‑ACCCCCCCACCGGGATCCGCCGTGCCTTCGCGCAGCATGTaaaa                                                                                          <  2:1267482/144‑1 (MQ=255)                                                          taaataatCACGATACCGGGT‑TTTTAAAAAAAGTCGACAAACCCCCCCCG‑GGCTCCCGGGCTTTTTTCCGGGGAAAA‑CACCCCA‑AAAATCATTTAAATGCGG‑ACCCACACACCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGa                                                                                      <  1:1066915/148‑1 (MQ=255)                                                            aaTGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATa                                                                                    >  2:1256473/1‑150 (MQ=255)                                                              tGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATaaa                                                                                  >  1:518179/1‑150 (MQ=255)                                                              tGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATaaa                                                                                  >  2:326507/1‑150 (MQ=255)                                                               gATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAc                                                                                >  2:519270/1‑151 (MQ=255)                                                                 accgcaagaCCGGGTTTTTGAAAAAACCCGGGAAACCCCGCCC‑‑GGGTTCGGGGGCACTTTTCCGGGGAAA‑ACGCCCGAAAAATCAGTACAAGATAG‑ACCCACAAAGCCGGACCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAc                                                                                <  2:1280795/142‑1 (MQ=255)                                                                  cACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGc                                                                             >  1:753382/1‑151 (MQ=255)                                                                   aCGATACCGGGT‑TTTTGAAAAAAGTCGACAAACCCCGCCCA‑GGCTCCCGGGCCTTTTTTCAGGCCAAA‑CAGCCCA‑AAAATCAGTTAAATGCTA‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATaaaa                                                                                 <  1:331888/146‑1 (MQ=255)                                                                      aTACCCGGG‑TTTTAAAAAAAGTCGAGGAACCCCGCCCA‑GGCTTCGGGGCCTTTTTTCAGGGAAAA‑CAGCCAA‑AAAATCAGTTCAATTCTA‑AACCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAAt                                                                          <  1:906367/150‑1 (MQ=255)                                                                       tAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAC‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGGGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAAttt                                                                        >  2:1168232/1‑151 (MQ=255)                                                                             gg‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTcc                                                                   >  2:18849/1‑150 (MQ=255)                                                                                 tttGTAAAAAGTCGCCGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGACGGAACAGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTccac                                                                 >  2:822206/1‑149 (MQ=255)                                                                                 tttGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTccacc                                                                >  1:324343/1‑150 (MQ=255)                                                                                  ttGAAAAAGGTCGCCGAACCCCGCCCG‑GGCTTCGGGCCCTTTTTTCAGGTAAAA‑AACCCCA‑AAAATCAGTCCAGGGCTA‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTccacca                                                               <  2:1207548/150‑1 (MQ=255)                                                                                     tAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTg                                                           >  1:1013895/1‑151 (MQ=255)                                                                                      aaaaaGCCGACGACCCCCCCCCA‑GGCTCCCGGGCCTTTTTTCAGGGCAAC‑CAGCCCA‑AAAATCAGTTCAATGCTA‑ACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTgg                                                          <  1:585064/151‑1 (MQ=255)                                                                                        aaaGTCGACGAACCCCCCCCA‑GGCTTCCGGGCCTTTTTTCAGGGCAAA‑CACCCCA‑AAAATCATTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGAc                                                        <  1:769745/151‑1 (MQ=255)                                                                                          aggccgccaacccccccCAGGGCTCCCGGGCCTTTTTTCGGGGAAAA‑CCGCCCA‑AAAATCAGTTAAATGCAA‑AACCACACACCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACg                                                       <  2:550308/140‑1 (MQ=255)                                                                                          aGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGt                                                      >  1:436653/1‑151 (MQ=255)                                                                                              gacgaATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCAG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGCCGTCAtt                                                  >  2:799068/1‑151 (MQ=255)                                                                                              acccactCCCGCCCA‑GCCTTCGGGGCCTTTTTCCGGGCCAAA‑GGGCCGA‑AAAATCATTTAAATATAG‑ACCCCCCCACCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCCCCAGGTGGACGTCAt                                                   <  2:175629/144‑1 (MQ=255)                                                                                                cgaATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTgg                                                >  2:398229/1‑151 (MQ=255)                                                                                                caaaccccGCCCA‑GGCTTCCGGCCTATTTTACAGGGAAAA‑CAGCCCA‑AAAATTATTGAAATGCAA‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCAtt                                                  <  1:744080/144‑1 (MQ=255)                                                                                                     cccGCCCT‑GG‑TTTCGGGACTTTTTGCAGGGGAAACCAGCCAA‑AAATGTAGTTCAATACTA‑AACCACACGGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTa                                           <  2:598321/151‑1 (MQ=255)                                                                                                         cccA‑GGCTTCCGGGCCTTTTTTCAGGGCAAA‑CACCCAA‑AAAATCAGTTCAATGCTG‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACAt                                        <  1:67138/150‑1 (MQ=255)                                                                                                                 ttCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGttt                                 >  1:701956/1‑150 (MQ=255)                                                                                                                 tgcggggCCATTTTT‑AGGGCAAACCAGCCCA‑AAATTCAGTAAAATGTGG‑AACCACAAGGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCTCCAGGTGGACGTCATTGGCCTTACATACGGt                                   <  2:1063193/144‑1 (MQ=255)                                                                                                                   ccGGCCCTTTTTTCAGGGCAAA‑CACCCGA‑AAAATCAGTTAAATGCTA‑ACCCACACACCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTa                                <  1:560320/149‑1 (MQ=255)                                                                                                                     gggCCTTTTTTCAGGGCAAA‑CACCCAA‑AAAATCAGTTAAATGCTG‑ACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACg                              <  1:86632/149‑1 (MQ=255)                                                                                                                        ccATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGttt                           >  2:822234/1‑149 (MQ=255)                                                                                                                         cATTTTTCCGGGGAAACCCCCCCA‑AAAAGCAGTAAAATCCTA‑ACCCCCACAGCCGGACCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAg                         <  2:1217765/151‑1 (MQ=255)                                                                                                                           tttttCGGGGCAAC‑CACCCCA‑AAAATCAGTTAAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgt                        <  2:393250/149‑1 (MQ=255)                                                                                                                           tttttCAGGGCAAA‑CCCCCCG‑AAAATCAGTTCAATTCTA‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgt                        <  2:1207444/149‑1 (MQ=255)                                                                                                                           tttttCAGGCAAAC‑CCGCCCA‑AAATTCAGTTCAATGCTA‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGttt                           <  2:1136385/146‑1 (MQ=255)                                                                                                                                  gggCAAA‑CACCCCA‑AAAATCATTTCAATGCTG‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATg                 <  1:1266786/149‑1 (MQ=255)                                                                                                                                  gggAAAA‑CACCCAA‑AAAATCAGTTCAATGCTG‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCAt                  <  1:53462/148‑1 (MQ=255)                                                                                                                                       aaccAGCCAA‑AAAATCTTTTAAATGCTA‑ACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATg                 <  1:782344/142‑1 (MQ=255)                                                                                                                                          cAGCCCG‑AAAA‑CATTTTAAAGGTGCACCCACACACCAGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACt          <  2:168616/149‑1 (MQ=255)                                                                                                                                                g‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCg  >  1:409238/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                         GGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGG‑TTTTGTAAAAAGTCGACGAATCCCGCCCA‑GGCTTCCGGGCCATTTTTCAGGGCAAA‑CAGCCCG‑AAAATCAGTTCAATGCTG‑AACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCG  >  NC_000913/4379498‑4379788

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