Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,227,322 Δ1 bp 5.6% coding (1785/2391 nt) PP_4606 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,227,3220A. 5.6% 209.8 / 7.0 71coding (1785/2391 nt)PP_4606ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base . (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

TTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCA  >  NC_002947/5227179‑5227469
                                                                                                                                               |                                                                                                                                                    
ccgatctAGACCGTCTCGCTTGTGACGGCGTCAAAGGTCGAGACCAGTACAGTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCTCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGcc                                                                                                                                                <  1:1046061/144‑1 (MQ=255)
  gCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCAc                                                                                                                                              <  1:438197/150‑1 (MQ=255)
      tAGACCGTCTCGCTTGTGGAGGCGACAAAGTTCGAGACCATTACACTGTGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGCCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgag                                                                                                                                         <  2:211577/150‑1 (MQ=255)
      tAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgag                                                                                                                                         <  1:1213100/150‑1 (MQ=255)
        gACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgagga                                                                                                                                       <  2:233181/150‑1 (MQ=255)
        gACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCACCCAGGTGCCACGCCTTGCCCTCCTCAGGTGCCaagatcgg                                                                                                                                        >  2:1046061/1‑143 (MQ=255)
             tctcGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTCGATCTTCCGAACTATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGa                                                                                                                                  <  2:1107345/150‑1 (MQ=255)
              ctcGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgaggagggcaa                                                                                                                                 <  2:325411/150‑6 (MQ=255)
                   ttGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGc                                                                                                                                >  1:528371/1‑146 (MQ=255)
                   ttGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGc                                                                                                                                <  2:528371/146‑1 (MQ=255)
                    tgtgACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCATTGCCGTCCTCAGGTGCCACCGAGGAGTCGAGCACGAc                                                                                                                           <  1:336445/150‑1 (MQ=255)
                          ggCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCGGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCGCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGcc                                                                                                                     <  1:608644/150‑1 (MQ=255)
                              tCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTc                                                                                                                 >  2:334221/1‑150 (MQ=255)
                                aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACTCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGa                                                                                                                            <  1:670551/137‑1 (MQ=255)
                                aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGa                                                                                                                            >  2:670551/1‑137 (MQ=255)
                                aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAg                                                                                                               >  1:159497/1‑150 (MQ=255)
                                aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCAcctgatga                                                                                                                                       <  2:355021/126‑7 (MQ=255)
                                aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgagga                                                                                                                                       >  1:355021/1‑126 (MQ=255)
                                 aaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAgc                                                                                                              <  1:334221/150‑1 (MQ=255)
                                   gTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAgcgg                                                                                                            >  2:951542/1‑150 (MQ=255)
                                    ttCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAgcggc                                                                                                           >  2:449492/1‑150 (MQ=255)
                                       gagaCCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGGGATGCGTAATGGCATGTCGCGGTTGGACCAAAAAAGATACCTACAACCCGTTTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGTGCAGCGGcgcg                                                                                                        >  1:1033313/1‑150 (MQ=255)
                                            cAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGt                                                                                                   >  1:731463/1‑150 (MQ=255)
                                             aGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTc                                                                                                  >  1:686511/1‑150 (MQ=255)
                                                     gcgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTggg                                                                                          <  2:522002/150‑1 (MQ=255)
                                                           ggTGATATCGAACAATGCCAGGGTGCCTGTGCTGAAGCATAACACGGCATACGTGGCCCGCAGCTCCCACACACCCACCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGaa                                                                                    <  1:1385769/150‑1 (MQ=255)
                                                            gTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGt                                                                                                   <  1:882678/134‑1 (MQ=255)
                                                            gTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGt                                                                                                   >  2:882678/1‑134 (MQ=255)
                                                             tGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACg                                                                                  >  1:1322258/1‑150 (MQ=255)
                                                              gATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGa                                                                                 <  1:951542/150‑1 (MQ=255)
                                                              gATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGa                                                                                 <  2:27011/150‑1 (MQ=255)
                                                                 atCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCga                                                                                                                                          >  1:1348383/1‑90 (MQ=255)
                                                                 atCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCga                                                                                                                                          <  2:1348383/90‑1 (MQ=255)
                                                                  tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                             >  2:1324633/1‑150 (MQ=255)
                                                                  tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                             >  2:72244/1‑150 (MQ=255)
                                                                  tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGACTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                             >  2:50430/1‑150 (MQ=255)
                                                                  tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGATTGGCGCCCAGCTCCAACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGCCGAGTGGGGCTTGGACGACTCg                                                                             >  2:631990/1‑150 (MQ=255)
                                                                  tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTAGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                             >  1:1363706/1‑150 (MQ=255)
                                                                    gAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCt                                                                           <  2:744697/150‑1 (MQ=255)
                                                                              aGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCa                                                                 >  1:1348946/1‑150 (MQ=255)
                                                                                 gTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAg                                                              >  2:964921/1‑150 (MQ=255)
                                                                                     cTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGa                                                                                                >  2:631498/1‑112 (MQ=255)
                                                                                     cTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGa                                                                                                <  1:631498/112‑1 (MQ=255)
                                                                                      tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGAt                                                                                                      >  2:625899/1‑105 (MQ=255)
                                                                                      tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGAt                                                                                                      <  1:625899/105‑1 (MQ=255)
                                                                                             cGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGACCGAGCGCTgg                                                  <  2:1348946/150‑1 (MQ=255)
                                                                                              ggCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGa                                                                                                >  1:344264/1‑103 (MQ=255)
                                                                                              ggCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGa                                                                                                <  2:344264/103‑1 (MQ=255)
                                                                                               gCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTgggc                                                <  1:50430/150‑3 (MQ=255)
                                                                                                        gTCTAGCTTGGCGCTCAGTTCCCACGCCTTGCCCTCTTCAGGTGCCACCGAGGAGTCGAGTACGACATTGCCGGTCAGCGGTTCAATGGTCGATTTGTTCTTGAACGACTCGCTTTATCTGCCTGAGAACGTTATCTGGTCTTCTACtgt                                       <  2:11683/150‑3 (MQ=255)
                                                                                                         tCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                      >  1:22106/1‑150 (MQ=255)
                                                                                                         tCGAGCTTGGCGCCCAGCTACCACGCCGTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                      >  2:837132/1‑150 (MQ=255)
                                                                                                             gCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTaaa                                  <  1:231346/150‑1 (MQ=255)
                                                                                                              cTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAAc                                 <  2:731463/150‑1 (MQ=255)
                                                                                                               ttGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACg                                >  1:458454/1‑150 (MQ=255)
                                                                                                                          cTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTg                     >  2:780179/1‑150 (MQ=255)
                                                                                                                           ccccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGCGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGa                    <  2:209451/149‑1 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGAc                                          >  1:766127/1‑126 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGAc                                          <  2:766127/126‑1 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                      >  1:1319813/1‑130 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                      <  2:1319813/130‑1 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAgg                  >  2:77449/1‑150 (MQ=255)
                                                                                                                             ccACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAgg                  >  2:1109805/1‑150 (MQ=255)
                                                                                                                                  ccTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTAccc             >  2:1236725/1‑150 (MQ=255)
                                                                                                                                   cTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGCCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTg                                                   <  2:96591/111‑1 (MQ=255)
                                                                                                                                   cTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGCCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTg                                                   >  1:96591/1‑111 (MQ=255)
                                                                                                                                      gCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGc         <  1:1109805/150‑1 (MQ=255)
                                                                                                                                       gactcctcGGTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGcc        >  2:1273259/3‑150 (MQ=255)
                                                                                                                                           cctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGc    <  1:1324633/150‑1 (MQ=255)
                                                                                                                                             tcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCa  >  1:955204/1‑150 (MQ=255)
                                                                                                                                             tcAGGTGCCACCGAGGAGTAGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTggacc  >  1:622480/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCA  >  NC_002947/5227179‑5227469

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: