Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,612,673 A→T 7.4% intergenic (+1635/‑472) 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,6730AT 7.4% 191.8 / 5.1 68intergenic (+1635/‑472)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base T (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

GCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACACCTTGA  >  NC_002947/4612529‑4612816
                                                                                                                                                |                                                                                                                                               
gcACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGTGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGcc                                                                                                                                             >  3:291782/1‑149 (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                                                                                                                                   >  1:189861/1‑148 (MQ=255)
           ggTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATc                                                                                                                                   <  2:189861/148‑1 (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)
                            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                                                                                                                               <  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)
                                                                                                                                        gtTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACACCt     >  5:240646/1‑149 (MQ=255)
                                                                                                                                         ttCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGGCAACGCCACCtt    >  4:287907/1‑149 (MQ=255)
                                                                                                                                            aTCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACACCTTGa  >  5:393142/1‑148 (MQ=255)
                                                                                                                                            aTCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACACCTTGa  >  7:127196/1‑148 (MQ=255)
                                                                                                                                             tCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGaa                                                         <  1:351562/92‑1 (MQ=255)
                                                                                                                                             tCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGaa                                                         >  2:351562/1‑92 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCACTTACTGCGGTGACCAACTGCCGCCTTACGATCTGCCCGATCACAAAACACGCATGACCATCAAGAGCCAGACCCACAAGGGGGAAGGTTTCAACGAGCTGCGTTTCGAGGATGAGCTGGGCAAGGAGGAAGTGTTCATCCATGCCCAGAGAGATCAGAACACTCTGGTCAATCACAACCAGTCCCTCTCTGTGGGCGTCGATCGCACCCAACAGGTTGGCCAGGACGAATCGGTCGCCATTGGCCGCAATCGTCTGCGCGTCGTCAAGGCCAACGACACCTTGA  >  NC_002947/4612529‑4612816

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