Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 610,207 C→G 7.6% intergenic (‑128/+47) entD ← / ← fepA phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913610,2070CG 7.6% 317.3 / 15.0 94intergenic (‑128/+47)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/ferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (46/38);  new base G (4/3);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 102 total aligned reads displayed.

GGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCC  >  NC_000913/610062‑610353
                                                                                                                                                 |                                                                                                                                                  
ggTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                                 <  3:191082/149‑1 (MQ=255)
   aTGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                              <  2:1722784/149‑1 (MQ=21)
                        cATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAac                                                                                                                                       >  3:1852906/1‑135 (MQ=255)
                        cATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAac                                                                                                                                       <  4:1852906/135‑1 (MQ=255)
                                     gatgaCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACg                                                                                                            >  4:1461406/1‑149 (MQ=255)
                                      atgaCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGt                                                                                                           >  2:1407184/1‑149 (MQ=255)
                                           aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCa                                                                                                      >  4:1854123/1‑149 (MQ=255)
                                                     cTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                                                            <  1:1407184/149‑1 (MQ=255)
                                                     cTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                                                            <  3:1854123/149‑1 (MQ=255)
                                                           gATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                      >  1:461921/1‑149 (MQ=255)
                                                            aTGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                                                     <  3:1461406/149‑1 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTc                                                                                                                        <  4:1502033/108‑1 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTc                                                                                                                        >  3:1502033/1‑108 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATc                                                                                                                 <  3:1340991/115‑1 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATc                                                                                                                 >  4:1340991/1‑115 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                               >  4:1576345/1‑149 (MQ=255)
                                                                    cGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata                                                                             >  2:1282370/1‑149 (MQ=255)
                                                                     gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                            >  2:1330500/1‑149 (MQ=255)
                                                                     gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                            >  3:998570/1‑149 (MQ=255)
                                                                     gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                            >  1:1639280/1‑149 (MQ=255)
                                                                      cTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCAtata                                                                           >  1:1238698/1‑149 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata                                                                             <  3:840780/141‑1 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata                                                                             >  4:840780/1‑141 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                                     >  4:994083/1‑149 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCa                                                                        >  1:1680436/1‑146 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCa                                                                        <  2:1680436/146‑1 (MQ=255)
                                                                             ggCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGa                                                                                                              <  3:414937/107‑1 (MQ=255)
                                                                             ggCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGa                                                                                                              >  4:414937/1‑107 (MQ=255)
                                                                                cTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgtacgt                                                                 >  3:386893/1‑149 (MQ=255)
                                                                                  aCATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                      <  2:1218034/126‑1 (MQ=255)
                                                                                  aCATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                      >  1:1218034/1‑126 (MQ=255)
                                                                                     tGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                                            >  2:596410/1‑149 (MQ=255)
                                                                                     tGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                                            <  2:1639280/149‑1 (MQ=255)
                                                                                        tCCTGCAATATTTTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGACGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                         >  3:1244398/1‑149 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             >  1:424251/1‑62 (MQ=17)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             <  2:424251/62‑1 (MQ=17)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          >  1:672695/1‑113 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          <  2:672695/113‑1 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                         <  4:188261/146‑1 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                         >  3:188261/1‑146 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                      >  3:660038/1‑149 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                      >  1:886231/1‑149 (MQ=255)
                                                                                            gCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                                                     <  4:1779026/149‑1 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                   <  1:1282370/149‑1 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                   <  3:994083/149‑1 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                   >  4:410779/1‑149 (MQ=255)
                                                                                               atatatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                               <  4:693917/120‑1 (MQ=255)
                                                                                               atatatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                               >  3:693917/1‑120 (MQ=255)
                                                                                                  tatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          >  2:1176024/1‑106 (MQ=255)
                                                                                                  tatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          <  1:1176024/106‑1 (MQ=255)
                                                                                                    ttGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGtata                                             >  3:1127568/1‑149 (MQ=255)
                                                                                                     tGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATAc                                            <  3:1576345/149‑1 (MQ=255)
                                                                                                       aaTTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGc                                          >  4:795736/1‑149 (MQ=255)
                                                                                                        aTTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcc                                         >  3:1496348/1‑149 (MQ=255)
                                                                                                         ttGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGccc                                        <  4:1127568/149‑1 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                           <  4:1205169/47‑1 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                           >  3:1205169/1‑47 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                      >  2:677980/1‑100 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                      <  1:677980/100‑1 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                   >  4:302766/1‑149 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                   >  4:171413/1‑149 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                   >  3:1694646/1‑149 (MQ=255)
                                                                                                                  tttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCgg                               >  1:1740174/1‑149 (MQ=255)
                                                                                                                  tttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCgg                               >  4:1818262/1‑149 (MQ=255)
                                                                                                                     ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGa                            <  2:1238698/149‑1 (MQ=255)
                                                                                                                      tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGat                           <  3:171413/149‑1 (MQ=255)
                                                                                                                       gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGata                          <  2:461921/149‑1 (MQ=255)
                                                                                                                        tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGatat                         <  3:410779/149‑1 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                         <  4:1510540/97‑1 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                         >  3:1510540/1‑97 (MQ=255)
                                                                                                                             cGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCAtata                                                                           <  1:1120032/94‑1 (MQ=255)
                                                                                                                             cGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCAtata                                                                           >  2:1120032/1‑94 (MQ=255)
                                                                                                                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTgg                   <  2:886231/149‑1 (MQ=255)
                                                                                                                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGACACCGGCGATATAGTTgg                   >  1:1011875/1‑149 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt              <  4:386893/149‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                         >  1:1106610/1‑89 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                         <  2:1106610/89‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTg             <  1:596410/149‑1 (MQ=255)
                                                                                                                                     gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGcaggggc            >  4:434713/1‑143 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTa                                                                                                         >  4:1107749/1‑53 (MQ=18)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTa                                                                                                         <  3:1107749/53‑1 (MQ=21)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAg                       <  3:1875152/135‑1 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAg                       >  4:1875152/1‑135 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt              >  4:1539939/1‑144 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt              <  3:1539939/144‑1 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa         >  4:48200/1‑149 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          >  4:223013/1‑68 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                          <  3:223013/68‑1 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa         >  2:1439802/1‑149 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt       <  3:302766/149‑1 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt       <  4:1809317/149‑1 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt       <  4:1694646/149‑1 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                    <  4:1056245/134‑1 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                    >  3:1056245/1‑134 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                       >  4:77718/1‑96 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                       <  3:77718/96‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc  <  1:1439802/149‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc  >  3:421021/1‑149 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc  >  4:1662879/1‑149 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATACGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc  >  4:628333/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCC  >  NC_000913/610062‑610353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: