Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 610,210 C→G 6.7% intergenic (‑131/+44) 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,2100CG 6.7% 375.4 / 14.9 105intergenic (‑131/+44)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/ferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (53/45);  new base G (4/3);  total (57/48)
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.44e-01

Only 100 of 106 total aligned reads displayed.

GGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGG  >  NC_000913/610062‑610357
                                                                                                                                                    |                                                                                                                                                   
ggTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                                     <  3:191082/149‑1 (MQ=255)
   aTGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                  <  2:1722784/149‑1 (MQ=21)
                        cATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAac                                                                                                                                           >  3:1852906/1‑135 (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                                                                                                                     >  4:1340991/1‑115 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATc                                                                                                                     <  3:1340991/115‑1 (MQ=255)
                                                                  tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                                   >  4:1576345/1‑149 (MQ=255)
                                                                    cGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata                                                                                 >  2:1282370/1‑149 (MQ=255)
                                                                     gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                                >  1:1639280/1‑149 (MQ=255)
                                                                     gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                                >  2:1330500/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                                                                            <  2:1680436/146‑1 (MQ=255)
                                                                            aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCa                                                                            >  1:1680436/1‑146 (MQ=255)
                                                                             ggCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGa                                                                                                                  >  4:414937/1‑107 (MQ=255)
                                                                             ggCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGa                                                                                                                  <  3:414937/107‑1 (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:1639280/149‑1 (MQ=255)
                                                                                     tGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                                                >  2:596410/1‑149 (MQ=255)
                                                                                        tCCTGCAATATTTTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGACGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                             >  3:1244398/1‑149 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                 <  2:424251/62‑1 (MQ=17)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                 >  1:424251/1‑62 (MQ=17)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                              >  1:672695/1‑113 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                             <  4:188261/146‑1 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                             >  3:188261/1‑146 (MQ=255)
                                                                                           tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                          >  3:660038/1‑149 (MQ=255)
                                                                                            gCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                                                         <  4:1779026/149‑1 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                       >  4:410779/1‑149 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                       <  1:1282370/149‑1 (MQ=255)
                                                                                              aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                       <  3:994083/149‑1 (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                                                                                                                                               >  3:1205169/1‑47 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                               <  4:1205169/47‑1 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                          >  2:677980/1‑100 (MQ=255)
                                                                                                            gCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                          <  1:677980/100‑1 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                       >  3:1694646/1‑149 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                       >  4:171413/1‑149 (MQ=255)
                                                                                                              aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                       >  4:302766/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                                                                             >  3:1510540/1‑97 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                             <  4:1510540/97‑1 (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                           >  4:1875152/1‑135 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAg                           <  3:1875152/135‑1 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt                  >  4:1539939/1‑144 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt                  <  3:1539939/144‑1 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa             >  4:48200/1‑149 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                              <  3:223013/68‑1 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                              >  4:223013/1‑68 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa             >  2:1439802/1‑149 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt           <  4:1809317/149‑1 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt           <  4:1694646/149‑1 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt           <  3:302766/149‑1 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                        >  3:1056245/1‑134 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                        <  4:1056245/134‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                           <  3:77718/96‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                           >  4:77718/1‑96 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc      >  3:421021/1‑149 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc      >  4:1662879/1‑149 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc      <  1:1439802/149‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATACGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc      >  4:628333/1‑149 (MQ=255)
                                                                                                                                                  aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTg   <  3:795736/149‑1 (MQ=255)
                                                                                                                                                  aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTg   <  4:1244398/149‑1 (MQ=255)
                                                                                                                                                  aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGACACCGGCGATATAGTTGGCCCCTGCCAAATAGCCCGTg   <  2:1011875/149‑1 (MQ=255)
                                                                                                                                                   tCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTgg  >  2:34115/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGG  >  NC_000913/610062‑610357

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