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

TCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAGGTCATCCCTGACA  >  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)
                aGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGGCTTTCATGACACGAATACGGGAAATAAACCCGCCTGGGGACACACCAGCGCAATCCTCATCCCCATAATGGCACACTgcc                                                                                                                                                                                                        >  2:169592/1‑228 (MQ=255)
                  gcagactgacgCTTTTTCAGCCAATGCCACGTGTTTGTTCAACGAAAATTGGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGCCATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCa                                                                                                                                                                                                      <  1:705668/219‑1 (MQ=255)
                           aGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGACTTTAATGACACGAATACGGGAAATAAACTCGCCTGGTGACAACCAAGCGCAAACCTGATTCCCATGATGCGACACCGCGCAAGCCCCtca                                                                                                                                                                                             >  2:692014/1‑227 (MQ=255)
                             cTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGGGGTCTTTAATGACACGAATACGGGCAATAAACTCGTCTGGTGACAACGCAACGCAACACCGATTCTCCTGAATAGACACTGCGCAAGCCACTTaac                                                                                                                                                                                           >  2:171491/1‑229 (MQ=255)
                                    aGCCAATGAAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCCCGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACg                                                                                                                                                                                    <  1:388298/229‑1 (MQ=255)
                                       cAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGCCACGCATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGAAATCCTGATTCTCATGATGCTCAACTGGGCAAGCAACTTCACAGGTCCGCCa                                                                                                                                                                                 >  1:297117/1‑229 (MQ=255)
                                          tGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGTGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGGGACACCGCAGCGCAATCCTGATTCTCATGATACGACACTACGCAAGCCACTTCCAAGGACCGGCAAc                                                                                                                                                                               >  1:739808/1‑228 (MQ=255)
                                             aaCGTGTTCTGTCAACTAAAATACGCCCATGTCCCCGGTTTAACGCGTAAACAGGGTAAACCATAGCCCATGGTGCATTTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATGCGGTAAATAAACTCTTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCTCAAGCCACTTCACAGGTACGGCAACCGAt                                                                                                                                                                           <  1:732439/229‑1 (MQ=255)
                                                gtgtTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGGGCAAGCCACTTCACAGGTACGGCACCCGATAc                                                                                                                                                                        >  1:446140/1‑228 (MQ=255)
                                                     tgttCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACCCTGATGACAGGCTACGGCCGAGGTCCAGCAGTGGTCTTTAATGACCCGAATCCGGGAAATCAACTCGTCAGTTGACACCGCAGCCAAATCCTTATTCGACAGATGCGACCCTGCGCCAGCCCATTCACAGGAACAGCAAACGAGAAATTtcg                                                                                                                                                                  >  2:695101/1‑229 (MQ=255)
                                                           aCGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTTATTCTCATGATGCGACACTCGGCAAGCCCCTTCCCAAGAAACGCAACCAGAACATTtcgtccccc                                                                                                                                                            >  2:724095/1‑225 (MQ=255)
                                                           aCGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGCACCTGCGCAACCCACTTCCCAGGTAGGGCAACCAATCAATTTCGCCGCac                                                                                                                                                            >  2:687264/1‑228 (MQ=255)
                                                            cGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCCCATGATGCGAACCTGCCCAAGCAACTTCCCAGGGACGCAAACCGAAACATTtcatcgacac                                                                                                                                                           >  2:494179/1‑225 (MQ=255)
                                                              aaaaTATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCCCGGTGCATCTTCAAACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCg                                                                                                                                                         <  1:363972/229‑1 (MQ=255)
                                                                         aTGTTCGTGGCTTAACGGGTACACAGGGACGGAATTTGCGCACGGTGCAGCTTACCACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAAc                                                                                                                                              <  1:153608/229‑1 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGGGTCAACAGGGCAGACATTCGCGCCCGGTCCATTTTCACACTGATGACAGGCTACTGCCGGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAAc                                                                                                                                              <  1:171491/229‑1 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGTTGCGaca                                                                                                                                                                                                               >  2:436444/1‑167 (MQ=255)
                                                                         aTGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGaca                                                                                                                                                                                                               <  1:436444/167‑1 (MQ=255)
                                                                          tGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCCCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCACTAATAAAc                                                                                                                                              >  1:48619/1‑228 (MQ=255)
                                                                              cgccgcTTAGCGCGTCAACAGGGCAGACTTTCGCGCACGGTGCATATCCACACGTATGACAGGCAACTGCCTGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGTAAATAAACGCGTCTGGTGACAACGCAGCGCAATCCTTATTCTCAGGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGTTACATTTCGTCGCATCCGCAATAATAAACCGGtt                                                                                                                                         <  1:111852/229‑1 (MQ=255)
                                                                               gcgGCTTATCGCCTCAACAGGGCAGCAATTCGCGCACGGTGCATCTTCACACTGATGACGGGCGACGGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACCCTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTc                                                                                                                                        <  2:446140/229‑1 (MQ=255)
                                                                               gcagtctATCGCGTCATACGGGCATCCATTGCCGCACGGTGCATCTTCACTCTAGTGACAGGCTCCTTCCTTGGTCCATCAGTGTTCTTTAATGACACGAATACGGGAAATAACCTCTTCTGGTAACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTc                                                                                                                                        <  1:722567/223‑1 (MQ=255)
                                                                                      aTCGCGTCAACAGGTCAGACATTCGCGCCCGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGCAATAAACTCGTCTGGGGACAACGCACCGCAATCCTGATTCTCATGATGCTCCACTGCGGAGGCCCCTTACCAGCAACGGCAACCATACCATTACTTACCATCCGCAATAAAAAACCGTGTCATAAtac                                                                                                                                 >  2:443573/1‑227 (MQ=255)
                                                                                             cAACAGGGCAGACATTCGCGCCAGGTGCATCTTAACACTGATGACAGGCTAATGCGGTGGTCCAGCAGTGGTCTTTACTGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCtt                                                                                                                           <  1:615023/228‑1 (MQ=255)
                                                                                                 agggcgagcagtCGCGCACGGGGCATCTGCACATGTATGACAGGCTACTGCCGTGGTCAAGTAGTGGTCTTCAATGACACGAACACGGGAAATAAACTCTTCTGGTGAAAACCCCGCCCAATCCTGGTTCCCATGATGCGACACTGCGCAAGCCCCTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  2:69936/218‑1 (MQ=255)
                                                                                                 aGGGCAGCCATTTGCGCAAGGTGCATCTTCCCGCGAAAGACGGGCTACGGCCGTGGTCCAGAAGTGGTCTTAAATGCCAACAATACGGGAAATAAACTCGTCTGGTGACACCGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  1:351550/229‑1 (MQ=255)
                                                                                                 aGGGCAGACATTCGGCCAGGGTGAATCTTCACACGGATGACAGGCTACGGCCGTGGCCCAGCAGGGGTTTGTAATGACCCGAATCCGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCAc                                                                                                                      <  2:590176/229‑1 (MQ=255)
                                                                                                             cgcgCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCaga                                                                                                          >  1:172783/1‑229 (MQ=255)
                                                                                                             cgcgCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGAGACACTGCGCAAGCCACTTCACAGGAACGGCAACCGCTCAATTTCGTCGCATCCGCAAAAATAAACCGGTTTTTACTCCCTCCTTTTCCCTAAAacaacaaa                                                                                                          >  2:594097/1‑227 (MQ=255)
                                                                                                              gcgcgcGGTGCATCTTCACACTGATGGCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAAGGCCACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:687265/226‑1 (MQ=255)
                                                                                                              gcgcACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCCCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:386403/226‑1 (MQ=255)
                                                                                                              ccgccaggTGCTTTTGCACACGGAAGACAGGCTAATGCCGTGGTCCGGCAGGGGTCTTTAATGACAAGAATCCGGGAAATAAACTCGTCGGGTGACAACGCAGCGAAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCaga                                                                                                          <  2:428728/222‑1 (MQ=255)
                                                                                                                     gTGCATCTTCACACTGAGGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCaga                                                                                                  <  1:518466/229‑1 (MQ=255)
                                                                                                                      tGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCat                                                                                                   >  1:395250/1‑226 (MQ=255)
                                                                                                                        cATCTTAACCCGGATGAAAAGGTACTGCCGTGTCCCAGCAGTGGTATTTAATTACACGATGACGGGAAATAAACTCGTCTGGTGACAACCCAGCGCAATCTTGATTACCATTATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCagagag                                                                                               <  2:33490/229‑1 (MQ=255)
                                                                                                                        cATATTCACACTGATGACAGGCTAACGCCGTGGTCCAGCAGTGGTCTTTAATGCCACGAATTCGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACCCTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTTATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCagagag                                                                                               <  1:594098/229‑1 (MQ=255)
                                                                                                                            ttCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATCATGCGACACTGCGCAAGCCACTTCCCACGGACGGCAACCGATACCTTTCGTCGCATCCGCAATACTAAACCGGTTCCTGATCCCCCCTTTGCACTCAAACCCCAAAGTTCCaagaaacgct                                                                                           >  2:680430/1‑220 (MQ=255)
                                                                                                                             tCACTCTGATGACCGGCTATGGCCTTGGTCCAGCATGGGTCTTTAATGACACGAATACGGGAAATAACCTCGTCTGGTGACAACGCAGCGCAATCCGGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGt                                                                                          <  1:169592/229‑1 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCCATAATAAACCGGTTCATGATCACTCCTATGCACTAAAACAGCAGACGGCCacagaaccgctcca                                                                                       >  2:728967/1‑220 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTccc                                                                                       >  2:172791/1‑229 (MQ=255)
                                                                                                                                cacTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTcc                                                                                        >  1:194954/1‑228 (MQ=255)
                                                                                                                                     atgaCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACGGCAGAGTGCCAGAGAACGCGTCCCTAATg                                                                                  >  1:311610/1‑229 (MQ=255)
                                                                                                                                     atgaCAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAAAGTGCCAGAGAGCGCGTCCCTACTg                                                                                  >  1:364068/1‑229 (MQ=255)
                                                                                                                                           ggCTACTGCCGGGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGa                                                                            <  1:443573/229‑1 (MQ=255)
                                                                                                                                                      tGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            >  2:499373/1‑186 (MQ=255)
                                                                                                                                                      tGGTACAGCAGTGGTTTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:499373/186‑1 (MQ=255)
                                                                                                                                                       ggTCCCGCAGGTGTTTTTAAAGAATCGAATACGGGAAATAAACCCCGCTGGTCACAACCGAGCGCAATCCTGATTCACATGATGCGACACGGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTATGCACTAAAACAGGAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAgg                                                                 <  1:728968/228‑1 (MQ=255)
                                                                                                                                                         tCCAGCAGTGGTCTTTAATGACACGAATCCGGAAAATAACCCCTGTTGGGGACACCGCACCGAAATCAGAATTCCCATGAGCGACACCTCCGCAAGCCCCTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGaa                                                               <  2:194954/228‑1 (MQ=255)
                                                                                                                                                              cagTAGTCTTTAATGACACGAATACGGGAAATAACCTCGTCTGGTGACAACGGAGCCCAATCCTGAGTCCCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGtt                                                          <  1:347044/228‑1 (MQ=255)
                                                                                                                                                                   gccTTTAATGACACGAGTACGGGAAATAAACTCGTCTGGTGACAACGCAGAGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCCACCGATACATTTCGTCGCATCCCCAATAATAAACCGGTTCATGATCACTCCTTTGAACTAAACCAGCAACGTGCCACAGAGCGCGTCCCTACTGCTTTAGTCCACCGATAGAAATGTTAAAGtt                                                    >  1:754343/3‑229 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGt                                                     >  2:754509/1‑228 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACGGCTTGGATCAAACAAAGGAAATGTTAAAGt                                                     >  2:672510/1‑228 (MQ=255)
                                                                                                                                                                   gTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCAGGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTCCTTTGTTCAAACGAAGGAAATGTTa                                                         >  1:730822/1‑224 (MQ=255)
                                                                                                                                                                             acacGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTg                                          >  1:345494/1‑229 (MQ=255)
                                                                                                                                                                                 gAAGACGGGAAATAACCTCGCCTGGTGACAACGCAGCCCAATCCTGATTCTCATGATGCGACACTGCCCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCa                                      <  1:724096/229‑1 (MQ=255)
                                                                                                                                                                                      cGGAAAAAAAACTCTGTTGGTAACAACCCCGCGCATGCCTGATTCTAATGAGGCCACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCg                                 <  2:364068/229‑1 (MQ=255)
                                                                                                                                                                                         gAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAAGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGtc                               >  1:727307/1‑227 (MQ=255)
                                                                                                                                                                                               aaCTCGTCTGGTTACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGAt                        <  1:172791/229‑1 (MQ=255)
                                                                                                                                                                                                aCTCTTCGGGGGACAACGCAGCGCAATCCTGATTCCCATGATGCGAAACGGCGAAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATc                       <  1:360856/229‑1 (MQ=255)
                                                                                                                                                                                                     tCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            >  2:502306/1‑139 (MQ=255)
                                                                                                                                                                                                     tCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAAcagca                                                                                                            <  1:502306/139‑1 (MQ=255)
                                                                                                                                                                                                         gtcaCAACGCACCGCAACCCTGATTCTCATTAAGCGAACTGGCCCAAGCCACTTCACAGGTTCGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTTCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAgg              <  1:692015/226‑1 (MQ=255)
                                                                                                                                                                                                                    gcgcAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACAGCTTTGTTCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGACCAAAGGGTTATGTGATCGCGCCCAGGACATCCCTTAc   >  1:688644/1‑229 (MQ=255)
                                                                                                                                                                                                                    gcgcAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCCCGTCCCTACTGCTTTGATCCAACGAAGGAAAGGTTAAAGTTTTGGCCCTTTTCCAAAGCGTTATGTAATCGCGACCAGGCAATcccgcac   >  2:338806/1‑225 (MQ=255)
                                                                                                                                                                                                                     cgcAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAACGGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTAATCGCGCCCAGGGCATCCCTGACa  >  1:608738/1‑229 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         
TCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCTGGTGACAACGCAGCGCAATCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTTCATGATCACTCCTTTGCACTAAAACAGCAGAGTGCCAGAGAGCGCGTCCCTACTGCTTTGATCCAACGAAGGAAATGTTAAAGTTTTGGCACTTGTCCAAAGCGTTATGTGATCGCGCCCAGGTCATCCCTGACA  >  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: