Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,420,452 G→C 5.5% A216G (GCC→GGC)  PP_4760 ← alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,420,4520GC 5.5% 234.3 / 2.9 72A216G (GCC→GGC) PP_4760alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base C (2/2);  total (35/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 = 9.01e-01

GGCCACCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGACACACCGCCGG  >  NC_002947/5420304‑5420593
                                                                                                                                                    |                                                                                                                                             
ggCCACCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAgg                                                                                                                                               >  6:5915/1‑149 (MQ=255)
 gCCACCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGc                                                                                                                                              >  8:132966/1‑149 (MQ=255)
   cACCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCg                                                                                                                                            <  6:375057/149‑1 (MQ=255)
    aCCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCgg                                                                                                                                           >  1:356334/1‑149 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          <  3:208084/149‑1 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          <  3:119263/149‑1 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          <  2:97581/149‑1 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          <  6:236300/149‑1 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          >  2:89796/1‑149 (MQ=255)
     ccccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCggg                                                                                                                                          <  5:5915/149‑1 (MQ=255)
      cccGGTAAGGATGCCGATCAGCGAAACATGCCCGCCTATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGAACGACACCACCGTTGCCAGTCAGCAAGCGCACCTTCTCACCCCAGGCCcggg                                                                                                                                         >  1:60333/1‑149 (MQ=255)
             aGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGctct                                                                                                                                   <  4:112349/148‑1 (MQ=255)
             aGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCtt                                                                                                                                  <  5:194381/149‑1 (MQ=255)
             aGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCtt                                                                                                                                  >  8:22211/1‑149 (MQ=255)
             aGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCtt                                                                                                                                  <  6:334793/149‑1 (MQ=255)
                  gCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGt                                                                                                                             >  7:198051/1‑149 (MQ=255)
                   ccGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGt                                                                                                                             <  1:374573/148‑1 (MQ=255)
                    cgttCAGCGTAACTTGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCGACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTGTCGCAGCTTCTCAGCGCAGGCCGGGGTGCTGTTGTAGTTg                                                                                                                           <  5:1840/146‑1 (MQ=255)
                       tCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATc                                                                                                                        <  5:27815/149‑1 (MQ=255)
                         aGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCa                                                                                                                       <  8:55477/148‑1 (MQ=255)
                          gCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGCCACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAgg                                                                                                                     >  8:45551/1‑149 (MQ=255)
                          gCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCGGCCTGGGGTgagaa                                                                                                                                   <  5:145932/135‑5 (MQ=255)
                          gCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCGGCCTGGGGTgagaa                                                                                                                                   >  6:145932/1‑131 (MQ=255)
                           cGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGctc                                                                                                                                    <  2:216256/133‑1 (MQ=255)
                           cGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGctc                                                                                                                                    >  1:216256/1‑133 (MQ=255)
                           cGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGAGCACCTTCTCCCCCCAGGCCGGGGGGCTCTTGTAGTTGATCAggt                                                                                                                    >  8:317417/1‑149 (MQ=255)
                             aaaCATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAggtgg                                                                                                                  >  7:289286/1‑149 (MQ=255)
                                cATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAggt                                                                                                                    >  6:343030/1‑144 (MQ=255)
                                     ccGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATAGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCAcc                                                                                                          >  8:319498/1‑149 (MQ=255)
                                         cAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAgg                                                                                                      <  5:380747/149‑1 (MQ=255)
                                               gcgcAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCtt                                                                                                 >  5:113291/1‑148 (MQ=255)
                                               gcgcAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTc                                                                                                <  3:341382/149‑1 (MQ=255)
                                               gcgcAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTc                                                                                                <  8:289286/149‑1 (MQ=255)
                                               gcgcAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCGCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCtt                                                                                                 >  4:126454/1‑148 (MQ=255)
                                                  cagcaATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAgg                                                                                             <  7:22211/149‑1 (MQ=255)
                                                     caATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGTTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGgcgc                                                                                          <  1:81126/149‑1 (MQ=255)
                                                        tcatcGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCGGCCTGGGGTgagaag                                                                                                                                  >  8:118370/1‑101 (MQ=255)
                                                        tcatcGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCGGCCTGGGGTgagaag                                                                                                                                  <  7:118370/106‑6 (MQ=255)
                                                        tcatcGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTc                                                                                        >  6:300482/1‑148 (MQ=255)
                                                         catcGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCa                                                                                                            >  4:221679/1‑127 (MQ=255)
                                                         catcGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCa                                                                                                            <  3:221679/127‑1 (MQ=255)
                                                             gACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCt                                                                                   <  1:89796/148‑1 (MQ=255)
                                                                  cTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGc                                                                               <  2:185445/147‑1 (MQ=255)
                                                                     cAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGc                                                                          >  2:242999/1‑149 (MQ=255)
                                                                     cAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGc                                                                          >  7:42/1‑149 (MQ=255)
                                                                                    gccCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGa                                                           <  6:284708/149‑1 (MQ=255)
                                                                                         cTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACg                                                      >  4:301136/1‑149 (MQ=255)
                                                                                             aaTCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTc                                                   <  1:242999/148‑1 (MQ=255)
                                                                                              aTCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGTATCGCTGGACGAGGTGGCGATCACGGTcgc                                                 <  8:198051/149‑1 (MQ=255)
                                                                                                   aTGATCGACACAACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGcgcccg                                            >  2:23091/1‑149 (MQ=255)
                                                                                                      aTCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCATGGTCGCGCCCGCCa                                         <  5:300482/149‑1 (MQ=255)
                                                                                                      aTCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGcgcccgc                                           >  2:18564/1‑147 (MQ=255)
                                                                                                       tCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAg                                        <  2:60333/149‑1 (MQ=255)
                                                                                                           caccacGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTg                                    <  2:267475/149‑1 (MQ=255)
                                                                                                                 ggTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAAc                              <  7:132966/149‑1 (MQ=255)
                                                                                                                     gCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGc                           <  7:45551/148‑1 (MQ=255)
                                                                                                                       cAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAgg                        <  6:113291/149‑1 (MQ=255)
                                                                                                                       cAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAgg                        >  5:190376/1‑149 (MQ=255)
                                                                                                                       cAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAgg                        >  4:228837/1‑149 (MQ=255)
                                                                                                                       cAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAgg                        >  3:238710/1‑149 (MQ=255)
                                                                                                                         gTCAGCTCGCGCACCTTCTCACCCCAGGCCCGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGc                      >  5:198982/1‑149 (MQ=255)
                                                                                                                            aGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGa                    >  8:77684/1‑148 (MQ=255)
                                                                                                                            aGCTCGCGCAACTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCGTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGa                    >  6:120921/1‑148 (MQ=255)
                                                                                                                             gCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGaa                   >  4:260165/1‑148 (MQ=255)
                                                                                                                                 gcgcACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGGCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATc               <  5:78178/148‑1 (MQ=255)
                                                                                                                                  cgcACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGa             <  7:317417/149‑1 (MQ=255)
                                                                                                                                  cgcACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGa             <  3:228837/149‑1 (MQ=255)
                                                                                                                                   gcACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGa             >  1:264845/1‑148 (MQ=255)
                                                                                                                                     aCCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCAACCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGaca           >  6:139232/1‑148 (MQ=255)
                                                                                                                                        ttCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGACACAcc       >  6:147128/1‑149 (MQ=255)
                                                                                                                                         tctcACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGACACAccg      >  7:217143/1‑149 (MQ=255)
                                                                                                                                             aCCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGACACAccgccg   <  3:260165/148‑1 (MQ=255)
                                                                                                                                             aCCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGCGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGAGAACTGCAGGGCGAAAATCGACACACCGCCgg  >  1:332379/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGCCACCCCGGTAAGGATGCCGATCAGCGAAACATGCCCGCCAATGCGCGCAGCAATCATCGACTGCTCCAGCGTCGCCGGGCCGCCCACTTCAATCACATGATCGACACCACGGTTGCCAGTCAGCTCGCGCACCTTCTCACCCCAGGCCGGGGTGCTCTTGTAGTTGATCAGGTGGTCGGCACCCAGGGCCTTCAGGCGCTCCAGCTTGGCATCGCTGGACGAGGTGGCGATCACGGTCGCGCCCGCCAGCTTGGCGAACTGCAGGGCGAAAATCGACACACCGCCGG  >  NC_002947/5420304‑5420593

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