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→CGG) menI1,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

AGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGA  >  NC_000913/1765114‑1765555
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
aGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCCTCAATCGTTGCGAACAACACAAACGCCCTTTCCCATCGGAGACTTCAATCTTCCCAACCTGGGGACGCGAAACCAG‑AATGAACCGTTTTCAAACACCCCGCAc                                                                                                                                                                                                                          >  1:637178/1‑228 (MQ=255)
              ccTCGCTTTTAGACAATGCCTGCACAGTATACCACCGTTTACTGTGCAGAAAATGACCAAAAGCAATAAGCGTCACACTTTTCTGGTGACCACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTg                                                                                                                                                                                                           <  2:314798/228‑1 (MQ=255)
                              gACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGACAATCGTGACGAACAGGACAAACCGCCCTTCTCCTCGAAGATTTTCATCTGCCAGAACTGGTGACCCGAAACGCA‑ATTCAAAGGTTTGCCAACGCCGCGCCCCCCCCCGTTCTCGTGcaggcgccactggt                                                                                                                                                                                             >  1:63524/1‑214 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCCGACCCGGTGACGCGAACCGAG‑ATGCCACGAGTTGCCTACGCCCCGCCCCCCCCC‑TTCTCCTGCCCAGCGGACGTGGTTAGCA‑TTTAtt                                                                                                                                                                                >  2:454579/1‑227 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCCCGCACCCGCCC‑TTCTCATGCGGCGTGGACGGGGTTAGCA‑TTTAtt                                                                                                                                                                                >  2:375810/1‑227 (MQ=255)
                                         tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAAACGTGCAGCACAGCACAAAAGCCCTTTCTCATCGAAAGTTTCAATCTCCCCGACCTGGTGACGCGAAACGAG‑ACGCACCGGCTTTCCTACGCCCCGGCCCCCCCC‑CTCTCCTGCCCAGccggaggggttcggcttttgt                                                                                                                                                                                 >  1:584122/1‑207 (MQ=255)
                                            accaccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTTCCAGACCCGGGTACCCAAACCACGTATCCACGGGTTTGTAAAC‑CCCCCCCCCCCCCC‑CTCTCTTGCCCGGCGGCAGCGGTTAACC‑TTGATTTc                                                                                                                                                                              >  1:490591/1‑226 (MQ=255)
                                              caccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAAAGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCCTATGCCCCGCACCCGCCC‑TTCTCGTGCCGAGGCGACGGGGTTTGCA‑TTGCTTTCCAGa                                                                                                                                                                          >  2:555987/1‑228 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCa                                                                                                                                                                                                                           <  1:559021/176‑1 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCAAACGCCGCGCa                                                                                                                                                                                                                           >  2:559021/1‑176 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACCCGAACCCAA‑ATGCAACGGTTTGCCTACGCCCCGCCCCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCATaccca                                                                                                                                                                      >  2:148186/1‑225 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAAAGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGGGACGCGAACCCAG‑ATGCAACGGGTTGCATACGCCCCGCAACCCCCC‑TTCTCGTGCCGCGCGGGCCGGGTTAGCA‑TTGACTTCACGacccca                                                                                                                                                                     >  2:282657/1‑225 (MQ=255)
                                                   tttACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAAACCTGACGAACAGCACAAAAGCCCTTTCTCATCCACGAGTTTAAACTTCCAGACCTCGTGGCGCGACCCGAG‑GTCCACAGGGTTGCTTACGCCCCCCACCCCCCC‑TCCCCCGTCCCAGCGGACGCGCGTGGCA‑TTTAGTTCCaagccca                                                                                                                                                                      >  1:440464/1‑221 (MQ=255)
                                                          tgCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACACCACAAACGCCCCTTCTCCTCGAAAATTTCAATCTGCCAGACCTGGTGACGCGAACCCAG‑ATGCAACGGGTTTCCTACGCCGCCCCCCCCCCC‑TTCCCCTTCCGACCGCACGGGGGTAGCA‑TGAGTTTCCACAACAACCCCcttt                                                                                                                                                              >  1:285774/1‑228 (MQ=255)
                                                              aaTAATGACCAAAAGCAATTGGCGTTCCCCTTTTGTGGGGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCCCAAGCGCCCTTTCCCATCGAATATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCt                                                                                                                                                          <  1:555987/227‑1 (MQ=255)
                                                               aTAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTTAATCGTGACGAACAGCAACAACCCCCTTTCTCATCGAAGATTTCAATCTGCCAGGCCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTc                                                                                                                                                         <  2:581429/228‑1 (MQ=255)
                                                                    tcccAAAAGCTATTTGCGTCACGTTTTTCTGGTGACAACGTCACAAAATGGCGGTTGTCAATCGTGACGAACAGACAAAACGCCCTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTGCCGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  1:100941/226‑1 (MQ=255)
                                                                    gtccAAAAGCACTATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGATCCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  2:570045/226‑1 (MQ=255)
                                                                    gACCAAAGGAAATATGCGTCCCCCATTTCTGTGGACAACGCCCAAAAATGGCGTCCGTCAATCTGGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCtt                                                                                                                                                    <  1:148186/228‑1 (MQ=255)
                                                                                  tGCGTCACACTTTTCTGGTGACAACGGCACAAAATGGCGGTAGTCCAGCGTGACGAACAGCACAATCGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCCCTTTTTGCTCACCTTCGGTACataaata                                                                                                                                       <  2:637178/227‑1 (MQ=255)
                                                                                        acacTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGCGACGAACAGCCCAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAAAACCGGGGACCCCAACCGAG‑ATTCAACGGTTTGCATACGCCCCCCACCCCCCC‑CTCTCCTTGCCAACGGGAGGGGGTAGCATTTGTTTACCAAACAACCCCCTTTTTGCTCCCCTTCTG‑ACataaaaaaaagaca                                                                                                                                >  1:569066/1‑222 (MQ=255)
                                                                                              ttCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTTTCATCGAAGATTTCAAACTGCCAAACCTGGTGAAGCGAACCGAG‑GAGCAACGGTTTGCATACGCCGCGCACCCGCCC‑CTCTCTGGCCGAGCGGAAGCGGTTAGCG‑TTGATTTCCCAACACACCCCTTTTTTCTCACACTCGGTACatcaaacacgccaccgaac                                                                                                                           >  1:462921/1‑213 (MQ=255)
                                                                                                      aCAACGTCACAAAATGGCGGTCGGCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGGTTGCATACGCCCCGCACCCGCCC‑CTCTCGTGCCGGGCGGACGTGGTTAGCA‑TTGTTTTCCAGACCAACCACCTTTTGCTCCCCTTCGGTCCATAAAAAACATGCAACGCAAccagaaca                                                                                                                  >  1:632530/1‑222 (MQ=255)
                                                                                                            gcacaAAATGGCGGTCGTCAATCGTGACCAACAGCACAAACGCCCTTCCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGc                                                                                                            <  2:4075/227‑1 (MQ=255)
                                                                                                                       ccggcgtcAATCTTGACGAACAGCCCAAACGCCCTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTTACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAA‑‑AACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               <  2:462921/224‑1 (MQ=255)
                                                                                                                         gtcgtcAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGCCCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               <  2:24934/228‑1 (MQ=255)
                                                                                                                               aattggtgcagaccggctaaccGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGATGCGAACCGAG‑ATGCAACGGTTTGCTTACGCCGCGCCCCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGctcct                                                                                         <  1:180836/207‑1 (MQ=255)
                                                                                                                                    tGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCACCCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGCACCGATACTTTCGGCCATTACCACAGATGCCCCCCCaac                                                                                    >  1:524128/1‑226 (MQ=255)
                                                                                                                                       cgaaaacaacaAACGCCCTTTCTCATCGAAGAGTACAAGCGGCCAGACCGGGGGACGCGAACCGAG‑ATGCAACTGTTGGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcag                                                                                 <  1:81869/220‑1 (MQ=255)
                                                                                                                                         aaCAGCACAAACGCCTTTTCCCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcagca                                                                               <  2:285774/228‑1 (MQ=255)
                                                                                                                                                   aCGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCgg                                                                                                             >  2:292622/1‑188 (MQ=255)
                                                                                                                                                    cGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGc                                                                                                            <  1:292622/188‑1 (MQ=255)
                                                                                                                                                    cGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAAAAACCGGCACCGGAACCGATCCTTTCGGCCGGAGCCACGAATGCCCTTCCAGCCGGACACCCGAAAgg                                                                     >  2:620961/1‑227 (MQ=255)
                                                                                                                                                        cTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTgg                                                                                                                                                                                             >  2:461098/1‑104 (MQ=255)
                                                                                                                                                         tttCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGt                                                                                                                                                                                            <  1:461098/104‑1 (MQ=255)
                                                                                                                                                             taaacGAAGATTTCAATCTGCCAGGCCGGGGGCCGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCCCCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTc                                                           <  2:490591/224‑1 (MQ=255)
                                                                                                                                                                agggaGATTTCAAATGGCCAGACCTGGTGACGCGAACCGGG‑AGGCAACGGGTTGCAAACGCCCGCCACCCGCCC‑TTCCCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTGCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCcgc                                                        <  1:454579/224‑1 (MQ=255)
                                                                                                                                                                 cGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCAATGCAGAAACCCGAAAGGCTGCTTTGTCcgcg                                                       >  2:507158/1‑228 (MQ=255)
                                                                                                                                                                       tttCAATCCGCCAGGCCTGGGTACGCGAACCGAA‑ATGCAACGGTTTGACTACCCCGCCCCCCCGCCC‑TTCTCGTTCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTa                                                 <  2:642359/228‑1 (MQ=255)
                                                                                                                                                                                gCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                                                               >  2:306179/1‑173 (MQ=255)
                                                                                                                                                                                 ccAGACCTGGGGACGCGAACCGAG‑ATGCAACGGTTTGCATACCCCGCGCCCCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATg                                                                                              <  1:306179/173‑1 (MQ=255)
                                                                                                                                                                                         ggTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGCCGTGGTCAGCA‑TTGATACCCAGCCCAACCCATTTTTGCTCACCTTCGGTACCTAAATAACCGGCAACGGAACCGATACTTTCGGCACCTACCACGGATGCCCCGCCATTCAGCACCCCTAACTGCAGCTTTTTCCGCGAGTGTCCAGCCATTGTCGCTTCa                                >  1:514508/1‑227 (MQ=255)
                                                                                                                                                                                                          aTGCAACGGTTTGCAAACGCCCACCACCCGCCC‑TACTGGTGCCCGGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCTCTTTTTGCTCCCCTTCGGTACATAAATAACCGGCAACGGAACCGATCCTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAat               <  2:17345/228‑1 (MQ=255)
                                                                                                                                                                                                                     tGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTTTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCCCCACGATGTTCAAAGc    >  2:147720/1‑228 (MQ=255)
                                                                                                                                                                                                                       actACGCCGCGCCCCCGCCC‑TTCTCGTGCCGCGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCTCTTGTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATCCTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGa  <  2:569775/226‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
AGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCGAG‑ATGCAACGGTTTGCATACGCCGCGCACCCGCCC‑TTCTCGTGCCGAGCGGACGTGGTTAGCA‑TTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACTGGCATTGTCGCTTCAAGGGTGTCATCACCAATATGTTCAAAGCGA  >  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: ‑