Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 5.3% coding (65/957 nt) fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,3501.G 5.3% 553.5 / 26.9 133coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (66/60);  new base G (5/2);  total (71/62)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 141 total aligned reads displayed.

TCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCA  >  NC_000913/4512201‑4512499
                                                                                                                                                      |                                                                                                                                                     
tCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                        >  2:1433888/1‑150 (MQ=255)
 cAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                      <  1:1040682/150‑1 (MQ=255)
 cAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAAGACATAATAATGCTCGTGTCCGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                      >  1:1942442/1‑150 (MQ=255)
 cAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGACAGCAGCAAGCGCGGCAGTCGGAACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCTGACAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                      >  1:1185623/1‑150 (MQ=255)
  aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑gg                                                                                                                                                     >  2:1877436/1‑150 (MQ=255)
   gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGa                                                                                                                                                    >  2:599812/1‑150 (MQ=255)
     aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGt                                                                                                                                                  <  2:798558/150‑1 (MQ=255)
      cgcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTg                                                                                                                                                 <  1:1206429/150‑1 (MQ=255)
        cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGAt                                                                                                                                               >  1:1508557/1‑150 (MQ=255)
        cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGAt                                                                                                                                               >  2:1468101/1‑150 (MQ=255)
         ccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATa                                                                                                                                              <  1:1877436/150‑1 (MQ=255)
           gccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAAAGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACT‑CCATATGGGAGTGAtccc                                                                                                                                            <  2:1948259/150‑4 (MQ=255)
            ccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACa                                                                                                                                           >  2:5404/1‑150 (MQ=255)
            ccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACT‑CCATATGGGAGTGAtcccc                                                                                                                                           <  1:688892/150‑5 (MQ=255)
               cGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGGCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgcg                                                                                                                                        <  1:1433812/150‑1 (MQ=255)
                tgCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGAGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTTCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgcgc                                                                                                                                       >  1:2208763/2‑150 (MQ=255)
                ggCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgcgc                                                                                                                                       >  2:830019/1‑150 (MQ=255)
                  cGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGcac                                                                                                                                     <  1:821186/150‑1 (MQ=255)
                  cGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGcac                                                                                                                                     <  2:440904/150‑1 (MQ=255)
                     gggCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGc                                                                                                                                  >  2:1124908/1‑150 (MQ=255)
                       gCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGccc                                                                                                                                <  2:303448/150‑1 (MQ=255)
                       gCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGccc                                                                                                                                >  2:1293603/1‑150 (MQ=255)
                         tGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCGGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGc                                                                                                                              >  1:1951459/1‑150 (MQ=255)
                           cACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCa                                                                                                                            >  1:1542110/1‑150 (MQ=255)
                                aCAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCa                                                                                                                       <  2:1047908/150‑1 (MQ=255)
                                         gCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                              >  2:1582119/1‑150 (MQ=255)
                                         gCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCGGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                              <  1:645860/150‑1 (MQ=255)
                                           cagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                            <  1:169415/150‑1 (MQ=255)
                                             gcagcaAGCGCGGCAGTCGTTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                          <  2:1949137/150‑1 (MQ=255)
                                             gcagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                          <  1:2099353/150‑1 (MQ=255)
                                               agcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                        >  2:782076/1‑150 (MQ=255)
                                                 caAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                      >  2:537278/1‑150 (MQ=255)
                                                   aGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCa                                                                                                    >  2:244399/1‑150 (MQ=255)
                                                   aGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCa                                                                                                    >  2:721706/1‑150 (MQ=255)
                                                    gcgcGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAg                                                                                                   <  2:712188/150‑1 (MQ=255)
                                                        ggCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCa                                                                                               <  2:439193/150‑1 (MQ=255)
                                                         gCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCTTAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                              >  2:1696400/1‑150 (MQ=255)
                                                            gTCGGGACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑TGAGTGATTACTGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAttt                                                                                           <  1:1428715/150‑1 (MQ=255)
                                                             tCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACCGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtttt                                                                                          >  1:2153872/1‑150 (MQ=255)
                                                             tCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtttt                                                                                          <  2:1241255/150‑1 (MQ=255)
                                                             tCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtttt                                                                                          >  2:728249/1‑150 (MQ=255)
                                                             tCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAAAACCAGCGCAAtatt                                                                                          >  2:1440260/1‑150 (MQ=255)
                                                               ggTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCTCCAAGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                        <  1:305457/150‑1 (MQ=255)
                                                               ggTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                        >  1:2058550/1‑150 (MQ=255)
                                                               ggTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                        <  1:1260298/150‑1 (MQ=255)
                                                               ggTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                        <  1:1338917/150‑1 (MQ=255)
                                                                gTACTCCATCAATACATAATAATGCTCGTGTCCGGCCCGCCAGTCGGTCAGCAGCGCGCACCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCTAGGGTGAAGAAATTAAACAGAGCAATTTTCAt                                                                                       >  1:370553/1‑150 (MQ=255)
                                                                  aCTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcc                                                                                     <  2:1698845/150‑1 (MQ=255)
                                                                   cTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCCGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                                                    >  1:356614/1‑150 (MQ=255)
                                                                     ccATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                                  <  2:2019636/150‑1 (MQ=255)
                                                                        tCAATACATAATAATGCTCGAGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcct                                                                               <  1:1270729/150‑1 (MQ=255)
                                                                         cAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                              >  2:1022252/1‑150 (MQ=255)
                                                                                 aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCTCCAGGGCAAGGGTGCTGAAAATAACCAGCGCAAATTTCATCCTCGCCTCCTCACaac                                                                      >  2:32420/1‑149 (MQ=255)
                                                                                 aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAg                                                                      >  2:885945/1‑150 (MQ=255)
                                                                                  ataatGCTCGTGTCCGGCCTGCCAGTCGGTAATCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGATTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGc                                                                     <  2:1185623/150‑1 (MQ=255)
                                                                                    aatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                                   <  1:1028461/150‑1 (MQ=255)
                                                                                    aatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCACACAAGCCa                                                                   >  1:576321/1‑150 (MQ=255)
                                                                                     atGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAg                                                                  <  2:82283/150‑1 (MQ=255)
                                                                                     atGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAg                                                                  <  2:626089/150‑1 (MQ=255)
                                                                                      tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                 <  1:1582119/150‑1 (MQ=255)
                                                                                       gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAc                                                                <  1:1296463/150‑1 (MQ=255)
                                                                                        cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                               >  1:1994533/1‑150 (MQ=255)
                                                                                        cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                               >  2:88137/1‑150 (MQ=255)
                                                                                               ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCCCCAGCCAGACAAAGCAAg                                                        >  1:1131979/1‑150 (MQ=255)
                                                                                                cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                       >  2:2057378/1‑150 (MQ=255)
                                                                                                cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                       <  1:632454/150‑1 (MQ=255)
                                                                                                  gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGc                                                     >  1:1426527/1‑150 (MQ=255)
                                                                                                      gCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                                 <  1:1597165/150‑1 (MQ=255)
                                                                                                      gCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                                 >  1:247349/1‑150 (MQ=255)
                                                                                                      gCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                                 <  1:1993465/150‑1 (MQ=255)
                                                                                                       ccAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCa                                                <  1:1556605/150‑1 (MQ=255)
                                                                                                          gTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATc                                             <  2:490197/150‑1 (MQ=255)
                                                                                                            cGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCATCGCAATTTTCTTCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCac                                           >  2:943052/1‑149 (MQ=255)
                                                                                                            cGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAg                                           >  2:1033257/1‑150 (MQ=255)
                                                                                                              gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgcg                                         <  2:338503/150‑1 (MQ=255)
                                                                                                               gcagcagCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgcgc                                        <  2:1542110/149‑1 (MQ=255)
                                                                                                                  gcagcGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAg                                     >  1:1184538/1‑150 (MQ=255)
                                                                                                                  gcagcGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAg                                     <  1:480246/150‑1 (MQ=255)
                                                                                                                  gcagcGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAg                                     >  2:106942/1‑150 (MQ=255)
                                                                                                                    agcGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAGAACCAGAGAAATTATCATCCTCGCCGCCTCACAAGCCAGCCAAAGCAGGGGCTGCCAATCAGCGACAgaa                                   >  1:331169/1‑148 (MQ=255)
                                                                                                                       gcgcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgc                                          <  1:1251032/140‑1 (MQ=255)
                                                                                                                       gcgcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgc                                          >  2:1251032/1‑140 (MQ=255)
                                                                                                                          cgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcgc                             >  1:611645/1‑150 (MQ=255)
                                                                                                                            ccACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCt                           >  1:1004319/1‑150 (MQ=255)
                                                                                                                            ccACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCt                           >  2:2100012/1‑150 (MQ=255)
                                                                                                                             cacggcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                          >  2:2119024/1‑150 (MQ=255)
                                                                                                                               cggcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                        <  1:88137/150‑1 (MQ=255)
                                                                                                                                 gcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCGAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAACCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCggg                      <  2:1899978/150‑1 (MQ=255)
                                                                                                                                  cacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGc                     >  2:1784884/1‑150 (MQ=255)
                                                                                                                                   gctgttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCa                    >  2:1067781/7‑150 (MQ=255)
                                                                                                                                   gctgttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCa                    >  2:1499156/7‑150 (MQ=255)
                                                                                                                                    cggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg                   <  2:2001109/150‑1 (MQ=255)
                                                                                                                                     ggggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGa                  <  1:388846/150‑1 (MQ=255)
                                                                                                                                     ggggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGa                  <  1:1184375/150‑1 (MQ=255)
                                                                                                                                      gggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt                 >  2:1194172/1‑150 (MQ=255)
                                                                                                                                            aCTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgg           >  1:1751723/1‑150 (MQ=255)
                                                                                                                                             cTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg          >  2:1335358/1‑150 (MQ=255)
                                                                                                                                             cTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg          <  1:137599/150‑1 (MQ=255)
                                                                                                                                             cTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg          <  1:1114163/150‑1 (MQ=255)
                                                                                                                                                     t‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACTAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGccc  >  1:435654/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCA  >  NC_000913/4512201‑4512499

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