Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,215,086 A→G 6.5% Y625H (TAC→CAC) ‡ PP_1061 ← ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,215,0860AG 6.5% 180.5 / 3.0 61Y625H (TAC→CAC) ‡PP_1061ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

GGCCAGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACT  >  NC_002947/1214940‑1215233
                                                                                                                                                  |                                                                                                                                                   
ggCCAGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGt                                                                                                                                                   <  5:338299/149‑1 (MQ=255)
    aGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAg                                                                                                                                               <  2:355966/149‑1 (MQ=255)
      gCCACAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGc                                                                                                                                              >  4:111309/1‑148 (MQ=255)
       ccccAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCaa                                                                                                                                            <  1:148031/149‑1 (MQ=255)
        cccAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAAc                                                                                                                                           >  4:114803/1‑149 (MQ=255)
              cgcgGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCTGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGccc                                                                                                                                     <  1:82443/149‑1 (MQ=255)
              cgcgGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGTTGCTCGACTACCTGGccc                                                                                                                                     <  3:227877/149‑1 (MQ=255)
                   ttGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGcac                                                                                                                                 >  4:345717/1‑148 (MQ=255)
                                     ggTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGTTGCTCGACTACCTGGCCCGcacccgtg                                                                                                                            >  4:111239/1‑130 (MQ=38)
                                     ggTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGTTGCTCGACTACCTGGCCCGcacccgtg                                                                                                                            <  3:111239/135‑6 (MQ=38)
                                        gaatgaatGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAgg                                                                                                            >  2:74751/1‑148 (MQ=255)
                                           tgaatgaTCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGc                                                                                                        <  3:111309/149‑1 (MQ=255)
                                            gaatgaTCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGcc                                                                                                       >  1:188447/1‑149 (MQ=255)
                                                      ttGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGtt                                                                                                                  >  2:393896/1‑128 (MQ=255)
                                                      ttGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGtt                                                                                                                  <  1:393896/128‑1 (MQ=255)
                                                           tCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGt                                                                                        <  1:74751/149‑1 (MQ=255)
                                                             ccgccgGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTAc                                                                                      >  8:73522/1‑149 (MQ=255)
                                                                    aCTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCAcc                                                                               <  2:188447/149‑1 (MQ=255)
                                                                        aTCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCgtgt                                                                           >  4:344751/1‑149 (MQ=255)
                                                                        aTCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCgt                                                                             <  2:59158/147‑1 (MQ=255)
                                                                        aTCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGAAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGc                                                                                            >  1:36842/1‑132 (MQ=255)
                                                                        aTCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGAAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGc                                                                                            <  2:36842/132‑1 (MQ=255)
                                                                         tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                      <  4:319248/121‑1 (MQ=255)
                                                                         tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                      >  3:319248/1‑121 (MQ=255)
                                                                         tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTc                                                                          <  3:344751/149‑1 (MQ=255)
                                                                         tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTc                                                                          <  3:345717/149‑1 (MQ=255)
                                                                          cGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTcg                                                                         >  2:173313/1‑149 (MQ=255)
                                                                           gAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTcgc                                                                        >  8:95873/1‑149 (MQ=255)
                                                                                aaCGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCt                                                                   >  8:191188/1‑149 (MQ=255)
                                                                                  cggtCCAGCACCATCGTATCCAGCGAGGGCAGGGCTCCATGCTCTTCCCGGGTGCGGGCCAGGTAGTCGAGCATCTGGCCCGCACTAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTAGCTCCACCGTGTCGCCGCCTGc                                                                 <  8:165667/146‑1 (MQ=255)
                                                                                    ttCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGc                                                                 >  5:381473/1‑147 (MQ=255)
                                                                                    ttCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCa                                                               <  2:31086/149‑1 (MQ=255)
                                                                                      ccctCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGGCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAgc                                                             >  6:291610/4‑149 (MQ=255)
                                                                                       cATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAgcg                                                            >  6:402452/1‑149 (MQ=255)
                                                                                        aTCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                      >  7:90378/1‑106 (MQ=255)
                                                                                        aTCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                      <  8:90378/106‑1 (MQ=255)
                                                                                             cAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGtt                                                      >  7:6413/1‑149 (MQ=255)
                                                                                                   aTCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGt                                                >  3:264549/1‑149 (MQ=255)
                                                                                                    tCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGt                                                <  4:264549/148‑1 (MQ=255)
                                                                                                      cTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGa                                             <  7:95873/149‑1 (MQ=255)
                                                                                                       tGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGAt                                            >  7:232517/1‑149 (MQ=255)
                                                                                                        gCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGAtt                                           >  2:40105/1‑149 (MQ=255)
                                                                                                           agggcagggcGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGtt                                                      <  8:328487/135‑1 (MQ=255)
                                                                                                           agggcagggcGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGtt                                                      >  7:328487/1‑135 (MQ=255)
                                                                                                               cagggcGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCTCCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCt                                    <  7:73522/149‑1 (MQ=255)
                                                                                                                                cACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCAcggc                    <  8:6413/148‑1 (MQ=255)
                                                                                                                                 aCGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGa                                             <  7:189980/122‑1 (MQ=255)
                                                                                                                                 aCGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGa                                             >  8:189980/1‑122 (MQ=255)
                                                                                                                                 aCGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCAcggcg                   <  1:40105/148‑1 (MQ=255)
                                                                                                                                    tgTGCGGGCCAGTTGCTCGACTACCTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCAcggcgg                  >  6:397686/2‑146 (MQ=39)
                                                                                                                                     gTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCCGCCGTGGCCCGGCTGCAgg              <  6:381473/149‑2 (MQ=255)
                                                                                                                                         gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTg                                                                                                     >  4:149584/1‑58 (MQ=255)
                                                                                                                                         gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTg                                                                                                     <  3:149584/58‑1 (MQ=255)
                                                                                                                                         gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTg                                                                  <  8:195191/93‑1 (MQ=255)
                                                                                                                                         gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTg                                                                  >  7:195191/1‑93 (MQ=255)
                                                                                                                                         gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGc          >  5:356888/1‑149 (MQ=255)
                                                                                                                                          ggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCt         >  6:175906/1‑149 (MQ=255)
                                                                                                                                           gCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTc        <  1:173313/149‑1 (MQ=255)
                                                                                                                                           gCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCACGCCTGCAGCCGGGCCACGGCGGCAGACAGCTc        <  7:191188/149‑1 (MQ=255)
                                                                                                                                               ggTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCa    <  5:291610/149‑1 (MQ=255)
                                                                                                                                                 tAGTCGAGCAACTGGCCCGCACAAGCTTTGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACt  >  2:40070/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCCAGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACT  >  NC_002947/1214940‑1215233

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