Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,506 C→A 5.2% T228T (ACG→ACT)  fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,5060CA 5.2% 358.0 / 6.3 124T228T (ACG→ACT) fepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (68/49);  new base A (5/2);  total (73/51)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 125 total aligned reads displayed.

TGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCA  >  NC_000913/621357‑621651
                                                                                                                                                     |                                                                                                                                                 
tgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGcc                                                                                                                                                   >  2:1939664/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)
         cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGCCTGCCCAGCGCCGTCGCGgtg                                                                                                                                          >  1:1303709/1‑150 (MQ=255)
         cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACAACCACAAAGCACGGTAATCGCCAGCAGCCCAATCAACTGAGTGCGCGCCCCGCGCCTGCCCAGCGCCGTCGCGgtg                                                                                                                                          >  2:1286628/1‑150 (MQ=255)
            accaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCg                                                                                                                                       <  1:519429/150‑1 (MQ=255)
            accaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCg                                                                                                                                       <  1:186924/150‑1 (MQ=255)
             ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCTCGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGc                                                                                                                                      <  2:1083744/150‑1 (MQ=255)
             ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTAGc                                                                                                                                      >  1:870076/1‑150 (MQ=255)
               aGCCAACGCGCCATATGCGGCATCATCAGGACAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCAAGGTAATCGCCAGCAGACCAATCAACAGTGTGCGCGCCACGCGACAGCAAAGCCCCGTCGCGGTGTCGCTg                                                                                                                                    >  1:1401238/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)
                                       atcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGCTCTCGCTGCCGAGGCGCAAACTTTTCAgcgcg                                                                                                            >  2:907740/1‑150 (MQ=255)
                                       atcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACCCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcg                                                                                                            >  1:1311934/1‑150 (MQ=255)
                                         caGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGcgcgcg                                                                                                          <  2:1057558/150‑1 (MQ=255)
                                           ggCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTTATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAc                                                                                                        >  1:1335312/1‑150 (MQ=255)
                                           ggCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAc                                                                                                        <  2:544854/150‑1 (MQ=255)
                                            gCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGCCt                                                                                                       >  2:1380271/1‑150 (MQ=255)
                                               aaTACAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAg                                                                                                    <  1:310087/150‑1 (MQ=255)
                                                   aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                                >  2:581857/1‑150 (MQ=255)
                                                   aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                                >  2:1615579/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)
                                                             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)
                                                                                 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)
                                                                                                    tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTTCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                               <  1:28905/150‑1 (MQ=255)
                                                                                                    tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                               <  1:636476/150‑1 (MQ=255)
                                                                                                    tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                               <  1:764848/150‑1 (MQ=255)
                                                                                                    acGCCAGCAGACCAATCACCTGTGTCCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCGCAAACTGATCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                               <  1:765629/149‑1 (MQ=255)
                                                                                                     cGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGt                                              >  2:929878/1‑150 (MQ=255)
                                                                                                      gCCAGCTGACCAATCTACTGTGTGGGCGACAGGCTACTGCCCTGCGCCGTCGCGGTGTCGCTGCCGAGGCGCACACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTa                                             <  2:339130/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)
                                                                                                               ccAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCCCCGGGATCAGCACCACTTTTAAGGTATGTAGCTTg                                    >  2:1503757/1‑150 (MQ=255)
                                                                                                                cAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGTACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGc                                   <  2:654589/150‑1 (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)
                                                                                                                      aCTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCTTTTAAGGTATGTAGCTTGCGAatat                             >  1:1976147/1‑150 (MQ=255)
                                                                                                                        tgtgtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATcc                           <  1:1081203/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)
                                                                                                                              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)
                                                                                                                                  ccACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCTTTTAAGGTATGTAGCGTGCGAAGATCCCGCGCACCgg                 >  1:345753/1‑150 (MQ=255)
                                                                                                                                  caccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCgg                 >  1:1113652/4‑150 (MQ=255)
                                                                                                                                   cACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGc                <  2:1640668/150‑1 (MQ=255)
                                                                                                                                    ccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCt               >  2:929885/2‑150 (MQ=255)
                                                                                                                                     cgcgACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCATCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCAAGCGAACCGGCtt              >  2:1019446/1‑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)
                                                                                                                                            gCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaa       >  2:736065/1‑150 (MQ=255)
                                                                                                                                            gCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaa       >  2:1135501/1‑150 (MQ=255)
                                                                                                                                            gCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaa       >  2:1326524/1‑150 (MQ=255)
                                                                                                                                              ccAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCCGCGCGCGACTCAGTAATAGCGCCGTTGCTCCGGCGGGCAGCACCGGGATCAGCCCCTCTTTTTAGGGAGGTAGAGTGCGAATCTCCAGCGCACCGGCTTGCCAgccac     >  1:1250963/1‑146 (MQ=255)
                                                                                                                                              ccAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTTAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAAc     <  1:451546/150‑1 (MQ=255)
                                                                                                                                                aGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGc   <  1:845810/150‑1 (MQ=255)
                                                                                                                                                 gcgcCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGGAGCTTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCa  >  1:480904/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCA  >  NC_000913/621357‑621651

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