Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,577,355 C→A 6.5% A359A (GCC→GCA)  PP_4906 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,577,3550CA 6.5% 197.3 / 2.5 61A359A (GCC→GCA) PP_4906hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (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 = 8.63e-01

CCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCCTGCCCA  >  NC_002947/5577208‑5577500
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                                <  5:57836/138‑1 (MQ=255)
             gCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACTGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCaa                                                                                                                                                >  6:57836/1‑138 (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                                                                 <  5:291341/135‑1 (MQ=255)
                                                                                                ccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGcatca                                                                 >  6:291341/1‑135 (MQ=255)
                                                                                                ccAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                                   <  5:64338/149‑1 (MQ=255)
                                                                                                 cagcagCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACACAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                                   >  3:2838/1‑148 (MQ=255)
                                                                                                   gcagcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                              <  2:225239/103‑1 (MQ=255)
                                                                                                   gcagcGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGc                                                                                              >  1:225239/1‑103 (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                                                   >  1:374885/1‑134 (MQ=255)
                                                                                                               tGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTg                                                   <  2:374885/134‑1 (MQ=255)
                                                                                                               tGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGTCAAGGCCAGCGAc                                    >  5:149339/1‑149 (MQ=255)
                                                                                                                   aCGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTAcc                                >  7:100279/1‑149 (MQ=255)
                                                                                                                     gCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCAtg                                                               >  1:329297/1‑116 (MQ=255)
                                                                                                                     gCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCAtg                                                               <  2:329297/116‑1 (MQ=255)
                                                                                                                         gcgcACATTGTGCCAGCACCTGCTGGCCCAAACCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGcgcg                          >  1:40469/1‑149 (MQ=255)
                                                                                                                               aTTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGa                                                 <  8:81522/120‑1 (MQ=255)
                                                                                                                               aTTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGa                                                 >  7:81522/1‑120 (MQ=255)
                                                                                                                                   tgCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCaaa                 <  7:50481/148‑1 (MQ=255)
                                                                                                                                    gCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAAtt               >  3:219455/1‑149 (MQ=255)
                                                                                                                                      cAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCg                                      >  6:234499/1‑124 (MQ=255)
                                                                                                                                      cAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCg                                      <  5:234499/124‑1 (MQ=255)
                                                                                                                                         cACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTgcgc                           <  7:323465/132‑1 (MQ=255)
                                                                                                                                           ccTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGcc        >  2:233034/1‑149 (MQ=255)
                                                                                                                                                  gCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCCTGCCCa  >  2:368527/1‑148 (MQ=255)
                                                                                                                                                    ccAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCCTGCCCa  >  6:143899/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCACGCCCGGCACGCCCTGTTCCACCCGGCCGAGGCCTTGATCTGGCCGTGGCACC‑AACGCCGCCCTCACCCGGACATGCTGGCCCCGCCACCGCCCAGCAGCGACGACCTGGCACGCTGGCGCACATTGTGCCAGCACCTGCTGGCCCAACCCAGCCCGTTCACCAACGCCGTGCACCTGGCCACCCAGGCCCAAGAAGCCCTGCTGGCCTGCGGCATGCAACGCATCATGCTGCTCAGCCTGGACAAGGCCAGCGACTACCTGCGCGTGCAGCAAATTGCCGGCCTGCCCA  >  NC_002947/5577208‑5577500

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