Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 610,202 C→T 7.4% intergenic (‑123/+52) 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,2020CT 7.4% 305.7 / 11.5 94intergenic (‑123/+52)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/ferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (47/40);  new base T (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 = 9.20e-01

GAGGGAGGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAT  >  NC_000913/610056‑610348
                                                                                                                                                  |                                                                                                                                                  
gagggaggTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                  <  3:410744/149‑1 (MQ=255)
      ggTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                            <  3:191082/149‑1 (MQ=255)
         aTGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         <  2:1722784/149‑1 (MQ=21)
                           gACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATGCGGCGTGAAcgcctt                                                                                                                                               >  4:257162/1‑123 (MQ=18)
                              cATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAac                                                                                                                                  <  4:1852906/135‑1 (MQ=255)
                              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                                                                                                                   >  3:1502033/1‑108 (MQ=255)
                                                                        tACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTc                                                                                                                   <  4:1502033/108‑1 (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                                                                       >  3:998570/1‑149 (MQ=255)
                                                                           gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                       >  2:1330500/1‑149 (MQ=255)
                                                                           gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                       >  3:1809317/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                                                                   <  2:1680436/146‑1 (MQ=255)
                                                                                  aGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCa                                                                   >  1:1680436/1‑146 (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                                                       <  4:998570/149‑1 (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)
                                                                                                 tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                     <  2:672695/113‑1 (MQ=255)
                                                                                                 tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                    <  4:188261/146‑1 (MQ=255)
                                                                                                 tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                    >  3:188261/1‑146 (MQ=255)
                                                                                                 tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                 >  1:886231/1‑149 (MQ=255)
                                                                                                 tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                 >  3:660038/1‑149 (MQ=255)
                                                                                                  gCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                                                <  4:1779026/149‑1 (MQ=255)
                                                                                                    aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                              >  4:410779/1‑149 (MQ=255)
                                                                                                    aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                              <  3:994083/149‑1 (MQ=255)
                                                                                                    aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                              <  1:1282370/149‑1 (MQ=255)
                                                                                                     atatatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                          >  3:693917/1‑120 (MQ=255)
                                                                                                     atatatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                                          <  4:693917/120‑1 (MQ=255)
                                                                                                        tatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                     <  1:1176024/106‑1 (MQ=255)
                                                                                                        tatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                     >  2:1176024/1‑106 (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                              >  4:171413/1‑149 (MQ=255)
                                                                                                                    aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                              >  4:302766/1‑149 (MQ=255)
                                                                                                                    aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                              >  3:1694646/1‑149 (MQ=255)
                                                                                                                        tttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCgg                          >  4:1818262/1‑149 (MQ=255)
                                                                                                                        tttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCgg                          >  1:1740174/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                                                                      >  2:1120032/1‑94 (MQ=255)
                                                                                                                                   cGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCAtata                                                                      <  1:1120032/94‑1 (MQ=255)
                                                                                                                                    ggATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTgg              <  2:886231/149‑1 (MQ=255)
                                                                                                                                    ggATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGACACCGGCGATATAGTTgg              >  1:1011875/1‑149 (MQ=255)
                                                                                                                                         aGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt         <  4:386893/149‑1 (MQ=255)
                                                                                                                                          ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                    <  2:1106610/89‑1 (MQ=255)
                                                                                                                                          ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATAcc                                                                    >  1:1106610/1‑89 (MQ=255)
                                                                                                                                          ggCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTg        <  1:596410/149‑1 (MQ=255)
                                                                                                                                           gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGcaggggc       >  4:434713/1‑143 (MQ=255)
                                                                                                                                              ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTa                                                                                                    <  3:1107749/53‑1 (MQ=21)
                                                                                                                                              ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTa                                                                                                    >  4:1107749/1‑53 (MQ=18)
                                                                                                                                              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                                                                                     >  4:223013/1‑68 (MQ=255)
                                                                                                                                              ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                     <  3:223013/68‑1 (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               <  4:1056245/134‑1 (MQ=255)
                                                                                                                                                  cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg               >  3:1056245/1‑134 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAGGGAGGTATGCGTAGTTTTCATATCGACCATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAT  >  NC_000913/610056‑610348

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: