Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,765,336 G→C 26.1% R99R (CGC→CGG menI ← 1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,765,3360GC26.1% 79.7 / 2.1 45R99R (CGC→CGGmenI1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base G (12/20);  new base C (11/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 7.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGA  >  NC_000913/1765114‑1765555
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
aGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCCTCAATCGTTGCGAACAACACAAACGCCCTTTCCCATCGGAGACTTCAATCTTCCCAACCTGGGGACGCGAAACCAGAATGAACCGTTTTCAAACACCCCGCAc                                                                                                                                                                                                                          >  1:637178/1‑228 (MQ=255)
              ccTCGCTTTTAGACAATGCCTGCACAGTATACCACCGTTTACTGTGCAGAAAATGACCAAAAGCAATAAGCGTCACACTTTTCTGGTGACCACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTg                                                                                                                                                                                                           <  2:314798/228‑1 (MQ=255)
                              gACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGACAATCGTGACGAACAGGACAAACCGCCCTTCTCCTCGAAGATTTTCATCTGCCAGAACTGGTGACCCGAAACGCAATTCAAAGGTTTGCCAACGCCGCGCCCCCCCCCGTTCTCGTGcaggcgccactggt                                                                                                                                                                                             >  1:63524/1‑214 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCCGACCCGGTGACGCGAACCGAGATGCCACGAGTTGCCTACGCCCCGCCCCCCCCCTTCTCCTGCCCAGCGGACGTGGTTAGCATTTAtt                                                                                                                                                                                >  2:454579/1‑227 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCCCGCACCCGCCCTTCTCATGCGGCGTGGACGGGGTTAGCATTTAtt                                                                                                                                                                                >  2:375810/1‑227 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAAACGTGCAGCACAGCACAAAAGCCCTTTCTCATCGAAAGTTTCAATCTCCCCGACCTGGTGACGCGAAACGAGACGCACCGGCTTTCCTACGCCCCGGCCCCCCCCCTCTCCTGCCCAGccggaggggttcggcttttgt                                                                                                                                                                                 >  1:584122/1‑207 (MQ=255)
                                            accaccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTTCCAGACCCGGGTACCCAAACCACGTATCCACGGGTTTGTAAACCCCCCCCCCCCCCCCTCTCTTGCCCGGCGGCAGCGGTTAACCTTGATTTc                                                                                                                                                                              >  1:490591/1‑226 (MQ=255)
                                              caccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAAAGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCCTATGCCCCGCACCCGCCCTTCTCGTGCCGAGGCGACGGGGTTTGCATTGCTTTCCAGa                                                                                                                                                                          >  2:555987/1‑228 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCa                                                                                                                                                                                                                           <  1:559021/176‑1 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCAAACGCCGCGCa                                                                                                                                                                                                                           >  2:559021/1‑176 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACCCGAACCCAAATGCAACGGTTTGCCTACGCCCCGCCCCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCATaccca                                                                                                                                                                      >  2:148186/1‑225 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAAAGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGGGACGCGAACCCAGATGCAACGGGTTGCATACGCCCCGCAACCCCCCTTCTCGTGCCGCGCGGGCCGGGTTAGCATTGACTTCACGacccca                                                                                                                                                                     >  2:282657/1‑225 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAAACCTGACGAACAGCACAAAAGCCCTTTCTCATCCACGAGTTTAAACTTCCAGACCTCGTGGCGCGACCCGAGGTCCACAGGGTTGCTTACGCCCCCCACCCCCCCTCCCCCGTCCCAGCGGACGCGCGTGGCATTTAGTTCCaagccca                                                                                                                                                                      >  1:440464/1‑221 (MQ=255)
                                                          tgCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACACCACAAACGCCCCTTCTCCTCGAAAATTTCAATCTGCCAGACCTGGTGACGCGAACCCAGATGCAACGGGTTTCCTACGCCGCCCCCCCCCCCTTCCCCTTCCGACCGCACGGGGGTAGCATGAGTTTCCACAACAACCCCcttt                                                                                                                                                              >  1:285774/1‑228 (MQ=255)
                                                              aaTAATGACCAAAAGCAATTGGCGTTCCCCTTTTGTGGGGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCCCAAGCGCCCTTTCCCATCGAATATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCt                                                                                                                                                          <  1:555987/227‑1 (MQ=255)
                                                               aTAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTTAATCGTGACGAACAGCAACAACCCCCTTTCTCATCGAAGATTTCAATCTGCCAGGCCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTc                                                                                                                                                         <  2:581429/228‑1 (MQ=255)
                                                                    tcccAAAAGCTATTTGCGTCACGTTTTTCTGGTGACAACGTCACAAAATGGCGGTTGTCAATCGTGACGAACAGACAAAACGCCCTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTGCCGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  1:100941/226‑1 (MQ=255)
                                                                    gtccAAAAGCACTATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGATCCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  2:570045/226‑1 (MQ=255)
                                                                    gACCAAAGGAAATATGCGTCCCCCATTTCTGTGGACAACGCCCAAAAATGGCGTCCGTCAATCTGGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  1:148186/228‑1 (MQ=255)
                                                                                  tGCGTCACACTTTTCTGGTGACAACGGCACAAAATGGCGGTAGTCCAGCGTGACGAACAGCACAATCGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCCCTTTTTGCTCACCTTCGGTACataaata                                                                                                                                       <  2:637178/227‑1 (MQ=255)
                                                                                        acacTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGCGACGAACAGCCCAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAAAACCGGGGACCCCAACCGAGATTCAACGGTTTGCATACGCCCCCCACCCCCCCCTCTCCTTGCCAACGGGAGGGGGTAGCATTTGTTTACCAAACAACCCCCTTTTTGCTCCCCTTCTGACataaaaaaaagaca                                                                                                                                >  1:569066/1‑222 (MQ=255)
                                                                                              ttCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTTTCATCGAAGATTTCAAACTGCCAAACCTGGTGAAGCGAACCGAGGAGCAACGGTTTGCATACGCCGCGCACCCGCCCCTCTCTGGCCGAGCGGAAGCGGTTAGCGTTGATTTCCCAACACACCCCTTTTTTCTCACACTCGGTACatcaaacacgccaccgaac                                                                                                                           >  1:462921/1‑213 (MQ=255)
                                                                                                      aCAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGGTTGCATACGCCCCGCACCCGCCCCTCTCGTGCCGGGCGGACGTGGTTAGCATTGTTTTCCAGACCAACCACCTTTTGCTCCCCTTCGGTCCATAAAAAACATGCAACGCAAccagaaca                                                                                                                  >  1:632530/1‑222 (MQ=255)
                                                                                                            gcacaAAATGGCGGTCGTCAATCGTGACCAACAGCACAAACGCCCTTCCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGc                                                                                                            <  2:4075/227‑1 (MQ=255)
                                                                                                                       ccggcgtcAATCTTGACGAACAGCCCAAACGCCCTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTTACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAA‑‑AACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               <  2:462921/224‑1 (MQ=255)
                                                                                                                         gtcgtcAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGCCCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               <  2:24934/228‑1 (MQ=255)
                                                                                                                               aattggtgcagaccggctaaccGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGATGCGAACCGAGATGCAACGGTTTGCTTACGCCGCGCCCCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGctcct                                                                                         <  1:180836/207‑1 (MQ=255)
                                                                                                                                    tGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCACCCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGCACCGATACTTTCGGCCATTACCACAGATGCCCCCCCaac                                                                                    >  1:524128/1‑226 (MQ=255)
                                                                                                                                       cgaaaacaacaAACGCCCTTTCTCATCGAAGAGTACAAGCGGCCAGACCGGGGGACGCGAACCGAGATGCAACTGTTGGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcag                                                                                 <  1:81869/220‑1 (MQ=255)
                                                                                                                                         aaCAGCACAAACGCCTTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcagca                                                                               <  2:285774/228‑1 (MQ=255)
                                                                                                                                                   aCGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCgg                                                                                                             >  2:292622/1‑188 (MQ=255)
                                                                                                                                                    cGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGc                                                                                                            <  1:292622/188‑1 (MQ=255)
                                                                                                                                                    cGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAAAAACCGGCACCGGAACCGATCCTTTCGGCCGGAGCCACGAATGCCCTTCCAGCCGGACACCCGAAAgg                                                                     >  2:620961/1‑227 (MQ=255)
                                                                                                                                                        cTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTgg                                                                                                                                                                                             >  2:461098/1‑104 (MQ=255)
                                                                                                                                                         tttCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGt                                                                                                                                                                                            <  1:461098/104‑1 (MQ=255)
                                                                                                                                                             taaacGAAGATTTCAATCTGCCAGGCCGGGGGCCGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCCCCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTc                                                           <  2:490591/224‑1 (MQ=255)
                                                                                                                                                                agggaGATTTCAAATGGCCAGACCTGGTGACGCGAACCGGGAGGCAACGGGTTGCAAACGCCCGCCACCCGCCCTTCCCGTGCCGAGCGGACGTGGTTAGCATTGATTGCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCcgc                                                        <  1:454579/224‑1 (MQ=255)
                                                                                                                                                                 cGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCAATGCAGAAACCCGAAAGGCTGCTTTGTCcgcg                                                       >  2:507158/1‑228 (MQ=255)
                                                                                                                                                                       tttCAATCCGCCAGGCCTGGGTACGCGAACCGAAATGCAACGGTTTGACTACCCCGCCCCCCCGCCCTTCTCGTTCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTa                                                 <  2:642359/228‑1 (MQ=255)
                                                                                                                                                                                gCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               >  2:306179/1‑173 (MQ=255)
                                                                                                                                                                                 ccAGACCTGGGGACGCGAACCGAGATGCAACGGTTTGCATACCCCGCGCCCCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATg                                                                                              <  1:306179/173‑1 (MQ=255)
                                                                                                                                                                                         ggTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGCCGTGGTCAGCATTGATACCCAGCCCAACCCATTTTTGCTCACCTTCGGTACCTAAATAACCGGCAACGGAACCGATACTTTCGGCACCTACCACGGATGCCCCGCCATTCAGCACCCCTAACTGCAGCTTTTTCCGCGAGTGTCCAGCCATTGTCGCTTCa                                >  1:514508/1‑227 (MQ=255)
                                                                                                                                                                                                          aTGCAACGGTTTGCAAACGCCCACCACCCGCCCTACTGGTGCCCGGCGGACGTGGTTAGCATTGATTTCCAGACCAACCTCTTTTTGCTCCCCTTCGGTACATAAATAACCGGCAACGGAACCGATCCTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAat               <  2:17345/228‑1 (MQ=255)
                                                                                                                                                                                                                     tGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTTTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCCCCACGATGTTCAAAGc    >  2:147720/1‑228 (MQ=255)
                                                                                                                                                                                                                       actACGCCGCGCCCCCGCCCTTCTCGTGCCGCGCGGACGTGGTTAGCATTGATTTCCAGACCAACCTCTTGTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATCCTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGa  <  2:569775/226‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
AGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAGATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGA  >  NC_000913/1765114‑1765555

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