Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,375,484 A→T 6.4% L160H (CTT→CAT) ‡ menC ← O‑succinylbenzoyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,375,4840AT 6.4% 216.1 / 5.3 62L160H (CTT→CAT) ‡menCO‑succinylbenzoyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base A (23/35);  new base T (2/2);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

GGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGC  >  NC_000913/2375336‑2375628
                                                                                                                                                    |                                                                                                                                                
ggcaatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                                                                  <  2:50117/149‑1 (MQ=255)
        cGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCaaa                                                                                                                                          >  2:111645/1‑149 (MQ=255)
          gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtc                                                                                                                                        <  1:262650/149‑2 (MQ=255)
          gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATg                                                                                                                                        >  2:346519/1‑149 (MQ=255)
               aCGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGAt                                                                                                                                   <  1:255155/149‑1 (MQ=255)
                     aaaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                                                             >  1:442132/1‑149 (MQ=255)
                      aaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAt                                                                                                                            <  1:353875/149‑1 (MQ=255)
                       aCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                                                                           <  1:111645/149‑1 (MQ=255)
                       aCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                                                                           >  3:157060/1‑149 (MQ=255)
                          cTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGcc                                                                                                                        <  2:240824/149‑1 (MQ=255)
                           tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCaaa                                                                                                                                          <  3:248844/130‑1 (MQ=255)
                           tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCaaa                                                                                                                                          >  4:248844/1‑130 (MQ=255)
                           tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                       >  3:468701/1‑149 (MQ=255)
                           tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                       >  3:351436/1‑149 (MQ=255)
                            cgcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCACATCCGGAATTGCCTc                                                                                                                      <  4:351436/149‑1 (MQ=255)
                             gcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTTACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTcc                                                                                                                     >  3:471719/1‑149 (MQ=255)
                             gcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTcc                                                                                                                     >  1:264336/1‑149 (MQ=255)
                              cgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCa                                                                                                                    <  2:348025/149‑1 (MQ=255)
                               gAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                                                                   <  1:346519/149‑1 (MQ=255)
                                      cgcgGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGa                                                                                                            <  3:270356/149‑1 (MQ=255)
                                      cgcgGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGa                                                                                                            >  1:104970/1‑149 (MQ=255)
                                        cgGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGGAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGAtt                                                                                                          <  3:138430/149‑1 (MQ=255)
                                        cgGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGAtt                                                                                                          <  4:391236/149‑1 (MQ=255)
                                         gggTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTc                                                                                                         <  3:142021/149‑1 (MQ=255)
                                             tttGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcacca                                                                                                     <  2:357693/149‑1 (MQ=255)
                                              ttGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                                                                    <  4:157060/149‑1 (MQ=255)
                                                  aCGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATg                                                                                                >  4:183293/1‑149 (MQ=255)
                                                    ggCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGc                                                                                                                                       <  4:465417/108‑1 (MQ=255)
                                                    ggCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGc                                                                                                                                       >  3:465417/1‑108 (MQ=255)
                                                         ttCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTcgc                                                                                         >  1:363030/1‑149 (MQ=255)
                                                              gaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCAcc                                                                                    >  4:354721/1‑149 (MQ=255)
                                                                   gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTc                                                                               <  4:200387/149‑1 (MQ=255)
                                                                   gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTc                                                                               <  1:153363/149‑1 (MQ=255)
                                                                              gcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCa                                                                    <  4:367717/149‑1 (MQ=255)
                                                                                gATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACt                                                                  <  1:85301/149‑1 (MQ=255)
                                                                                 aTAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACtt                                                                 <  4:121576/149‑1 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGc                                                                                                                         >  2:239564/1‑89 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGc                                                                                                                         <  1:239564/89‑1 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                             >  4:181561/1‑149 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                             >  1:65597/1‑149 (MQ=255)
                                                                                              aCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCAc                                                    <  3:415196/149‑1 (MQ=255)
                                                                                              aCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAAGCCCACTTTGACCTTCGCCAc                                                    <  2:264336/149‑1 (MQ=255)
                                                                                                 tATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTa                                                                            <  1:292988/122‑1 (MQ=255)
                                                                                                 tATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTa                                                                            >  2:292988/1‑122 (MQ=255)
                                                                                                  aTTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTtc                                                <  3:71259/149‑1 (MQ=255)
                                                                                                       gCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGcc                                           >  1:323297/1‑149 (MQ=255)
                                                                                                          aaCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTgg                                        <  1:453766/149‑1 (MQ=255)
                                                                                                             tgctgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCat                                     >  1:247319/1‑149 (MQ=255)
                                                                                                              gctgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCata                                    <  2:442132/149‑1 (MQ=255)
                                                                                                                tgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCACGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATc                                  <  4:251379/149‑1 (MQ=255)
                                                                                                                  aCCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTg                                >  1:210762/1‑149 (MQ=255)
                                                                                                                           cGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                       <  4:469106/149‑1 (MQ=255)
                                                                                                                                 ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                 <  4:471719/149‑1 (MQ=255)
                                                                                                                                 ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                 <  3:183293/149‑1 (MQ=255)
                                                                                                                                 ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                 <  2:247319/149‑1 (MQ=255)
                                                                                                                                    ggCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGAt              >  1:174755/1‑149 (MQ=255)
                                                                                                                                     gcgcGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATc             <  2:104970/149‑1 (MQ=255)
                                                                                                                                          aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg        >  3:264485/1‑149 (MQ=255)
                                                                                                                                          aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg        >  2:423158/1‑149 (MQ=255)
                                                                                                                                          aTTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg        >  4:150730/1‑149 (MQ=255)
                                                                                                                                             tGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGAt     <  1:423158/149‑1 (MQ=255)
                                                                                                                                                gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGc  <  2:363030/149‑1 (MQ=255)
                                                                                                                                                gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGc  <  4:468701/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGC  >  NC_000913/2375336‑2375628

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑