Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,612,659 G→T 7.1% intergenic (+1621/‑486) 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,6590GT 7.1% 220.5 / 6.4 71intergenic (+1621/‑486)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base T (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

GCCGATGATCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAAC  >  NC_002947/4612512‑4612806
                                                                                                                                                   |                                                                                                                                                   
gCCGATGATCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAaggag                                                                                                                                                     >  3:230219/1‑148 (MQ=255)
    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                                                                                                                                          >  6:90220/1‑142 (MQ=255)
                 gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                          <  5:90220/142‑1 (MQ=255)
                 gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                          <  7:109103/142‑1 (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                                                                                                                                          <  7:152079/118‑1 (MQ=255)
                                         cgccTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCAt                                                                                                                                          >  8:152079/1‑118 (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:198680/149‑1 (MQ=255)
                                                              acaAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCctctct                                                                                      <  3:122412/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                                                                                                                     >  6:244194/1‑106 (MQ=255)
                                                                          tGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGaa                                                                                                                     <  5:244194/106‑1 (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                                         <  4:331246/120‑1 (MQ=38)
                                                                                                                                       tcTGGGCATGGATGAACACTTCCTCCTTGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTc                                         >  3:331246/2‑121 (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)
                                                                                                                                                   |                                                                                                                                                   
GCCGATGATCACCGGGCGCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAAC  >  NC_002947/4612512‑4612806

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: