Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,561,972 A→T 7.1% intergenic (‑149/‑41) glpE ← / → glpD thiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,561,9720AT 7.1% 300.1 / 11.9 84intergenic (‑149/‑41)glpE/glpDthiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (39/39);  new base T (3/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-02

TCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCC  >  NC_000913/3561831‑3562096
                                                                                                                                             |                                                                                                                            
tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                           >  4:624068/1‑145 (MQ=255)
tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                           <  3:624068/145‑1 (MQ=255)
    tctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                    <  4:517172/148‑1 (MQ=255)
       ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                           <  2:1080526/138‑1 (MQ=255)
       ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                           >  1:1080526/1‑138 (MQ=255)
                aGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACg                                                                                                        >  2:570167/1‑148 (MQ=255)
                  tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     >  1:1027280/1‑133 (MQ=255)
                  tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     <  2:1027280/133‑1 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                        <  3:771178/126‑1 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                        >  4:771178/1‑126 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     >  3:938994/1‑129 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     <  4:938994/129‑1 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                      >  4:469169/1‑128 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                      <  3:469169/128‑1 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTg                                                                                                  >  4:160527/1‑148 (MQ=255)
                      cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTg                                                                                                  <  3:160527/148‑1 (MQ=255)
                        aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGaa                                                                                                <  1:570167/148‑1 (MQ=255)
                                  acaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcagc                                                                                      >  3:553794/1‑148 (MQ=255)
                                       ttACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGa                                                                                 <  4:553794/148‑1 (MQ=255)
                                               aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                       >  3:966165/1‑102 (MQ=255)
                                               aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                       <  4:966165/102‑1 (MQ=255)
                                                             cGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                       <  1:266785/88‑1 (MQ=255)
                                                             cGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                       >  2:266785/1‑88 (MQ=255)
                                                                  tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                    <  3:695998/102‑1 (MQ=255)
                                                                  tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                    >  4:695998/1‑102 (MQ=255)
                                                                  tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    >  4:782590/1‑134 (MQ=255)
                                                                  tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    <  3:782590/134‑1 (MQ=255)
                                                                   aaTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                       >  4:3009/1‑82 (MQ=255)
                                                                   aaTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                       <  3:3009/82‑1 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                          <  2:63027/74‑1 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                          <  3:398154/74‑1 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                          >  1:63027/1‑74 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                          >  4:398154/1‑74 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     <  2:953570/79‑1 (MQ=255)
                                                                        tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                     >  1:953570/1‑79 (MQ=255)
                                                                           aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                    <  2:166563/77‑1 (MQ=255)
                                                                           aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                    >  1:166563/1‑77 (MQ=255)
                                                                           aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                      >  2:592853/1‑91 (MQ=255)
                                                                           aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                      <  1:592853/91‑1 (MQ=255)
                                                                             atatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                      <  3:9130/89‑1 (MQ=255)
                                                                             atatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                      >  4:9130/1‑89 (MQ=255)
                                                                              tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                        >  2:323796/1‑118 (MQ=255)
                                                                              tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                        <  1:323796/118‑1 (MQ=255)
                                                                                tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    <  1:976638/120‑1 (MQ=255)
                                                                                tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    >  2:976638/1‑120 (MQ=255)
                                                                                     ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                    >  3:1111225/1‑67 (MQ=255)
                                                                                     ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                    <  4:1111225/67‑1 (MQ=255)
                                                                                     ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAgg                                                                                           >  4:224180/1‑92 (MQ=255)
                                                                                     ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAgg                                                                                           <  3:224180/92‑1 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGTGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                           <  3:772900/106‑1 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                           >  4:772900/1‑106 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAggggg                                                      >  3:343603/1‑127 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAggggg                                                      <  4:343603/127‑1 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcgg                                                   >  3:332300/1‑130 (MQ=255)
                                                                                       aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcgg                                                   <  4:332300/130‑1 (MQ=255)
                                                                                         aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    >  2:298512/1‑111 (MQ=255)
                                                                                         aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    <  1:298512/111‑1 (MQ=255)
                                                                                                  tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                      <  2:1075895/100‑1 (MQ=255)
                                                                                                  tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                      >  1:1075895/1‑100 (MQ=255)
                                                                                                   cAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                              <  1:396891/107‑1 (MQ=255)
                                                                                                   cAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                              >  2:396891/1‑107 (MQ=255)
                                                                                                    aaTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    >  3:888680/1‑100 (MQ=255)
                                                                                                    aaTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                    <  4:888680/100‑1 (MQ=255)
                                                                                                       gTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCt                 >  3:1112901/1‑148 (MQ=255)
                                                                                                           ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCCATg                                           <  2:1014166/118‑1 (MQ=255)
                                                                                                           ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                           >  1:1014166/1‑118 (MQ=255)
                                                                                                              aaaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgcg          <  4:1112901/148‑1 (MQ=255)
                                                                                                               aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                       <  2:1099468/86‑1 (MQ=255)
                                                                                                               aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                       >  1:1099468/1‑86 (MQ=255)
                                                                                                                      aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                       <  1:182097/95‑1 (MQ=255)
                                                                                                                      aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCt                                      <  3:85565/112‑1 (MQ=255)
                                                                                                                      aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCt                                      >  4:85565/1‑112 (MQ=255)
                                                                                                                      aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATcc  >  2:808641/1‑148 (MQ=255)
                                                                                                                      aTTCTTTGCTAATATGTTCGATAACGAACATGTATGAGCTGTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                       >  2:182097/1‑95 (MQ=255)
                                                                                                                               tAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATa                                                           <  4:651870/82‑1 (MQ=255)
                                                                                                                               tAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATa                                                           >  3:651870/1‑82 (MQ=255)
                                                                                                                                aaTATGTTCGATAACGAACCTTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                               <  3:908231/61‑1 (MQ=255)
                                                                                                                                aaTATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                               >  4:908231/1‑61 (MQ=255)
                                                                                                                                  aatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                       <  2:266904/66‑1 (MQ=39)
                                                                                                                                  aatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                       >  1:266904/2‑67 (MQ=39)
                                                                                                                                         gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAcc                                                                             >  1:647514/1‑54 (MQ=255)
                                                                                                                                         gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAcc                                                                             <  2:647514/54‑1 (MQ=255)
                                                                                                                                          aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                        <  1:239312/58‑1 (MQ=255)
                                                                                                                                          aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                        >  2:239312/1‑58 (MQ=255)
                                                                                                                                             |                                                                                                                            
TCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCC  >  NC_000913/3561831‑3562096

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: