Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,496 T→G 5.2% S232R (AGT→CGT)  fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,4960TG 5.2% 352.8 / 5.2 120S232R (AGT→CGT) fepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (63/48);  new base G (4/3);  total (69/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCA  >  NC_000913/621349‑621643
                                                                                                                                                   |                                                                                                                                                   
aCCAGCGATGCGCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGGCAATACAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGc                                                                                                                                                   <  2:389674/150‑1 (MQ=255)
  cAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGccc                                                                                                                                                 >  2:1745268/1‑150 (MQ=255)
   aGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCa                                                                                                                                                <  2:2026071/150‑1 (MQ=255)
    gCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAg                                                                                                                                               >  1:57727/1‑150 (MQ=255)
      gatgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCATTAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcg                                                                                                                                             <  1:2052479/150‑1 (MQ=255)
          atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGt                                                                                                                                         <  2:627795/150‑1 (MQ=255)
            cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCAAGCAGACCAATCAACTGTGTGCGCGCCCCGCGACTGCCCAGCGCCGTcg                                                                                                                                       >  2:1062764/1‑150 (MQ=255)
             ggCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCCCTACCCCAAAGCACGGTACTCGCCAGCAGACCAATCCCCTGTGTGCGCGCCCCGCGACTGCCCAGCGCCGTcgc                                                                                                                                      >  1:137042/1‑150 (MQ=255)
               cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                                                                                                                    <  2:20579/150‑1 (MQ=255)
                 cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCCATACAGGCAATCGTGCGAACTATTTCCGTCGCACTACCAAAAAGCACGCTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCCCGCGACGGCGCTGGGCAGTcgcgtgg                                                                                                                                  <  1:1002639/150‑4 (MQ=255)
                 cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTCTTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtgg                                                                                                                                  <  1:479475/150‑4 (MQ=255)
                 cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACAACCACAAAGCACGGTAATCGCCAGCAGCCCAATCAACTGAGTGCGCGCCCCGCGCCTGCCCAGCGCCGTCGCGgtg                                                                                                                                  >  2:1286628/1‑150 (MQ=255)
                    accaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCg                                                                                                                               <  1:186924/150‑1 (MQ=255)
                    accaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCg                                                                                                                               <  1:519429/150‑1 (MQ=255)
                    accaGCCAACGCGACATATGCGGCATCATCAGGCCAAAAAAGGCAATCGCGCCAACTATTGCCGTCGCACTACCACAACGCCCGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGAGTCg                                                                                                                               >  1:2022984/1‑150 (MQ=255)
                     ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCTCGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGc                                                                                                                              <  2:1083744/150‑1 (MQ=255)
                        gCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGc                                                                                                                           >  2:845809/1‑150 (MQ=255)
                             cgagcCATATGCGGCATCATCAGGCCAAAAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCCACAGCCCAGCGCCGTCGCGTTGACGCTGCCGAgg                                                                                                                      >  2:1835115/4‑150 (MQ=255)
                              gcgcCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACCGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGc                                                                                                                     <  2:59752/150‑1 (MQ=255)
                                       gCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGt                                                                                                            >  2:1081202/1‑150 (MQ=255)
                                         ggCATCATCAGGCCACTAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCCGACCAATCAACTGTGGGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                          <  2:406281/150‑1 (MQ=255)
                                           catcatCAGGCCAATAAATGCAATCGGGCCAACTATTGCCGTCGCACTACCACACAGCACGGTAATCGCCAGCAGACCAATCATCTGTGTGCGCGCCACGCGACGGCCCATCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAg                                                                                                        <  1:872875/150‑1 (MQ=255)
                                              catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCTACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                     <  2:421795/150‑1 (MQ=255)
                                               atcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcg                                                                                                    >  2:924319/1‑150 (MQ=255)
                                               atcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACCCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcg                                                                                                    >  1:1311934/1‑150 (MQ=255)
                                                 caGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCcgcgcgcg                                                                                                  >  1:930245/1‑150 (MQ=255)
                                                 caGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGcgcgcg                                                                                                  <  2:1057558/150‑1 (MQ=255)
                                                   ggCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTTATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAc                                                                                                >  1:1335312/1‑150 (MQ=255)
                                                           aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                        >  2:1615579/1‑150 (MQ=255)
                                                           aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                        >  2:581857/1‑150 (MQ=255)
                                                           aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                        >  1:471392/1‑150 (MQ=255)
                                                             ggCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAg                                                                                      <  2:1823004/150‑1 (MQ=255)
                                                             ggCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAg                                                                                      >  2:654536/1‑150 (MQ=255)
                                                              gCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCTACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGGCGCTGCCGAGGCTCAAACTGTTCCGCGCGCGGCTCAGTAATAgc                                                                                     >  1:99875/1‑150 (MQ=255)
                                                               cAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTTAGTAATAgcg                                                                                    >  2:451546/1‑150 (MQ=255)
                                                               cAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGCTGTCGCTGCCGAGGCTCCCACTGTTCAGCGCGCGACTCAGTAATAgcg                                                                                    >  1:1361784/1‑150 (MQ=255)
                                                                    gggCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGCCTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAGCTGGTCAGCGCGCGCGTCAGTAATAGCGCAGtt                                                                               >  1:1266230/1‑150 (MQ=255)
                                                                     ggCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTg                                                                              <  1:1582728/150‑1 (MQ=255)
                                                                        cAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTc                                                                           >  2:1398529/1‑150 (MQ=255)
                                                                            tATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCGGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGc                                                                       <  1:99147/150‑1 (MQ=255)
                                                                             aTTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCg                                                                      >  1:1023119/1‑150 (MQ=255)
                                                                             aTTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCg                                                                      <  1:1303058/150‑1 (MQ=255)
                                                                                gCCGTCGCACTACCCCACACCACGGTAATCGCCCGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGCGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                   <  2:97099/150‑1 (MQ=255)
                                                                                gCCGTCGCACTACCACAAAGAACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGACACGCGACTGCCCAGCGCCGTCGCGGTGTCGCGGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                   <  2:1303221/150‑1 (MQ=255)
                                                                                gCCGTCGCACTAACACAAAGCACGGTAATCGCCAGAAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                   <  2:887869/150‑1 (MQ=255)
                                                                                 ccGTCGCACTACCACACAGCACGGTACTCGCCAGCCGACCAATCAACAGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGGGTCGCTGCCGAGGCGCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTAGCTCCGGCGATCa                                                                  <  1:147429/150‑1 (MQ=255)
                                                                                   gTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTTTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGc                                                                >  1:1325776/1‑150 (MQ=255)
                                                                                         cTACCACACAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGTCTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCggg                                                          <  1:1286629/150‑1 (MQ=255)
                                                                                          tACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGa                                                         >  2:1458257/1‑150 (MQ=255)
                                                                                             cacaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCa                                                      <  2:939572/150‑1 (MQ=255)
                                                                                             cacaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCa                                                      <  2:2004556/150‑1 (MQ=255)
                                                                                              acaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAg                                                     >  1:906635/1‑150 (MQ=255)
                                                                                               caAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGc                                                    >  2:1590946/1‑150 (MQ=255)
                                                                                                  aGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGGCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGcacc                                                 <  1:1835116/150‑1 (MQ=255)
                                                                                                  aGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGcacc                                                 >  1:72013/1‑150 (MQ=255)
                                                                                                      cGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCTCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCtt                                             >  2:1524665/1‑150 (MQ=255)
                                                                                                        gTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACtttt                                           >  1:1632966/1‑150 (MQ=255)
                                                                                                         tAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTa                                          >  1:1442761/1‑150 (MQ=255)
                                                                                                          aaTCGCCAGCAGACCAATCAACTGTGTGTGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCCCCGGGATCAGCACCACTTTTat                                         >  2:1456636/1‑149 (MQ=255)
                                                                                                          aaTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTaa                                         >  2:1454095/1‑150 (MQ=255)
                                                                                                            tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTTCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                       <  1:28905/150‑1 (MQ=255)
                                                                                                            tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                       <  1:636476/150‑1 (MQ=255)
                                                                                                            tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                       <  1:764848/150‑1 (MQ=255)
                                                                                                            tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                       >  2:1572886/1‑150 (MQ=255)
                                                                                                            acGCCAGCAGACCAATCACCTGTGTCCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCGCAAACTGATCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                       <  1:765629/149‑1 (MQ=255)
                                                                                                             cGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGt                                      >  2:929878/1‑150 (MQ=255)
                                                                                                              gCCAGCTGACCAATCTACTGTGTGGGCGACAGGCTACTGCCCTGCGCCGTCGCGGTGTCGCTGCCGAGGCGCACACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTa                                     <  2:339130/150‑1 (MQ=255)
                                                                                                              gCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTTTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTa                                     <  2:1325775/150‑1 (MQ=255)
                                                                                                              gCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTa                                     <  1:1764049/150‑1 (MQ=255)
                                                                                                                 agcagACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCTGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGt                                  <  2:1219905/150‑1 (MQ=255)
                                                                                                                    agACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGa                               >  2:899492/1‑150 (MQ=255)
                                                                                                                    agACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGa                               >  2:891831/1‑150 (MQ=255)
                                                                                                                     gACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCAGTTTTAAGGTATGTAGAt                              >  2:1253407/1‑150 (MQ=255)
                                                                                                                       ccAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGCTTg                            >  2:244681/1‑150 (MQ=255)
                                                                                                                        cAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGc                           >  2:1749694/1‑150 (MQ=255)
                                                                                                                         aaTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCCCCACTTTTAAGGTATGTAGATTGCg                          >  1:1961426/1‑150 (MQ=255)
                                                                                                                         aaTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCg                          <  2:622583/150‑1 (MQ=255)
                                                                                                                          aTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCCCCACATTTAAGGTATGTAGATTGCGa                         >  1:1032827/1‑150 (MQ=255)
                                                                                                                           tCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGaa                        <  1:1981546/150‑1 (MQ=255)
                                                                                                                                tgtgtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATcc                   >  2:1585301/1‑150 (MQ=255)
                                                                                                                                   gtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGc                <  2:1311933/150‑1 (MQ=255)
                                                                                                                                   gtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGACGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCTTTTAAGGTATGTAGCTTGCGAATATCCAGc                >  2:1453030/1‑150 (MQ=255)
                                                                                                                                    ggcgcgcCACGCGACTGCCCAGAGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCg               <  1:178062/149‑1 (MQ=255)
                                                                                                                                     gcgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCGGAGGCTCAACCTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGCTCAGCACCACGTTTAAGGTATGTAGATTGCGAATATCCAGCGa              >  1:76628/1‑150 (MQ=255)
                                                                                                                                     gcgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGCCTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTGTTAAGGTATGTAGATTGCGAATATCCAGCGa              >  1:789157/1‑150 (MQ=255)
                                                                                                                                      cgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGCCTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGac             >  1:1998600/1‑149 (MQ=255)
                                                                                                                                      cgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCCCCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGaa             >  1:761850/1‑150 (MQ=255)
                                                                                                                                      cgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGaa             >  2:2002422/1‑150 (MQ=255)
                                                                                                                                        gacaccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAAcc           >  1:1106386/6‑150 (MQ=255)
                                                                                                                                        cgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTCAGAGCGCAGTTGCTCCGGCGATCCGCGCCGGGATCAGCCCCCCGTTTAAGGTATGTAGATGGCGAATATCCCGCGCAcc           >  1:330511/1‑150 (MQ=255)
                                                                                                                                        cgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAAcc           >  1:1423433/1‑150 (MQ=255)
                                                                                                                                          ccACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCTTTTAAGGTATGTAGCGTGCGAAGATCCCGCGCACCgg         >  1:345753/1‑150 (MQ=255)
                                                                                                                                          caccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCgg         >  1:1113652/4‑150 (MQ=255)
                                                                                                                                           cACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGc        <  1:1037826/150‑1 (MQ=255)
                                                                                                                                            ccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCt       >  1:448823/2‑150 (MQ=255)
                                                                                                                                            ccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCt       >  2:929885/2‑150 (MQ=255)
                                                                                                                                               cgACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAAGTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGc    <  2:743296/150‑1 (MQ=255)
                                                                                                                                               cgACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGc    >  1:606313/1‑150 (MQ=255)
                                                                                                                                                 aCTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTGCAGCGCGCGACTCAGTAATAGCGCAGTTGCGCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCa  <  1:1029091/150‑1 (MQ=255)
                                                                                                                                                 aCTGCCCAGCGCCGTCGCGGTGTCGCGGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCa  <  2:1868912/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCA  >  NC_000913/621349‑621643

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