Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,215,092 G→C 6.7% L623V (CTC→GTC)  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,0920GC 6.7% 194.6 / 4.4 60L623V (CTC→GTC) PP_1061ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

AGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACTGCGCC  >  NC_002947/1214944‑1215238
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                 <  2:36842/132‑1 (MQ=255)
                                                                    aTCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGAAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGc                                                                                                 >  1:36842/1‑132 (MQ=255)
                                                                     tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                           >  3:319248/1‑121 (MQ=255)
                                                                     tCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                           <  4:319248/121‑1 (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                                                                                                           <  8:90378/106‑1 (MQ=255)
                                                                                    aTCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCt                                                                                                           >  7:90378/1‑106 (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                                                                       >  7:195191/1‑93 (MQ=255)
                                                                                                                                     gggCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTg                                                                       <  8:195191/93‑1 (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)
                                                                                                                                                  gAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACTGCGcc  >  4:205832/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGCCCCAGGCGCGGTTGATGCGGCTACCGTAGGGTGAATGAATGATCAGTTGGGTCCCGCCGGACTCATCGAAAAAACGTTCCATCACCAGCGTATCCTGCGAGGGCAGGGCGCCAAGCACTTCACGGGTGCGGGCCAGGTAGTCGAGCAACTGGCCCGCACAAGCTTCGTCCAGTTCGAAGGTCGCCTGCAACCAGGCCAGTACCTCCACCGTGTCGCCGCCTGCCAGCGCCAGTTGCTGGTCGATTCGCGCCTGCAGCCGGGCCACGGCGGCAGACAGCTCGTCACTGCGCC  >  NC_002947/1214944‑1215238

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: ‑