Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,627 C→G 7.4% P318P (CCG→CCC)  flhB ← flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,6270CG 7.4% 210.4 / 6.8 68P318P (CCG→CCC) flhBflagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base G (3/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-02

GTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCG  >  NC_002947/4945485‑4945769
                                                                                                                                              |                                                                                                                                              
gTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAt                                                                                                                                          <  6:153603/149‑1 (MQ=255)
   gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTc                                                                                                                                        <  1:241330/148‑1 (MQ=255)
   gCTCT‑GGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCg                                                                                                                                       <  6:267546/148‑1 (MQ=255)
      catggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAgg                                                                                                                                    <  1:114544/146‑1 (MQ=255)
       tGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                     <  7:160648/147‑1 (MQ=255)
         ggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCACAAGTTCGGTGGAGTAATAGATAGCCCGCCCCAGCGCCGGCGATTCGAGGAt                                                                                                                                  >  2:278248/1‑148 (MQ=255)
                  ccttgccGGCGCGGTATTGGCGAATCTGGACCACATACGCCAGCACCTGCGCCACCGCCATGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtggtg                                                                                                                        <  3:220234/149‑1 (MQ=255)
                  ccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgttg                                                                                                                        <  7:288838/149‑1 (MQ=255)
                         ggCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAgg                                                                                                                                    >  7:87879/1‑130 (MQ=255)
                         ggCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAgg                                                                                                                                    <  8:87879/130‑1 (MQ=255)
                         ggCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCa                                                                                                                 >  4:106854/1‑149 (MQ=255)
                            gcgGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTgg                                                                                                                              >  4:183604/1‑133 (MQ=255)
                            gcgGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTgg                                                                                                                              <  3:183604/133‑1 (MQ=255)
                                     gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgt                                                                                                                          >  1:146472/1‑128 (MQ=255)
                                     gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgt                                                                                                                          <  2:146472/128‑1 (MQ=255)
                                     gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAt                                                                                                                            >  4:289687/1‑126 (MQ=255)
                                     gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAt                                                                                                                            <  3:289687/126‑1 (MQ=255)
                                     gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCAc                                                                                                     >  3:253362/1‑149 (MQ=255)
                                     gCGAATCTCGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGTCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                      <  8:167273/148‑1 (MQ=255)
                                       gAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                       <  1:14303/113‑4 (MQ=255)
                                       gAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                       >  2:14303/1‑110 (MQ=255)
                                       gAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACg                                                                                                    >  3:50583/1‑148 (MQ=255)
                                          tCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATc                                                                                                >  2:74075/1‑149 (MQ=255)
                                             ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCtt                                                                                              >  5:155511/1‑148 (MQ=255)
                                              gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGa                                                                                                   <  8:261666/142‑1 (MQ=255)
                                              gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGa                                                                                                   >  7:261666/1‑142 (MQ=255)
                                              gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                             <  6:155511/148‑1 (MQ=255)
                                                acacATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAAc                                                                                          >  2:143399/1‑149 (MQ=255)
                                                        gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                             <  2:97898/138‑1 (MQ=255)
                                                        gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                             >  1:97898/1‑138 (MQ=255)
                                                        gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGa                                                                                   >  5:362259/1‑148 (MQ=255)
                                                         ccAGCTCCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCTGTTGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGata                                                                                 <  1:278248/149‑3 (MQ=255)
                                                          cAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCCCGCCAAGTTCACGAATCGTCAACGCCATGAAAt                                                                                >  8:121864/1‑149 (MQ=255)
                                                           aGCACCTGCGCCACCGACAGGTACAGCCCAGCCGTGATATCGTGCTCTAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGTTCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATc                                                                               <  5:325728/149‑1 (MQ=255)
                                                                  gcgcCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAccc                                                                        >  4:370457/1‑149 (MQ=255)
                                                                        ccGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGcc                                                                  <  1:143399/149‑1 (MQ=255)
                                                                        ccGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGcc                                                                  <  3:370457/149‑1 (MQ=255)
                                                                            cAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagca                                                              <  2:155982/149‑1 (MQ=255)
                                                                               gTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATAGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCagcagc                                                            >  3:178873/1‑148 (MQ=255)
                                                                                 aCAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAAc                                                         >  4:30269/1‑149 (MQ=255)
                                                                                        cgccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGCTCTGGATGTTGTGCTCCACGCCAATTTCACTAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGcgccg                                                  >  8:313453/1‑149 (MQ=255)
                                                                                                  tcCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCttt                                        >  5:47359/1‑149 (MQ=255)
                                                                                                    cTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCt                                          <  4:50583/145‑1 (MQ=255)
                                                                                                      gCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTcc                                    >  2:374386/1‑149 (MQ=255)
                                                                                                        tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCGACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                  >  7:149825/1‑149 (MQ=255)
                                                                                                        tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                  >  5:220316/1‑149 (MQ=255)
                                                                                                          aaGTTCGGTGGAGTAATAGCTCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGt                                <  6:362259/149‑1 (MQ=255)
                                                                                                          aaGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                  <  6:220316/147‑1 (MQ=255)
                                                                                                                    gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTc                               <  8:320357/140‑1 (MQ=255)
                                                                                                                    gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTc                               >  7:320357/1‑140 (MQ=255)
                                                                                                                    gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCa                       >  1:273004/1‑148 (MQ=255)
                                                                                                                      gTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAg                    >  7:148165/1‑149 (MQ=255)
                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAATTCGCTACCCTTGGGCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                 >  7:201863/1‑149 (MQ=255)
                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagcag                                                             >  8:148778/1‑105 (MQ=255)
                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagcag                                                             <  7:148778/105‑1 (MQ=255)
                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                 >  1:164857/1‑149 (MQ=255)
                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                 <  2:273004/149‑1 (MQ=255)
                                                                                                                             aTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCg              >  8:358521/1‑148 (MQ=255)
                                                                                                                              tCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGc             >  4:7740/1‑148 (MQ=255)
                                                                                                                                 cccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATaa         <  1:74075/149‑1 (MQ=255)
                                                                                                                                 cccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATaa         >  3:5438/1‑149 (MQ=255)
                                                                                                                                  cggcgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTTCTCCAAGCTAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCTACTTTCTCCGGTTCGTATTGCAAAGCCACCTCATAAt        <  4:178873/147‑1 (MQ=255)
                                                                                                                                  ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCGACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAAt        <  8:149825/149‑1 (MQ=255)
                                                                                                                                  ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAAt        <  6:47359/149‑1 (MQ=255)
                                                                                                                                     cgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTg                         <  3:179849/129‑1 (MQ=255)
                                                                                                                                     cgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTg                         >  4:179849/1‑129 (MQ=255)
                                                                                                                                        cAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCg  >  4:322359/1‑149 (MQ=255)
                                                                                                                                           cgcCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTgg                                                                    <  2:208752/80‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCG  >  NC_002947/4945485‑4945769

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