Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,577,346 C→G 6.1% H356Q (CAC→CAG)  PP_4906 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,577,3460CG 6.1% 213.5 / 3.9 66H356Q (CAC→CAG) PP_4906hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GCGCCCGCCGCCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCC  >  NC_002947/5577198‑5577494
                                                                                                                                                     |                                                                                                                                                    
gcgcCCGCCGCCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑ACCGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCCGCAGCGACGCCCTGGCACGCTGGCGCACATTGTGCCAGCAc                                                                                                                                                      >  7:100814/1‑149 (MQ=255)
  gcccgccGCCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACct                                                                                                                                                    <  2:62792/149‑1 (MQ=255)
  gcacgccgcCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACct                                                                                                                                                    >  7:311932/4‑149 (MQ=255)
   cccgccGCCACGCCCGGCCCGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACCTGCTGGCCCCGCCCCCGCCCCGCAGCGACGACCTGGCACCCGGGCGCCCATTGTGCCAGCACctg                                                                                                                                                   >  5:82620/1‑149 (MQ=255)
   cccgccGCCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACctg                                                                                                                                                   <  7:275211/149‑1 (MQ=255)
       ccgccACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTgg                                                                                                                                               >  2:187818/1‑149 (MQ=255)
          ccACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGccc                                                                                                                                            <  8:109587/149‑1 (MQ=255)
            aCGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACCAAACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCCACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGccc                                                                                                                                            >  1:103686/1‑148 (MQ=255)
              gCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAAcc                                                                                                                                        <  1:187818/149‑1 (MQ=255)
                       gCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACTGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCaa                                                                                                                                          >  6:57836/1‑138 (MQ=255)
                       gCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACTGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCaa                                                                                                                                          <  5:57836/138‑1 (MQ=255)
                         ccTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCa                                                                                                                             >  8:2622/1‑149 (MQ=255)
                          cTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCAc                                                                                                                            <  8:311932/149‑1 (MQ=255)
                                 aCCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGcc                                                                                                                     >  7:77343/1‑149 (MQ=255)
                                 aCCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGGGCCAGCAGGTGCTGGCACaat                                                                                                                                         <  7:215254/129‑2 (MQ=255)
                                 aCCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGGGCCAGCAGGTGCTGGCACaat                                                                                                                                         >  8:215254/1‑128 (MQ=255)
                                   ccGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGt                                                                                                                   >  8:50481/1‑149 (MQ=255)
                                    cGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTg                                                                                                                  <  2:103686/149‑1 (MQ=255)
                                       ccGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCAc                                                                                                               <  4:238751/149‑1 (MQ=255)
                                        cGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCGGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCAcc                                                                                                              <  7:18612/149‑1 (MQ=255)
                                         gAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGtt                                                                                                                               <  2:196404/131‑1 (MQ=255)
                                         gAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGtt                                                                                                                               >  1:196404/1‑131 (MQ=255)
                                                  aTCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCaac                                                                                                      >  2:165288/1‑145 (MQ=255)
                                                    cTGGCCGTGGCACC‑AAGGCCGCCCTCACCCGGACATGCTGGCCTCACCACCGCCAAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCTGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACTCCGTGCACCTGGCCACCCAg                                                                                                   <  8:100814/148‑1 (MQ=255)
                                                      ggCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGcc                                                                                                >  7:93993/1‑149 (MQ=255)
                                                       gCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGccc                                                                                               <  8:294290/149‑1 (MQ=255)
                                                              cACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGGCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                        <  6:82620/149‑1 (MQ=255)
                                                              cACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                        >  3:311180/1‑149 (MQ=255)
                                                               aCC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGcc                                                                                       <  1:170799/149‑1 (MQ=255)
                                                               aCC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGGGCCAGCAGGTGCTGGCACaa                                                                                                                                          <  8:297864/98‑1 (MQ=255)
                                                               aCC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGGGCCAGCAGGTGCTGGCACaa                                                                                                                                          >  7:297864/1‑98 (MQ=255)
                                                                 c‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCct                                                                                     >  5:248218/1‑149 (MQ=255)
                                                                      gccgccCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCctgctg                                                                                 <  7:2622/149‑1 (MQ=255)
                                                                          cccTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCt                                                                             <  5:256659/149‑1 (MQ=255)
                                                                           ccTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTCGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTg                                                                            <  1:188221/149‑1 (MQ=255)
                                                                                 ccGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCAt                                                                      <  8:77343/149‑1 (MQ=255)
                                                                                  cGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATg                                                                     <  5:231853/149‑1 (MQ=255)
                                                                                   ggACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTGACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATg                                                                     >  7:355378/1‑148 (MQ=255)
                                                                                          cTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGcat                                                             >  3:94740/1‑149 (MQ=255)
                                                                                          cTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGca                                                              >  6:275055/1‑148 (MQ=255)
                                                                                                  gCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCAtgctgc                                                     >  4:110355/1‑149 (MQ=255)
                                                                                                         cccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCt                                              <  1:157370/149‑1 (MQ=255)
                                                                                                          ccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGcatca                                                           >  6:291341/1‑135 (MQ=255)
                                                                                                          ccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGcatca                                                           <  5:291341/135‑1 (MQ=255)
                                                                                                          ccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                             <  5:64338/149‑1 (MQ=255)
                                                                                                           cagcagCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACACAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                             >  3:2838/1‑148 (MQ=255)
                                                                                                             gcagcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                        >  1:225239/1‑103 (MQ=255)
                                                                                                             gcagcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                        <  2:225239/103‑1 (MQ=255)
                                                                                                                gcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCAc                                                                                                       <  6:229204/85‑1 (MQ=255)
                                                                                                                gcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCAc                                                                                                       >  5:229204/1‑85 (MQ=255)
                                                                                                                  gacgacCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTgg                                                                                                           >  5:3668/1‑79 (MQ=255)
                                                                                                                  gacgacCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTgg                                                                                                           <  6:3668/79‑1 (MQ=255)
                                                                                                                    cgacCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCa                                   <  4:311180/149‑1 (MQ=255)
                                                                                                                         tGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                             <  2:374885/134‑1 (MQ=255)
                                                                                                                         tGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                             >  1:374885/1‑134 (MQ=255)
                                                                                                                         tGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGTCAAGGCCAGCGAc                              >  5:149339/1‑149 (MQ=255)
                                                                                                                             aCGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTAcc                          >  7:100279/1‑149 (MQ=255)
                                                                                                                               gCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCAtg                                                         <  2:329297/116‑1 (MQ=255)
                                                                                                                               gCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCAtg                                                         >  1:329297/1‑116 (MQ=255)
                                                                                                                                   gcgcACATTGTGCCAGCACCTGCTGGCCCAAACCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGcgcg                    >  1:40469/1‑149 (MQ=255)
                                                                                                                                         aTTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGa                                           >  7:81522/1‑120 (MQ=255)
                                                                                                                                         aTTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGa                                           <  8:81522/120‑1 (MQ=255)
                                                                                                                                             tgCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCaaa           <  7:50481/148‑1 (MQ=255)
                                                                                                                                              gCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAAtt         >  3:219455/1‑149 (MQ=255)
                                                                                                                                                cAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCg                                <  5:234499/124‑1 (MQ=255)
                                                                                                                                                cAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCg                                >  6:234499/1‑124 (MQ=255)
                                                                                                                                                   cACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTgcgc                     <  7:323465/132‑1 (MQ=255)
                                                                                                                                                     ccTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGcc  >  2:233034/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCCCGCCGCCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCC  >  NC_002947/5577198‑5577494

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