Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,514,335 A→T 6.7% intergenic (‑27/+18) fecB ← / ← fecA ferric citrate ABC transporter periplasmic binding protein/ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,3350AT 6.7% 276.0 / 6.3 88intergenic (‑27/+18)fecB/fecAferric citrate ABC transporter periplasmic binding protein/ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (36/46);  new base T (2/4);  total (38/50)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCA  >  NC_000913/4514189‑4514483
                                                                                                                                                  |                                                                                                                                                    
aaTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc                                                                                                                                                   <  2:323058/150‑1 (MQ=255)
         tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCtttt                                                                                                                                          <  1:1027485/150‑2 (MQ=255)
         tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTa                                                                                                                                          <  1:929265/150‑1 (MQ=255)
         tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTa                                                                                                                                          >  1:566943/1‑150 (MQ=255)
         tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTa                                                                                                                                          <  1:42015/150‑1 (MQ=255)
           gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAac                                                                                                                                        <  2:1542485/150‑1 (MQ=255)
            cGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAATGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttagc                                                                                                                                       <  2:1301531/150‑5 (MQ=255)
            cGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttt                                                                                                                                           >  1:1819238/1‑146 (MQ=255)
            cGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttt                                                                                                                                           <  2:1819238/146‑1 (MQ=255)
             gTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttagct                                                                                                                                      <  1:625618/150‑6 (MQ=255)
             gTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                      <  1:1938846/150‑1 (MQ=255)
              ttttttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAAcaca                                                                                                                                     <  1:613164/150‑1 (MQ=255)
                   tCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGa                                                                                                                                <  1:1033477/150‑1 (MQ=255)
                   tCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGa                                                                                                                                >  2:165048/1‑150 (MQ=255)
                     gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                              <  2:470398/150‑1 (MQ=255)
                      aGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACt                                                                                                                             <  1:1705785/150‑1 (MQ=255)
                         gtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCa                                                                                                                          <  2:1026734/150‑1 (MQ=255)
                         gtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCa                                                                                                                          <  1:1438252/150‑1 (MQ=255)
                         gtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCa                                                                                                                          <  1:1167669/150‑1 (MQ=255)
                           aaaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAc                                                                                                                        <  2:1258925/150‑1 (MQ=255)
                            aaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACg                                                                                                                       <  2:302280/150‑1 (MQ=255)
                              cGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAc                                                                                                                     <  1:818526/150‑1 (MQ=255)
                                    gtgtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGc                                                                                                               >  2:1165120/1‑150 (MQ=255)
                                     tgtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                              <  1:139813/150‑1 (MQ=255)
                                     tgtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                              <  2:228360/150‑1 (MQ=255)
                                      gtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCat                                                                                                             <  2:404056/150‑1 (MQ=255)
                                          gTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATAc                                                                                                         <  2:1631254/150‑1 (MQ=255)
                                           tCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACa                                                                                                        >  1:941479/1‑150 (MQ=255)
                                              tGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCg                                                                                                     <  2:73016/150‑1 (MQ=255)
                                                aaCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTg                                                                                                   <  2:781122/150‑1 (MQ=255)
                                                 aCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTgc                                                                                                  >  1:1610136/1‑150 (MQ=255)
                                                      ggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGc                                                                                             >  1:127556/1‑150 (MQ=255)
                                                         ggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTgg                                                                                          >  1:1763876/1‑150 (MQ=255)
                                                           cAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCATCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGcc                                                                                        >  2:1399463/1‑150 (MQ=255)
                                                            aaaGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGCCCCCTGCATATACAGCGTGCGCGGCTGGCCt                                                                                       >  2:1900642/1‑150 (MQ=255)
                                                               ggCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAACGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCa                                                                                    >  1:772849/1‑150 (MQ=255)
                                                                gCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCAt                                                                                   <  2:1620892/150‑1 (MQ=255)
                                                                 cGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATa                                                                                  <  2:1573608/150‑1 (MQ=255)
                                                                      cTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGc                                                                             >  1:587039/1‑150 (MQ=255)
                                                                       tGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAACGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGcc                                                                            >  1:864167/1‑150 (MQ=255)
                                                                       tGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGcc                                                                            >  2:977916/1‑150 (MQ=255)
                                                                       tGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGcc                                                                            >  1:1351105/1‑150 (MQ=255)
                                                                        gATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGcc                                                                                        <  1:739093/137‑1 (MQ=255)
                                                                        gATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGcc                                                                                        >  2:739093/1‑137 (MQ=255)
                                                                         aTCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCACCGACCCATGCATATACAGCGAGCGCGGCAGGCCTGCATAGATGCCtt                                                                          >  1:1077853/1‑150 (MQ=255)
                                                                           cACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttg                                                                        >  1:944683/1‑150 (MQ=255)
                                                                            aCCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTg                                                                                                                <  2:1882788/109‑1 (MQ=255)
                                                                            aCCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTg                                                                                                                >  1:1882788/1‑109 (MQ=255)
                                                                            aCCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgt                                                                       >  2:167255/1‑150 (MQ=255)
                                                                               agcagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgt                                                                    <  2:1351105/150‑1 (MQ=255)
                                                                                 cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                  <  2:941479/150‑1 (MQ=255)
                                                                                 cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                  <  2:1577320/150‑1 (MQ=255)
                                                                                   gcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGCCGTTAACACATCAGAAATTCAACTACCAATGCATAAACAGAGTGCGAGGCTGACCTTCAAAGATACCTTTGTTGTTgtc                                                                >  2:1720502/1‑150 (MQ=255)
                                                                                        ccTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCAt                                                           >  1:1065175/1‑150 (MQ=255)
                                                                                         cTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATa                                                          >  1:954321/1‑150 (MQ=255)
                                                                                          tGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCt                                                                           <  1:1916632/132‑1 (MQ=255)
                                                                                          tGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCt                                                                           >  2:1916632/1‑132 (MQ=255)
                                                                                            cAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                  <  1:1821981/139‑1 (MQ=255)
                                                                                            cAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                  >  2:1821981/1‑139 (MQ=255)
                                                                                            cAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAaga                                                       >  2:614906/1‑150 (MQ=255)
                                                                                              aaaaGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGc                                                     <  1:167255/150‑1 (MQ=255)
                                                                                                  gAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGAt                                                 <  2:809417/150‑1 (MQ=255)
                                                                                                   aaaaCGGATAAATGCCATCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATg                                                <  1:1399463/150‑1 (MQ=255)
                                                                                                     aaCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGaa                                              >  1:615690/1‑150 (MQ=255)
                                                                                                       cGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagt                                            <  1:1165120/150‑1 (MQ=255)
                                                                                                           tAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTc                                        >  1:1638764/1‑150 (MQ=255)
                                                                                                           tAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTc                                        >  2:450051/1‑150 (MQ=255)
                                                                                                            aaaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTcc                                       <  1:86896/150‑1 (MQ=255)
                                                                                                             aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCt                                      >  2:497977/1‑150 (MQ=255)
                                                                                                                gCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGt                                   <  1:977916/150‑1 (MQ=255)
                                                                                                                  cAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCg                                 >  2:5671/1‑150 (MQ=255)
                                                                                                                       tAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGAt                            <  2:1610136/150‑1 (MQ=255)
                                                                                                                         aTCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGt                          <  1:450051/150‑1 (MQ=255)
                                                                                                                         aTCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGt                          >  1:1325555/1‑150 (MQ=255)
                                                                                                                            acaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGtttt                       >  2:1163502/1‑150 (MQ=255)
                                                                                                                              aTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                     <  1:614906/150‑1 (MQ=255)
                                                                                                                                  cAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCacac                 <  1:64872/150‑1 (MQ=255)
                                                                                                                                      tAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGaa             >  1:1216873/1‑150 (MQ=255)
                                                                                                                                        aacgCCCGGCTTGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg           >  1:1885960/4‑150 (MQ=255)
                                                                                                                                        aaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg           <  2:429244/150‑1 (MQ=255)
                                                                                                                                        aaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg           >  2:1283208/1‑150 (MQ=255)
                                                                                                                                        aaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg           <  2:772849/150‑1 (MQ=255)
                                                                                                                                         aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGc          >  1:562711/1‑150 (MQ=255)
                                                                                                                                          aGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGcc         <  2:1763876/150‑1 (MQ=255)
                                                                                                                                            cccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAg       <  1:497977/150‑1 (MQ=255)
                                                                                                                                              cGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGt     <  1:1730155/150‑1 (MQ=255)
                                                                                                                                               ggCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGtt    >  2:1754025/1‑150 (MQ=255)
                                                                                                                                                 cAAGCCGGGCGTTAAGACATCAGAACTTCAACGACCGCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa  <  1:1720502/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCA  >  NC_000913/4514189‑4514483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: