Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,741,500 T→C 13.1% E28G (GAG→GGG)  ysaA ← putative hydrogenase, 4Fe‑4S ferredoxin‑type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,741,5000TC13.1% 150.5 / 2.3 68E28G (GAG→GGG) ysaAputative hydrogenase, 4Fe‑4S ferredoxin‑type component
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (3/6);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAGGTCATCCCTGACA  >  NC_000913/3741273‑3741714
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         
tCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATtctc                                                                                                                                                                                                                          >  1:590177/1‑229 (MQ=255)
 cGCACCAAACGGGCAAGCCAGCATTCAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATAGGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCa                                                                                                                                                                                                                         <  1:348020/229‑1 (MQ=255)
 cGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCCATGCAACGTGTTTGTTCAACCAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCa                                                                                                                                                                                                                         <  1:242312/229‑1 (MQ=255)
             gCAAGCCAGCATAAAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGGCTTTAATGACACGAATACGGGAAATAAAATCGTCTGGTGGCAACGCAGCGCAATCCTGATTCACAGAATGCGaaca                                                                                                                                                                                                              >  1:609373/1‑225 (MQ=255)
                aGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGGCTTTCATGACACGAATACGGGAAATAAACCCGCCTGGGGACACACCAGCGCAATCCTCATCCCCATAATGGCACA‑CT‑gcc                                                                                                                                                                                                        >  2:169592/1‑228 (MQ=255)
                  gcagactgacgCTTTTTCAGCCAATGCCACGTGTTTGTTCAACGAAAATTGGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGCCATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCa                                                                                                                                                                                                      <  1:705668/219‑1 (MQ=255)
                           aGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGACTTTAATGACACGAATACGGGAAATAAACTCGCCTGGTGACAACCAAGCGCAAACCTGATTCCCATGATGCGACA‑CC‑GCGCAAGCCCCtca                                                                                                                                                                                             >  2:692014/1‑227 (MQ=255)
                             cTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGGGGTCTTTAATGACACGAATACGGGCAATAAACTCGTCTGGTGACAACGCAACGCAACACCGATTCTCCTGAATAGACA‑CT‑GCGCAAGCCACTTaac                                                                                                                                                                                           >  2:171491/1‑229 (MQ=255)
                                    aGCCAATGAAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCCCGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACg                                                                                                                                                                                    <  1:388298/229‑1 (MQ=255)
                                       cAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGCCACGCATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGAAATCCTGATTCTCATGATGCTCAA‑CT‑GGGCAAGCAACTTCACAGGTCCGCCa                                                                                                                                                                                 >  1:297117/1‑229 (MQ=255)
                                          tGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGTGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGGGACACCGCAGCGCAATCCTGATTCTCATGATACGACA‑CT‑ACGCAAGCCACTTCCAAGGACCGGCAAc                                                                                                                                                                               >  1:739808/1‑228 (MQ=255)
                                             aaCGTGTTCTGTCAACTAAAATACGCCCATGTCCCCGGTTTAACGCGTAAACAGGGTAAACCATAGCCCATGGTGCATTTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATGCGGTAAATAAACTCTTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCTCAAGCCACTTCACAGGTACGGCAACCGAt                                                                                                                                                                           <  1:732439/229‑1 (MQ=255)
                                                gtgtTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GGGCAAGCCACTTCACAGGTACGGCACCCGATA‑c                                                                                                                                                                        >  1:446140/1‑228 (MQ=255)
                                                     tgttCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACCCTGATGACAGGCTACGGCCGAGGTCCAGCAGTGGTCTTTAATGACCCGAATCCGGGAAATCAACTCGTCAGTTGACACCGCAGCCAAATCCTTATTCGACAGATGCGACC‑CT‑GCGCCAGCCCATTCACAGGAACAGCAAACGAGA‑AATTtcg                                                                                                                                                                  >  2:695101/1‑229 (MQ=255)
                                                           aCGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTTATTCTCATGATGCGACA‑CT‑CGGCAAGCCCCTTCCCAAGAAACGCAACCAGAA‑CATTtcgtccccc                                                                                                                                                            >  2:724095/1‑225 (MQ=255)
                                                           aCGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGCAC‑CT‑GCGCAACCCACTTCCCAGGTAGGGCAACCAATC‑AATTTCGCCGCac                                                                                                                                                            >  2:687264/1‑228 (MQ=255)
                                                            cGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCCCATGATGCGAAC‑CT‑GCCCAAGCAACTTCCCAGGGACGCAAACCGAAA‑CATTtcatcgacac                                                                                                                                                           >  2:494179/1‑225 (MQ=255)
                                                              aaaaTATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCCCGGTGCATCTTCAAACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCg                                                                                                                                                         <  1:363972/229‑1 (MQ=255)
                                                                         aTGTTCGTGGCTTAACGGGTACACAGGGACGGAATTTGCGCACGGTGCAGCTTACCACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAAc                                                                                                                                              <  1:153608/229‑1 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGGGTCAACAGGGCAGACATTCGCGCCCGGTCCATTTTCACACTGATGACAGGCTACTGCCGGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAAc                                                                                                                                              <  1:171491/229‑1 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGTTGCGaca                                                                                                                                                                                                               >  2:436444/1‑167 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGaca                                                                                                                                                                                                               <  1:436444/167‑1 (MQ=255)
                                                                          tGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCCCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCACTAATAAAc                                                                                                                                              >  1:48619/1‑228 (MQ=255)
                                                                              cgccgcTTAGCGCGTCAACAGGGCAGACTTTCGCGCACGGTGCATATCCACACGTATGACAGGCAACTGCCTGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGTAAATAAACGCGTCTGGTGACAACGCAGCGCAATCCTTATTCTCAGGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGTTA‑CATTTCGTCGCATCCGCAATAATAAACCGGtt                                                                                                                                         <  1:111852/229‑1 (MQ=255)
                                                                               gcgGCTTATCGCCTCAACAGGGCAGCAATTCGCGCACGGTGCATCTTCACACTGATGACGGGCGACGGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACC‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTc                                                                                                                                        <  2:446140/229‑1 (MQ=255)
                                                                               gcagtctATCGCGTCATACGGGCATCCATTGCCGCACGGTGCATCTTCACTCTAGTGACAGGCTCCTTCCTTGGTCCATCAGTGTTCTTTAATGACACGAATACGGGAAATAACCTCTTCTGGTAACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTc                                                                                                                                        <  1:722567/223‑1 (MQ=255)
                                                                                      aTCGCGTCAACAGGTCAGACATTCGCGCCCGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGCAATAAACTCGTCTGGGGACAACGCACCGCAATCCTGATTCTCATGATGCTCCA‑CT‑GCGGAGGCCCCTTACCAGCAACGGCAACC‑ATACCATTACTTACCATCCGCAATAAAAAACCGTGTCATAAtac                                                                                                                                 >  2:443573/1‑227 (MQ=255)
                                                                                             cAACAGGGCAGACATTCGCGCCAGGTGCATCTTAACACTGATGACAGGCTAATGCGGTGGTCCAGCAGTGGTCTTTACTGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCtt                                                                                                                           <  1:615023/228‑1 (MQ=255)
                                                                                                 agggcgagcagtCGCGCACGGGGCATCTGCACATGTATGACAGGCTACTGCCGTGGTCAAGTAGTGGTCTTCAATGACACGAACACGGGAAATAAACTCTTCTGGTGAAAACCCCGCCCAATCCTGGTTCCCATGATGCGACA‑CT‑GCGCAAGCCCCTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  2:69936/218‑1 (MQ=255)
                                                                                                 aGGGCAGCCATTTGCGCAAGGTGCATCTTCCCGCGAAAGACGGGCTACGGCCGTGGTCCAGAAGTGGTCTTAAATGCCAACAATACGGGAAATAAACTCGTCTGGTGACACCGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  1:351550/229‑1 (MQ=255)
                                                                                                 aGGGCAGACATTCGGCCAGGGTGAATCTTCACACGGATGACAGGCTACGGCCGTGGCCCAGCAGGGGTTTGTAATGACCCGAATCCGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  2:590176/229‑1 (MQ=255)
                                                                                                             cgcgCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCaga                                                                                                          >  1:172783/1‑229 (MQ=255)
                                                                                                             cgcgCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGAGACA‑CT‑GCGCAAGCCACTTCACAGGAACGGCAACCGCTC‑AATTTCGTCGCATCCGCAAAAATAAACCGGTTTTTACTCCCTCCTTTTCCCTAAAacaacaaa                                                                                                          >  2:594097/1‑227 (MQ=255)
                                                                                                              gcgcgcGGTGCATCTTCACACTGATGGCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAAGGCCACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:687265/226‑1 (MQ=255)
                                                                                                              gcgcACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCCCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:386403/226‑1 (MQ=255)
                                                                                                              ccgccaggTGCTTTTGCACACGGAAGACAGGCTAATGCCGTGGTCCGGCAGGGGTCTTTAATGACAAGAATCCGGGAAATAAACTCGTCGGGTGACAACGCAGCGAAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCaga                                                                                                          <  2:428728/222‑1 (MQ=255)
                                                                                                                     gTGCATCTTCACACTGAGGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCaga                                                                                                  <  1:518466/229‑1 (MQ=255)
                                                                                                                      tGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCat                                                                                                   >  1:395250/1‑226 (MQ=255)
                                                                                                                        cATCTTAACCCGGATGAAAAGGTACTGCCGTGTCCCAGCAGTGGTATTTAATTACACGATGACGGGAAATAAACTCGTCTGGTGACAACCCAGCGCAATCTTGATTACCATTATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCagagag                                                                                               <  2:33490/229‑1 (MQ=255)
                                                                                                                        cATATTCACACTGATGACAGGCTAACGCCGTGGTCCAGCAGTGGTCTTTAATGCCACGAATTCGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACC‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTTATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCagagag                                                                                               <  1:594098/229‑1 (MQ=255)
                                                                                                                            ttCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATCATGCGACA‑CT‑GCGCAAGCCACTTCCCACGGACGGCAACCGATA‑CCTTTCGTCGCATCCGCAATACTAAACCGGTTCCTGATCCCCCCTTTGCACTCAAACCCCAAAGTTCCaagaaacgct                                                                                           >  2:680430/1‑220 (MQ=255)
                                                                                                                             tCACTCTGATGACCGGCTATGGCCTTGGTCCAGCATGGGTCTTTAATGACACGAATACGGGAAATAACCTCGTCTGGTGACAACGCAGCGCAATCCGGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGt                                                                                          <  1:169592/229‑1 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCCATAATAAACCGGTTCATGATCACTCCTATGCACTAAAACAGCAGACGGCCacagaaccgctcca                                                                                       >  2:728967/1‑220 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTccc                                                                                       >  2:172791/1‑229 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTcc                                                                                        >  1:194954/1‑228 (MQ=255)
                                                                                                                                     atgaCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACGGCAGAGTGCCAGAGAACGCGTCCCTAATg                                                                                  >  1:311610/1‑229 (MQ=255)
                                                                                                                                     atgaCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAAAGTGCCAGAGAGCGCGTCCCTACTg                                                                                  >  1:364068/1‑229 (MQ=255)
                                                                                                                                           ggCTACTGCCGGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGa                                                                            <  1:443573/229‑1 (MQ=255)
                                                                                                                                                      tGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            >  2:499373/1‑186 (MQ=255)
                                                                                                                                                      tGGTACAGCAGTGGTTTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:499373/186‑1 (MQ=255)
                                                                                                                                                       ggTCCCGCAGGTGTTTTTAAAGAATCGAATACGGGAAATAAACCCCGCTGGTCACAACCGAGCGCAATCCTGATTCACATGATGCGACA‑CG‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTATGCACTAAAACAGGAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAgg                                                                 <  1:728968/228‑1 (MQ=255)
                                                                                                                                                         tCCAGCAGTGGTCTTTAATGACACGAATCCGGAAAATAACCCCTGTTGGGGACACCGCACCGAAATCAGAATTCCCATGA‑GCGACACCT‑CCGCAAGCCCCTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGaa                                                               <  2:194954/228‑1 (MQ=255)
                                                                                                                                                              cagTAGTCTTTAATGACACGAATACGGGAAATAACCTCGTCTGGTGACAACGGAGCCCAATCCTGAGTCCCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGtt                                                          <  1:347044/228‑1 (MQ=255)
                                                                                                                                                                   gccTTTAATGACACGAGTACGGGAAATAAACTCGTCTGGTGACAACGCAGAGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCCACCGATA‑CATTTCGTCGCATCCCCAATAATAAACCGGTTCATGATCACTCCTTTGAACTAAACCAGCAACGTGCCACAGAGCGCGTCCCTACTGCTTTAGTCCACCGATAGAAATGTTAAAGtt                                                    >  1:754343/3‑229 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGt                                                     >  2:754509/1‑228 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACGGCTTGGATCAAACAAAGGAAATGTTAAAGt                                                     >  2:672510/1‑228 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCAGGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTCCTTTGTTCAAACGAAGGAAATGTTa                                                         >  1:730822/1‑224 (MQ=255)
                                                                                                                                                                             acacGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTg                                          >  1:345494/1‑229 (MQ=255)
                                                                                                                                                                                 gAAGACGGGAAATAACCTCGCCTGGTGACAACGCAGCCCAATCCTGATTCTCATGATGCGACA‑CT‑GCCCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCa                                      <  1:724096/229‑1 (MQ=255)
                                                                                                                                                                                      cGGAAAAAAAACTCTGTTGGTAACAACCCCGCGCATGCCTGATTCTAATGAGGCCACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCg                                 <  2:364068/229‑1 (MQ=255)
                                                                                                                                                                                         gAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAAGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGtc                               >  1:727307/1‑227 (MQ=255)
                                                                                                                                                                                               aaCTCGTCTGGTTACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGAt                        <  1:172791/229‑1 (MQ=255)
                                                                                                                                                                                                aCTCTTCGGGGGACAACGCAGCGCAATCCTGATTCCCATGATGCGAAA‑CG‑GCGAAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATc                       <  1:360856/229‑1 (MQ=255)
                                                                                                                                                                                                     tCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            >  2:502306/1‑139 (MQ=255)
                                                                                                                                                                                                     tCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:502306/139‑1 (MQ=255)
                                                                                                                                                                                                         gtcaCAACGCACCGCAACCCTGATTCTCATTAAGCGA‑A‑CTGGCCCAAGCCACTTCACAGGTTCGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTTCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAgg              <  1:692015/226‑1 (MQ=255)
                                                                                                                                                                                                                    gcgcAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACAGCTTTGTTCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGACCAAAGGGTTATGTGATCGCGCCCAGGACATCCCTTAc   >  1:688644/1‑229 (MQ=255)
                                                                                                                                                                                                                    gcgcAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCCCGTCCCTACTGCTTTGATCCAACGAAGGAAAGGTTAAAGTTTTGGCCCTTTTCCAAAGCGTTATGTAATCGCGACCAGGCAATcccgcac   >  2:338806/1‑225 (MQ=255)
                                                                                                                                                                                                                     cgcAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAACGGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTAATCGCGCCCAGGGCATCCCTGACa  >  1:608738/1‑229 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         
TCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACA‑CT‑GCGCAAGCCACTTCACAGGTACGGCAACCGATA‑CATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAGGTCATCCCTGACA  >  NC_000913/3741273‑3741714

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