Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,612,664 T→A 6.9% intergenic (+1626/‑481) hcpC‑II → / → PP_4084 hemolysin‑coregulated protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,612,6640TA 6.9% 210.8 / 5.1 72intergenic (+1626/‑481)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base A (2/3);  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 = 3.67e-01

ATGATCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACAC  >  NC_002947/4612516‑4612811
                                                                                                                                                    |                                                                                                                                                   
atgatCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAgt                                                                                                                                                     >  7:185914/1‑149 (MQ=255)
  gatCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACACACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAgtgt                                                                                                                                                   >  8:304179/1‑149 (MQ=255)
    tCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTc                                                                                                                                                 <  1:123559/149‑1 (MQ=255)
    tCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTc                                                                                                                                                 <  5:256142/149‑1 (MQ=255)
      aCCGGGCGCACTTACTGCGGTGACAAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCa                                                                                                                                                <  4:331803/148‑1 (MQ=255)
       ccGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATc                                                                                                                                              <  7:304179/149‑1 (MQ=255)
        cGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATc                                                                                                                                              >  1:217091/1‑148 (MQ=255)
             gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGTGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGcc                                                                                                                                        >  3:291782/1‑149 (MQ=255)
             gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               <  5:90220/142‑1 (MQ=255)
             gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               <  7:109103/142‑1 (MQ=255)
             gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               >  6:90220/1‑142 (MQ=255)
             gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               >  8:109103/1‑142 (MQ=255)
               aCTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCATGGATGAACACTTCCTCCTTGCCCa                                                                                                                                      <  5:151225/149‑1 (MQ=21)
               aCTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGcc                                                                                                                                        >  5:118229/1‑147 (MQ=255)
               aCTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCa                                                                                                                                      <  4:230219/149‑1 (MQ=255)
                cTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCag                                                                                                                                     <  6:118229/149‑1 (MQ=255)
                   aCTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCaga                                                                                                                                    <  1:147037/147‑1 (MQ=255)
                      gCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCaga                                                                                                                                    >  4:122412/1‑144 (MQ=255)
                       cGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATc                                                                                                                              <  3:248945/149‑1 (MQ=255)
                        ggTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATc                                                                                                                              <  2:189861/148‑1 (MQ=255)
                        ggTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATc                                                                                                                              >  1:189861/1‑148 (MQ=255)
                         gTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTGATCCATGCCCAGAGAGAt                                                                                                                               >  8:369119/1‑146 (MQ=255)
                              cAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAacac                                                                                                                       >  4:222387/1‑149 (MQ=255)
                                  tGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTg                                                                                                                   >  5:228789/1‑149 (MQ=255)
                                    ccgccTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGTGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGt                                                                                                                 <  4:291782/149‑1 (MQ=255)
                                    ccgccTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTg                                                                                                                   <  6:344646/147‑1 (MQ=255)
                                     cgccTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               >  8:152079/1‑118 (MQ=255)
                                     cgccTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                               <  7:152079/118‑1 (MQ=255)
                                         ttACGATCTGTCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAAt                                                                                                             >  1:364548/1‑148 (MQ=255)
                                         ttACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAAt                                                                                                             >  7:38526/1‑148 (MQ=255)
                                           aagATCTGCCCGATCACAAAACACGCATGAACATCAAGAGCAAGACCCACAAGGGGGAAGGTTTCAACGACCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATca                                                                                                           >  4:189466/3‑148 (MQ=255)
                                           aCGATCTGCCCGATCACATAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATca                                                                                                           >  2:68591/1‑148 (MQ=255)
                                             gATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACaa                                                                                                        >  6:134530/1‑149 (MQ=255)
                                                 tGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAg                                                                                                    >  4:272404/1‑149 (MQ=255)
                                                   cccGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATc                                                                                                            <  2:248383/139‑1 (MQ=255)
                                                        tCAAAAAACACTCATGACCATCAAGATCCAGACCCACAAGTGGGAATGATTCAACGAGCTGCGTTTCGAGGATTAGATGTGCAAGGACGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCct                                                                                               <  5:274542/147‑1 (MQ=255)
                                                          acaAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCctctct                                                                                           <  3:122412/149‑1 (MQ=255)
                                                          acaAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCctctct                                                                                           <  3:198680/149‑1 (MQ=255)
                                                            aaaaCACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCtgt                                                                                         >  3:327508/1‑149 (MQ=255)
                                                             aaaCACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCtgtg                                                                                        >  7:201014/1‑149 (MQ=255)
                                                              aaCACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTgg                                                                                       <  8:201014/149‑1 (MQ=255)
                                                              aaCACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTgg                                                                                       <  8:314052/149‑1 (MQ=255)
                                                                   gCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTc                                                                                  >  2:214006/1‑149 (MQ=255)
                                                                     aTGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGa                                                                                >  7:204348/1‑149 (MQ=255)
                                                                      tGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGaa                                                                                                                          <  5:244194/106‑1 (MQ=255)
                                                                      tGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGaa                                                                                                                          >  6:244194/1‑106 (MQ=255)
                                                                         ccATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCGCAACCAGTCCCTCTCTGTGGGCGTCGATCg                                                                             >  4:86451/1‑148 (MQ=255)
                                                                                gagCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAAc                                                                     <  8:185914/149‑1 (MQ=255)
                                                                                  gCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGAACCTCTCTGTGGGCGTCGATCGCCCacaacat                                                                   >  4:21819/1‑148 (MQ=255)
                                                                                     aGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGtt                                                                >  5:334231/1‑149 (MQ=255)
                                                                                      gACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTGATCCATGCCCaga                                                                                                                                    >  6:274883/1‑80 (MQ=255)
                                                                                          cacaAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCTGAACACTCTGGTCATTCACAACCAGTCGCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGc                                                             >  2:274219/1‑147 (MQ=255)
                                                                                            caAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTTGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAgg                                                         <  3:295666/149‑1 (MQ=255)
                                                                                             aaGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGa                                                        <  6:170678/149‑1 (MQ=255)
                                                                                                 gggAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTATGGTCAAACACACCCAGACCGTATATGTGGCGGTCGAACGCACCCAACAGGTTGGCCAGGCCGa                                                     >  2:126829/1‑148 (MQ=255)
                                                                                                    aaGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGATGGCCAGGACGAATCg                                                 <  6:228789/149‑1 (MQ=255)
                                                                                                          tCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGc                                            <  3:222387/148‑1 (MQ=255)
                                                                                                           cAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGcc                                           <  6:334231/148‑1 (MQ=255)
                                                                                                               gttcttggTTTCGATGTTGATCTGGGCAAGGAGGAAGTGTTCATCCATGCCCTTAGAGATCAGAACACTCTCGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTgg                                      <  1:126829/142‑1 (MQ=255)
                                                                                                                 gCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGcc                                    <  5:134530/149‑1 (MQ=255)
                                                                                                                  gtGCGGTTCGCGGATGAGCTGGGCAAGGAGGAGGTGTTGGTCCATGGCCAGAGTGATCAGAACCCTCTGGTCAGTCACGACGAGTGCCTCTGTGTGGGCGTCGCTCGCAGCCAGCGGGTTGGGCAGGACGACTCGGTCGCCATGGGGCg                                   >  4:48413/2‑149 (MQ=255)
                                                                                                                  cTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCg                                   <  3:189466/149‑1 (MQ=255)
                                                                                                                             ggATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTgc                        <  1:68591/149‑1 (MQ=255)
                                                                                                                                ctctcTGGGCATGGATGAACACTTCCTCCTTGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGcgcg                     >  1:17248/5‑149 (MQ=38)
                                                                                                                                  aGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGcgtc                   >  7:270369/1‑149 (MQ=255)
                                                                                                                                   tcTGGGCATGGATGAACACTTCCTCCTTGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTc                                              >  3:331246/2‑121 (MQ=38)
                                                                                                                                   tcTGGGCATGGATGAACACTTCCTCCTTGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTc                                              <  4:331246/120‑1 (MQ=38)
                                                                                                                                   tcTGGGCATGGATGAACACTTCCTCCTTGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGcgcg                     <  2:17248/145‑1 (MQ=255)
                                                                                                                                      gggCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCgtcgtc                >  4:117408/1‑148 (MQ=255)
                                                                                                                                       ggCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCgtcgtc                >  7:127753/1‑147 (MQ=255)
                                                                                                                                          aaGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGc           <  6:264340/149‑1 (MQ=255)
                                                                                                                                              aggaAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAAc       <  7:369119/149‑1 (MQ=255)
                                                                                                                                               ggaAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAAc       >  8:124247/1‑148 (MQ=255)
                                                                                                                                                  aGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGaca   <  4:327508/149‑1 (MQ=255)
                                                                                                                                                   gtgtTCATCCATGCCCAGAGAGATCAGAACACTCTGGGCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGacac  <  3:48413/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGATCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACAC  >  NC_002947/4612516‑4612811

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