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,8230GA95.6% 129.2 / ‑5.6 45L78L (CTG→TTG) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (27/16);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
*NC_0009131,123,8240TG97.7% 139.9 / ‑5.2 45L77F (TTA→TTC) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (27/16);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123582‑1124056
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                          
aTCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAAAGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGttt                                                                                                                                                                                                                                    >  1:67848/1‑251 (MQ=255)
              cgcACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTAt                                                                                                                                                                                                                      <  1:379526/251‑1 (MQ=255)
                           tcGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCCGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACCGCAGTTTGCTGGCGTAGCCATTGCCCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTCCAG‑AACGGTTTGCCGTCGAGTTTATCGCCATactgcac                                                                                                                                                                                                         >  1:93476/1‑245 (MQ=255)
                            cGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAAc                                                                                                                                                                                                         <  1:214769/250‑1 (MQ=255)
                             gCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCCGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGCGTTTATCGACATAGTGAACg                                                                                                                                                                                                        >  1:475240/1‑250 (MQ=255)
                              cTGATATTCTTAATCAGACCGTCAGAGGGGGCCAGCAGGTCGGCAAAATTACCCACGCTGGCCTCTTTCTTTTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCCAATGCCGCTTTGGCAAACCTCGGCAACCGCAGTTTGTTTACGTATCCATTTACCTTGCGCCGCCACTCTTCAATTGCCTGGGTCCGTTTTCCCGTCCAGTA‑CGGTTTGCCATCGAGTTTATCGACTTAGTgtactat                                                                                                                                                                                                      >  1:75671/1‑245 (MQ=255)
                                   atTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTcc                                                                                                                                                                                                  <  1:55953/250‑1 (MQ=255)
                                        tAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCtt                                                                                                                                                                                            >  1:535659/1‑251 (MQ=255)
                                        tAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGGTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCtt                                                                                                                                                                                            >  1:51739/1‑251 (MQ=255)
                                                       agTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCgg                                                                                                                                                                             >  1:386016/1‑251 (MQ=255)
                                                                 cagGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCCACGAAATATTAGCGCGCGGCGGGCGTAGAAAACTCATCACATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGCGTTTATCGACATAGTGAACGATGTCCAAGCTTTCTGGCATATAGCGGCTG‑Tcttctt                                                                                                                                                                  >  1:259208/1‑251 (MQ=255)
                                                                                        tGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCCTTTAGCGGCTGGTC‑TTTTTTTTCAGTATTGGAACCTGtttt                                                                                                                                           >  2:108615/1‑251 (MQ=255)
                                                                                        gggtctttttcttctcagcgataattttGCGCGGGGCGGGAGTAAAAAACTCACCAAATGCCGATTTGTCAAACCGCGGCAAAAGCATTTTGTTGGCGTAGCCATTGACCTACCCCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGtttt                                                                                                                                           <  2:326435/226‑1 (MQ=255)
                                                                                           ccTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCGTCTTTTTGCAAAATGGGAACCTGTTTTTgt                                                                                                                                        >  1:122653/1‑250 (MQ=255)
                                                                                           ccTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTga                                                                                                                                        <  1:108615/251‑1 (MQ=255)
                                                                                            ctctTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCCATGGCAGGGGCACCTTTGCCGGCC‑AAGACCGGTTTGCCGTCGAGTTTATTGACATAGGGAACGAAGACCCTGCTCTCTGGCATATAGAGGCTG‑TCCCTTTTTTGTAGCATGGGAAACTGgtttgta                                                                                                                                        >  2:504621/1‑247 (MQ=255)
                                                                                                           gaAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCCTCCgggggg                                                                                                                        >  1:376284/1‑248 (MQ=255)
                                                                                                            aaaTATTGGCGCGCGGCGGGAGTAGAAAACTCATCAAAGGCGTATTGGGCAAAACGCGGCAACAGCAGTTTTTTGGCGTAGCCATTGACCTGGCGCAGCCACTCTTCAATGGCAGGGGAACGGTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCgggtggg                                                                                                                       <  2:475240/251‑1 (MQ=255)
                                                                                                                     cgcgcgGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGGGTTTCTGc                                                                                                              >  1:282203/1‑251 (MQ=255)
                                                                                                                                    gAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGc                                                                                                 >  1:6836/1‑249 (MQ=255)
                                                                                                                                      aaaCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTAAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TTATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGGGTTTCTGCGTCGCCGTTGAGCAgt                                                                                              >  2:552842/1‑249 (MQ=255)
                                                                                                                                        aCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGAAATAGTGAACGATGTCCCTGCTTTCTGGCAGATTAAGGATG‑TCATTTTTTTGCCGGAAGGGGACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTTCGTCGGTGTTGTGCAGaca                                                                                            >  2:423237/1‑248 (MQ=255)
                                                                                                                                          tcatcaAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGAGGCTG‑TCATCTTTTTGCAGAATGGGAACATGTTTTTGGCCGACCATCCGGGTGGGGGTTTCTGAGTCGGCGGGGAGCAGA‑ACa                                                                                          >  2:367924/1‑249 (MQ=255)
                                                                                                                                              aaaaagccattttGGCAAAACGCGGCAACAGCAGTTTGTTGGGGTAGCCATTGATCTGGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTAGCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAAGGGGAACCTGTTGTTGACCTACCATCCTGGTGGGTGTTTCTGCGTCGGCGTTGAGCAGA‑TCATGTAAt                                                                                    <  2:386016/241‑1 (MQ=255)
                                                                                                                                                 aTGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCg                                                                                 <  1:585762/251‑1 (MQ=255)
                                                                                                                                                                          cagcagTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCCTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATCGTGAACGCTGTCCATGCTTTCTGGCATGTAGAGGCTG‑TCATCTTTTTGCAGAATGGGCACCTGTTTTTTACCGACCCACCGGGAGGGGTTTTCTGCGTCGTAGTTGAGCAAC‑ACATGTAACTCGTCGGGGATATTTTTCCTGCCGaaa                                                         >  2:228850/1‑250 (MQ=255)
                                                                                                                                                                           agcagTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAAACTGTTTTTGACCGACACTCCGGGGGGGGGTTTCTGCGTCGTCGTTGAGCTGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAAt                                                       >  1:66174/1‑251 (MQ=255)
                                                                                                                                                                              ggTTTGTTGGCGTGGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCAt                                                    <  2:6836/250‑1 (MQ=255)
                                                                                                                                                                                     tgGCGTAGACATTGCCCTTGAGCAGCCACTCTTCAATGACAGGGGAACGGTTGCAGGTC‑AAGACCGGTTTGCCGTCGAGTTTATCGACCTGATGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCCTTTTGCAGAATGCGAACCCGTTTTTGGACGACCATCCGGGAGGGGGTTTCAGCGTGGGAGTGGTGCGGA‑ACATGTAATTCGGTGGGGGTATTTTTCAGGCCGaaacaccttcagga                                              >  2:655174/1‑239 (MQ=255)
                                                                                                                                                                                       gCGTAGCCTTTGACTTTGCCCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATAATTTTCAGGCCGAAAATCATGCGGTCttt                                           <  2:40040/251‑1 (MQ=255)
                                                                                                                                                                                          tAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTg                                          >  1:563896/1‑249 (MQ=255)
                                                                                                                                                                                            cacattaccctttgtcggcctcttttCAAGGGCAGGAGAACGTTGGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTAGCGACATAGTTAACGATTTCCATGAGTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGAACGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGc                                      <  2:304605/227‑1 (MQ=255)
                                                                                                                                                                                              cATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGc                                      >  1:294528/1‑249 (MQ=255)
                                                                                                                                                                                               aTTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTAGCGACATAGTGAACGAGGTCCATGCTTTCTGGCATATAGCGGCTG‑TCAACTTGTTGCAGAATGGGAACCTGTTTTTGGCCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGt                                   >  1:545407/1‑251 (MQ=255)
                                                                                                                                                                                                  gtcctgccccAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTGTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTaa                                 <  2:563896/241‑1 (MQ=255)
                                                                                                                                                                                                  gACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGGGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGGCGGGGATATTTTTCAGGCCGGAAATAATGCGGGCTTTGAGTCAGTa                                  >  2:529703/1‑249 (MQ=255)
                                                                                                                                                                                                   aCCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAgg                               >  1:514420/1‑251 (MQ=255)
                                                                                                                                                                                                       tGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTGGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGTGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCa                            <  1:655174/250‑1 (MQ=255)
                                                                                                                                                                                                         cgcAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTg                         >  1:93355/1‑251 (MQ=255)
                                                                                                                                                                                                          gcAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAg                                        <  1:460257/235‑1 (MQ=255)
                                                                                                                                                                                                           cAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGAt                       <  1:610729/251‑1 (MQ=255)
                                                                                                                                                                                                                     tCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGCACCTGTTTTTGACCGACCATCCGGGTGGGGGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGt              >  1:393841/1‑250 (MQ=255)
                                                                                                                                                                                                                            cAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAgg                                       <  1:503145/218‑1 (MQ=255)
                                                                                                                                                                                                                              ggggAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGCATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCAc    >  1:625171/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggtccGTTTGCCAGTC‑AAGAACGGATTGACGTCGAGTTTATCGACATAGCGACCGAAGACCATGCTTTCTGGCCTATAGCGGCTG‑TCATCTTTTTGCAGAATTGGCACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGATGCGGAGAGAC‑TGTAATTCGACGGGGATATTTTTCAGTCCGAAAATCATGCGGGCTTTGAGGCATTAAGGGCATTGCTCGAAAATGTATCGATTCACGt  >  2:324028/5‑251 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                          
ATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTG‑TCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGA‑ACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123582‑1124056

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