Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,123,823 2 bp→AG coding (231‑232/648 nt) grxB ← reduced glutaredoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,123,8230GA97.4% 117.5 / ‑5.3 39L78L (CTG→TTG) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/24);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
*NC_0009131,123,8240TG97.4% 116.6 / ‑4.2 39coding (231/648 nt)grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (14/24);  minor base . (1/0);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGTGACTCCTCC  >  NC_000913/1123580‑1124065
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                 
cGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGt                                                                                                                                                                                                                                             >  2:135305/1‑251 (MQ=255)
           cccAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAg                                                                                                                                                                                                                                  <  1:254693/250‑1 (MQ=255)
               gcgcACGTCAATCATCGCCGATATACTTAATCAGCGCGTCAGAGAGGGCCAACAGCTCGGCAAAATTAACAGCGCTGGCCTCTTTCTTGTCGACGCAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACCGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTCCCGGTCTAG‑AACGGTTTGCCGTAGAGttt                                                                                                                                                                                                                               >  2:77719/1‑250 (MQ=255)
               gcgcACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTa                                                                                                                                                                                                                              <  1:543891/251‑1 (MQ=255)
                               gCTGATCTTCTTAATCCGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCATCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGACGATTTGGCAAAAAGCGGCAACAGCAGTTTGAT‑GGCGGAGCCATTGACGTTGCGCAGCCACTCTTCAATGGCAGGGGAAAGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCggactaggggacgg                                                                                                                                                                                                              >  2:76123/1‑238 (MQ=255)
                                           aaTCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCttt                                                                                                                                                                                                  <  1:87001/251‑1 (MQ=255)
                                                 aCCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATAAAATGACGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGCGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGc                                                                                                                                                                                             >  1:532872/1‑250 (MQ=255)
                                                    gTCAGAGTGGGCCAGCAGGTCGGCAAAATTCCCCGCTCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGGGTATAAACCTCATCAAATGCCGATTTGGCAAAACTCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGCCCTTGCGCAGCCACTCTTCAATGGCTGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGTCatat                                                                                                                                                                                         <  2:211134/251‑1 (MQ=255)
                                                      cAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGGGTTTATCGACATAGTGAACGATGTCCCTGCTTTCTGGCatat                                                                                                                                                                                         >  2:474491/1‑249 (MQ=255)
                                                        gagTGGGCCAGCAGGTCGGCAAAATTTCCCGAGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCACCAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCg                                                                                                                                                                                     <  2:430091/251‑1 (MQ=255)
                                                            gggCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCCGTTTGTT‑GGCGAAGCCATTGACCTTGAGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTg                                                                                                                                                                                 >  1:241647/1‑251 (MQ=255)
                                                                     ggTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAAGGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCTTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCtttt                                                                                                                                                                        <  2:28879/251‑1 (MQ=255)
                                                                                   cccGCGCTTGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGGCATCTTTTGGCAGAATGGGAAc                                                                                                                                                          <  1:511697/251‑1 (MQ=255)
                                                                                      gcgcTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTg                                                                                                                                                       <  1:129992/251‑1 (MQ=255)
                                                                                                         cgacgaAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCg                                                                                                                                     <  1:442682/250‑1 (MQ=255)
                                                                                                                      gcgcgcGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTTACCGACCATCCGGGTGGGTGTTTCt                                                                                                                        <  1:173068/250‑1 (MQ=255)
                                                                                                                      gcgcgcGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTg                                                                                                                       <  1:160315/251‑1 (MQ=255)
                                                                                                                       cgcgcgGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGc                                                                                                                      <  2:171227/251‑1 (MQ=255)
                                                                                                                          gcggcgGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTGCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGcgtc                                                                                                                   <  1:95637/251‑1 (MQ=255)
                                                                                                                            ggcggGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGcgtcgt                                                                                                                 >  1:136572/1‑251 (MQ=255)
                                                                                                                                       aaaaCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGGGTTTCTGCGTCGCCGTTGAGc                                                                                                         >  1:12529/1‑248 (MQ=255)
                                                                                                                                         acctaaacaatggCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCAGTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTGTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAAc                                                                                                    <  2:355792/239‑1 (MQ=255)
                                                                                                                                            ccagcaAAGGCCGATTTGGCAAAACGCTGCAACAGCAGTTTGTT‑GGCGTAGCCATTGCCTTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGGATGGGAACCTGTTGTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATg                                                                                                 <  2:91663/247‑1 (MQ=255)
                                                                                                                                                 aaaTGCCGTTTTGGCAAAACGCGGCAACAGCA‑TTTGTTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTTCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAAtt                                                                                            <  1:330864/251‑1 (MQ=255)
                                                                                                                                                        aaTTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCATGACCTGTAATTCGACggg                                                                                     >  2:185143/2‑251 (MQ=255)
                                                                                                                                                                   ccgcgGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCATTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAAGTCGACGGGGATATTTTTCa                                                                          <  2:474857/250‑1 (MQ=255)
                                                                                                                                                                      cgGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGGTATTTTTCAGGc                                                                       >  2:286646/1‑251 (MQ=255)
                                                                                                                                                                       ggCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTCTGCCGTCGAGTTTATCGACATAGTGACCGCTGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGACTTCGAACCTGTTTGTGACCGACCATCCGGGTGGGGGTTTCTGCGTCGTCGTTGCGCATTACATGTAATTCGCCAGGGTTCTTTTTCAGGcc                                                                      >  2:131577/1‑251 (MQ=255)
                                                                                                                                                                        gCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCg                                                                     >  1:183259/1‑251 (MQ=255)
                                                                                                                                                                                    tgtt‑gGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCg                                                          <  1:135305/250‑1 (MQ=255)
                                                                                                                                                                                          gCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCttt                                                    <  1:474491/251‑1 (MQ=255)
                                                                                                                                                                                                      aCCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGCGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAgg                                        >  1:313217/1‑251 (MQ=255)
                                                                                                                                                                                                                            tGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGGGATCGTAAATGTATAGc                  >  2:10573/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggCAGGGGAACTGTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACAGAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTGGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGGGATCGTAAATGTATAGCt                 <  2:409589/251‑1 (MQ=255)
                                                                                                                                                                                                                                aGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACAt                                                                                                  >  1:75565/1‑167 (MQ=255)
                                                                                                                                                                                                                                aGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGGCGTCGTTGAGCAGAACAt                                                                                                  <  2:75565/167‑1 (MQ=255)
                                                                                                                                                                                                                                 ggggAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCAc             <  1:76123/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGTTTGCCGGTC‑AAGAACGGTTTGCCGGCGAGTGTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGGTTCTGCGTCGGCGTTGAGCAGAACATGTAAGTCGACGGGGATATTGTTCAGGCCGAAAATCATGCGGGCTTTGAGTCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGTgtct       <  2:559081/251‑4 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTTAACGATGTCCTTGCGTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATTCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGTGACtcctcc  <  1:146932/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                 
CGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTT‑GGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGTGACTCCTCC  >  NC_000913/1123580‑1124065

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: ‑