Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 610,215 A→G 6.2% intergenic (‑136/+39) 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,2150AG 6.2% 382.9 / 11.1 112intergenic (‑136/+39)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/ferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (58/47);  new base G (4/3);  total (62/50)
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.80e-01

Only 100 of 116 total aligned reads displayed.

CATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGG  >  NC_000913/610086‑610362
                                                                                                                                 |                                                                                                                                                   
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                                                                                                     <  3:1854123/149‑1 (MQ=255)
                                   gATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                                                               >  1:461921/1‑149 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                                                              <  3:1461406/149‑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                                                                                     >  2:1330500/1‑149 (MQ=255)
                                             gCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                                                                     >  3:998570/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                                                                                 >  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)
                                                          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)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                                   >  1:672695/1‑113 (MQ=255)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                                   <  2:672695/113‑1 (MQ=255)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                                  >  3:188261/1‑146 (MQ=255)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTc                                                                  <  4:188261/146‑1 (MQ=255)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                               >  1:886231/1‑149 (MQ=255)
                                                                   tGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                                                               >  3:660038/1‑149 (MQ=255)
                                                                    gCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                                                              <  4:1779026/149‑1 (MQ=255)
                                                                      aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                            <  3:994083/149‑1 (MQ=255)
                                                                      aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                            <  1:1282370/149‑1 (MQ=255)
                                                                      aaTATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtata                                                            >  4:410779/1‑149 (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                                                                                                                                                    <  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                                            >  3:1694646/1‑149 (MQ=255)
                                                                                      aaGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCa                                            >  4:171413/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)
                                                                                               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)
                                                                                                      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)
                                                                                                             gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGcaggggc                     >  4:434713/1‑143 (MQ=255)
                                                                                                                ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAg                                >  4:1875152/1‑135 (MQ=255)
                                                                                                                ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAg                                <  3:1875152/135‑1 (MQ=255)
                                                                                                                ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt                       <  3:1539939/144‑1 (MQ=255)
                                                                                                                ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCt                       >  4:1539939/1‑144 (MQ=255)
                                                                                                                ttCATGCCGGATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa                  >  4:48200/1‑149 (MQ=255)
                                                                                                                ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                                   >  4:223013/1‑68 (MQ=255)
                                                                                                                ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                                                                   <  3:223013/68‑1 (MQ=255)
                                                                                                                  cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt                <  3:302766/149‑1 (MQ=255)
                                                                                                                  cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAAt                <  4:1694646/149‑1 (MQ=255)
                                                                                                                    cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                             >  3:1056245/1‑134 (MQ=255)
                                                                                                                    cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTg                             <  4:1056245/134‑1 (MQ=255)
                                                                                                                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                                >  4:77718/1‑96 (MQ=255)
                                                                                                                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                                                <  3:77718/96‑1 (MQ=255)
                                                                                                                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc           <  1:1439802/149‑1 (MQ=255)
                                                                                                                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc           >  4:1662879/1‑149 (MQ=255)
                                                                                                                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc           >  3:421021/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)
                                                                                                                             cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTc     <  3:1662879/149‑1 (MQ=255)
                                                                                                                             cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTc     <  3:628333/149‑1 (MQ=255)
                                                                                                                              ggCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCt    >  2:31828/1‑149 (MQ=255)
                                                                                                                               gCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTgg                            >  4:588018/1‑124 (MQ=255)
                                                                                                                               gCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTgg                            <  3:588018/124‑1 (MQ=255)
                                                                                                                                cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTgg  >  3:154770/1‑149 (MQ=255)
                                                                                                                                cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTgg  <  4:1496348/149‑1 (MQ=255)
                                                                                                                                cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTgg  >  2:1817635/1‑149 (MQ=255)
                                                                                                                                cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTgg  >  4:987827/1‑149 (MQ=255)
                                                                                                                                 aTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                                                           <  2:183385/91‑1 (MQ=255)
                                                                                                                                 aTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                                                           >  1:183385/1‑91 (MQ=255)
                                                                                                                                 aTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGata                                   <  2:604481/115‑1 (MQ=255)
                                                                                                                                 aTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGata                                   >  1:604481/1‑115 (MQ=255)
                                                                                                                                 |                                                                                                                                                   
CATATTCGAGACTGATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACATGGTCCTGCAATATATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGG  >  NC_000913/610086‑610362

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: