Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,123,556 C→T E167K (GAA→AAA)  grxB ← reduced glutaredoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,123,5560CT100.0% 163.2 / NA 52E167K (GAA→AAA) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/32);  total (20/32)

ACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTC  >  NC_000913/1123310‑1123796
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
acataaCGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAGCCCGCGAAGACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCa                                                                                                                                                                                                                                               <  2:466158/247‑1 (MQ=255)
 cacaACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCAt                                                                                                                                                                                                                                              <  1:166063/251‑1 (MQ=255)
      cGCAGCATGTGGCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACg                                                                                                                                                                                                                                         <  2:42762/251‑1 (MQ=255)
                    ccAGCGATAAACTCGCCAGCAGAAAAAACGGAAAACGGCGGCAAACCGCGAATTCAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGGGTCTGTTTCGCCATATTATCGCGGCAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTGCGGTTTGAc                                                                                                                                                                                                                           <  2:437468/251‑1 (MQ=255)
                           tAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGAATAAAACGCCATTGATGATAAAAAATTTATTTTTGTCTGTTTCGCCCAATTATCGCGGGAATCAGCAAAGCGGCTTGGCCCGGTAATTCCGGCTACCCGCGGCAGATTTCGCCGTAGCGGGGAAAACTGGATATCCTCTTCCCAAAATTTTCCATTCACCGCGGTCgggttttcgatcatg                                                                                                                                                                                                                    >  2:364647/1‑238 (MQ=255)
                            aaaCTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGt                                                                                                                                                                                                                    <  1:697311/250‑1 (MQ=255)
                                   ccAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAg                                                                                                                                                                                                            <  1:364647/251‑1 (MQ=255)
                                        agAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCAc                                                                                                                                                                                                       <  1:589631/251‑1 (MQ=255)
                                                        ggctgcAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGGCCAGCGCACGTAAATCATCGCTGat                                                                                                                                                                                       <  1:455606/247‑1 (MQ=255)
                                                          cggcATCCCGCGAATACAGGCTGCCGCGTCGGGTCAGGATTAAATCGCCATTGAGGATAACAAAGTGATTTGTGTCTGTTGCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGatat                                                                                                                                                                                     <  2:103612/251‑1 (MQ=255)
                                                                 ccGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAAt                                                                                                                                                                              <  1:522629/251‑1 (MQ=255)
                                                                  cgcgAATCCAGGCTGCCGCGGCGGGTCAGGATTAAACTGCCATTGATGATAACAAATTTATTTGTGTCTTTTTCGCCTTATTATCGCGTTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATAATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAAta                                                                                                                                                                             <  2:163590/251‑2 (MQ=255)
                                                                  cgcgAATACAGGCTGCCGCGGCGGGTCAAGATTAAATCGCCATTGGTGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTTGCCAGTTAATTTCGGCTACCAGCGGCAGAATACGCAGTAGCTGGAATAGCCTAATATCATCTTCCCAAAGTTTTCCATTCACGGGGGTCGGGTTGGCGATCCGTTTTTTCCGCGCACGGAAAACATTGCTGatttttttcaa                                                                                                                                                                              >  1:703349/1‑241 (MQ=255)
                                                                    cgAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAAAAAATTGATTTGTGGCTGTTTTGCCATATTATCGAGGGAATCAGCAAAGCGGCTTGGGCAGTTAATTCCGGCTACCAACGTCAGAGTACGCAGTAGCGGGGAAAGATGAATATCATCTTCCGAAAGGTTGGCAATCACGGCGGTCGGGTTGGCGGACAGGTTGTCCAAGGCACGGAAATCATAGCTGATATTCTTAAta                                                                                                                                                                             >  1:278565/1‑248 (MQ=255)
                                                                                  cgcgGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGcc                                                                                                                                                                                                                                                <  1:376332/168‑1 (MQ=255)
                                                                                   gcggcgGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGc                                                                                                                                                           <  1:18329/251‑1 (MQ=255)
                                                                                    gggcgtgtcatgattaaaccgaacttGATGAAATCATATTGATTATTGGCTGTTTCGACTTATTAGCGCGGGAAACTGCTACGCGTCTTGGGAAGTTAATTCCTGCTACCAGCGTCATATTACGCCGTAGCGGGAATAGCGGAATATAATCTCACGAAAGTGTGCCACTCACTGCGTTCGGTTTGATGACCAGTTTGTCCAGCGTACGTAAATCATCGGGGATATTCTTAATCAGACCGTCAG‑AGTGGGcc                                                                                                                                                          <  2:697809/227‑1 (MQ=255)
                                                                                            cAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTc                                                                                                                                                  >  1:208721/1‑251 (MQ=255)
                                                                                                    aaTCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAAt                                                                                                                                          >  1:496030/1‑251 (MQ=255)
                                                                                                                  atAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGcc                                                                                                                            >  1:126474/1‑251 (MQ=255)
                                                                                                                   tAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCct                                                                                                                           >  1:349464/1‑251 (MQ=255)
                                                                                                                   aaaCAAATTGATTTGGGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCct                                                                                                                           <  1:362167/250‑1 (MQ=255)
                                                                                                                           tGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGCATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGACtctttcttt                                                                                                                   >  2:298620/1‑250 (MQ=255)
                                                                                                                           tGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCtctttctt                                                                                                                    <  1:40304/250‑1 (MQ=255)
                                                                                                                               ttGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTCCGCAGTAGCGGGAATAACTGAATATCATCTTCCGAAAGTTTGCCCTTCACGGCGTTCGGTTTGACGATCATTTTGTCCAGCGAACGTAAATCATCGCTGATATTCTTAATCGGAACGTCAGTAG‑GGGCCAGCAGGCCGGCAAAATTACCCGCGCTGGCCTCCTTCTTTTcgc                                                                                                               >  2:690232/1‑250 (MQ=255)
                                                                                                                               ttGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTcg                                                                                                                >  1:860480/1‑250 (MQ=255)
                                                                                                                                  tgtCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGGTGGCCTCTTTCTTTTcgacgg                                                                                                            >  2:865709/1‑250 (MQ=255)
                                                                                                                                     cTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGTAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAat                                                                                                         >  2:721270/1‑251 (MQ=255)
                                                                                                                                          tCGCTATATTATAGCGGTAATCAGCAACGCGGCTTGGCCAGGTAATTCCGGCTACCAGCGTCAGATGACGCAGTATCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTGGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATAttt                                                                                                     <  2:382683/250‑1 (MQ=255)
                                                                                                                                                    atCGCGGTAATCTGCAACGCGGTTGGGCCAGTTAATTATGGCTACAAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCAGCGTCCGAAAGTTTTCCCTTCACGGCTTCCGGTTTGACGTGCTGTTTGTCCAGCGCACGGAAATCCTCGCTGATATTCTTAATCATACCGTCAG‑AGTGGGGCAGTAGGTCGGCAAAATTACCCGCGCTGGCCTCTTGCTTGTCGACGAAATATTTGCGCgcggcg                                                                                           <  2:394115/250‑1 (MQ=255)
                                                                                                                                                         ggTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCCGGTCg                                                                                                                                                 >  2:790070/1‑191 (MQ=255)
                                                                                                                                                         ggTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCg                                                                                                                                                 <  1:790070/191‑1 (MQ=255)
                                                                                                                                                            aaTCAGCAACGAGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCCTCTTCCGCAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGGTTGTCCAGAGCACGTAAACCATCGCTGCTATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTAGGCAAAATTACCCGCGCTTGCCTCTTTCTTGTCGACGAACTATTTGCGCGCGGCGGGAGgaga                                                                                   >  2:603991/1‑250 (MQ=255)
                                                                                                                                                             aTCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCCCGGCGTTCGGTTTGACGATCAGTTTGTACAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCCCTTTCTTGTCGGCGAAATATTTGCGCGCGGCGGGAGGAGaa                                                                                  >  2:442866/1‑250 (MQ=255)
                                                                                                                                                               cAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAAc                                                                               <  1:630011/251‑1 (MQ=255)
                                                                                                                                                                  cAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGACATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGGGTAGAAAAcccc                                                                            >  2:183454/1‑248 (MQ=255)
                                                                                                                                                                   aaCGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACtcat                                                                           <  1:375171/251‑1 (MQ=255)
                                                                                                                                                                        ggCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAAt                                                                      >  1:405098/1‑251 (MQ=255)
                                                                                                                                                                               ccAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGAtt                                                               <  1:8976/251‑1 (MQ=255)
                                                                                                                                                                               ccAGTTAATTCCGGCTACCAGCGTAAGATTACGCAGTAGCGGGAATAGCTTAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGAtt                                                               <  1:729036/251‑1 (MQ=255)
                                                                                                                                                                                    tAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCCCGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGAAGAAGACAAATCATCAAAttcccgtttggcc                                                          >  2:267265/1‑239 (MQ=255)
                                                                                                                                                                                               aCCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACtc                                                                             <  1:359375/221‑1 (MQ=255)
                                                                                                                                                                                                 cAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAAca                                             <  1:865709/251‑1 (MQ=255)
                                                                                                                                                                                                  aGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAAcag                                            <  2:208721/251‑1 (MQ=255)
                                                                                                                                                                                                         gATTACGCAGTAGCGGGAAGTGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTttg                                     <  2:351924/251‑1 (MQ=255)
                                                                                                                                                                                                           ttACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGGAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTttgtt                                   >  1:546079/1‑251 (MQ=255)
                                                                                                                                                                                                                cAGTAGCGGGAATAGCTGAATATCATCTTCCAAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGt                              <  1:863989/251‑1 (MQ=255)
                                                                                                                                                                                                                       gggAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCAtt                       <  1:464309/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gAATATCATTTTCCGAAAGTTTGCCATTCACGGCTTTCGGTTTTCCGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTgc               <  2:682269/249‑1 (MQ=255)
                                                                                                                                                                                                                                   atatCATCTTCCGAAAGTTTGCCATGCACGGCGTGCGGGTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAAGTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAg           <  2:500900/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCAccct    >  2:683500/1‑248 (MQ=255)
                                                                                                                                                                                                                                             ccGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTc  <  1:621797/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAG‑AGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTC  >  NC_000913/1123310‑1123796

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