Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 8.2% intergenic (+19/+24) ampE → / ← aroP ampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913120,1540T. 8.2% 296.6 / 19.2 85intergenic (+19/+24)ampE/aroPampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base . (5/2);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

CTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCA  >  NC_000913/120008‑120302
                                                                                                                                                  |                                                                                                                                                     
ctctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGt                                                                                                                                                     <  1:889655/149‑1 (MQ=255)
 tctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTa                                                                                                                                                    >  2:1056842/1‑149 (MQ=255)
    ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gg                                                                                                                                                <  3:158563/149‑1 (MQ=255)
    ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gg                                                                                                                                                >  1:1014074/1‑149 (MQ=255)
    ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gg                                                                                                                                                >  1:247110/1‑149 (MQ=255)
      cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                                                              <  3:247290/149‑1 (MQ=255)
      cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                                                              <  1:312838/149‑1 (MQ=255)
        cgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                                                                            <  3:708354/149‑1 (MQ=255)
         gTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                                                                           <  2:680942/149‑1 (MQ=255)
            aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACATCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGTGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGc                                                                                                                                        >  3:144531/1‑149 (MQ=255)
                 cATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAAGCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTaga                                                                                                                                   >  3:98130/1‑149 (MQ=255)
                                  gaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                  <  3:411112/149‑1 (MQ=255)
                                  gaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                  <  1:569271/149‑1 (MQ=255)
                                  gaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                  <  2:247110/149‑1 (MQ=255)
                                  gaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                  <  2:1014074/149‑1 (MQ=255)
                                   aCGCCGAAGGCAGCGGTTTCAACGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTTATGCGCTTTTa                                                                                                                 <  3:652141/149‑1 (MQ=255)
                                     gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAAt                                                                                                                          >  2:279019/1‑138 (MQ=255)
                                     gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAAt                                                                                                                          <  1:279019/138‑1 (MQ=255)
                                     gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAg                                                                                                                                           >  4:416772/1‑121 (MQ=255)
                                     gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAg                                                                                                                                           <  3:416772/121‑1 (MQ=255)
                                        gAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCt                                                                                                            >  3:1112746/1‑149 (MQ=255)
                                        gAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCt                                                                                                            >  1:69814/1‑149 (MQ=255)
                                         aaGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                                                                 >  1:188743/1‑127 (MQ=255)
                                         aaGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                                                                 <  2:188743/127‑1 (MQ=255)
                                                 ggTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                   <  2:872229/149‑1 (MQ=255)
                                                 ggTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                   <  4:1058216/149‑1 (MQ=255)
                                                 ggTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                   <  4:762312/149‑1 (MQ=255)
                                                   tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                 <  2:270765/149‑1 (MQ=255)
                                                      cAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                                                              <  4:144531/149‑1 (MQ=255)
                                                         tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                                                                   <  2:172901/94‑1 (MQ=255)
                                                         tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                                                                   >  1:172901/1‑94 (MQ=255)
                                                         tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                        >  4:280450/1‑120 (MQ=255)
                                                         tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                        <  3:280450/120‑1 (MQ=255)
                                                             gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                       >  2:750394/1‑149 (MQ=255)
                                                             gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                       >  4:1044261/1‑149 (MQ=255)
                                                                 aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAt                                                                                   >  1:714291/1‑149 (MQ=255)
                                                                  aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTCAACAGATa                                                                                  <  4:98130/149‑1 (MQ=255)
                                                                   aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                   <  4:848541/131‑1 (MQ=255)
                                                                   aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                   >  3:848541/1‑131 (MQ=255)
                                                                   aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAg                                                                                 >  3:1111931/1‑149 (MQ=255)
                                                                    aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                                                                                >  2:406222/1‑101 (MQ=255)
                                                                    aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                                                                                <  1:406222/101‑1 (MQ=255)
                                                                      cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                                                                            <  1:805278/135‑1 (MQ=255)
                                                                      cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                                                                            >  2:805278/1‑135 (MQ=255)
                                                                               ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                                                                     >  1:56645/1‑149 (MQ=255)
                                                                                gtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                                                              <  4:422751/123‑1 (MQ=255)
                                                                                gtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                                                              >  3:422751/1‑123 (MQ=255)
                                                                                    tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGAt                                                                <  4:1112746/149‑1 (MQ=255)
                                                                                        ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc                                                            <  3:1044261/149‑1 (MQ=255)
                                                                                        ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc                                                            <  1:1056842/149‑1 (MQ=255)
                                                                                           gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGa                                                                                                                                            <  4:1125066/66‑1 (MQ=38)
                                                                                           gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGa                                                                                                                                            >  3:1125066/1‑66 (MQ=38)
                                                                                            aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAaca                                                                    >  4:487340/1‑137 (MQ=255)
                                                                                            aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAaca                                                                    <  3:487340/137‑1 (MQ=255)
                                                                                                cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                               <  2:18706/90‑1 (MQ=255)
                                                                                                cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                               >  1:18706/1‑90 (MQ=255)
                                                                                                     ctGACGATTTACGGGGCGTTGGAGTAAAGATTATTGCCCTCACCCTGTAC‑GGGGGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                               >  2:711248/1‑149 (MQ=255)
                                                                                                      tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                                                                                <  3:45101/115‑1 (MQ=255)
                                                                                                      tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                                                                                >  4:45101/1‑115 (MQ=255)
                                                                                                       gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                           >  1:198610/1‑103 (MQ=255)
                                                                                                       gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                           <  2:198610/103‑1 (MQ=255)
                                                                                                            tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATAc                                        <  1:750394/149‑1 (MQ=255)
                                                                                                             ttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATAcc                                       <  2:714291/149‑1 (MQ=255)
                                                                                                                 ggggCGTTGATGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGaaa                                   >  4:736321/1‑149 (MQ=255)
                                                                                                                       ttGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCa                             <  2:69814/149‑1 (MQ=255)
                                                                                                                           tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTcc                         >  1:424954/1‑149 (MQ=255)
                                                                                                                            gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                 >  3:977348/1‑140 (MQ=255)
                                                                                                                            gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                 <  4:977348/140‑1 (MQ=255)
                                                                                                                             tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                               >  3:888447/1‑125 (MQ=255)
                                                                                                                             tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                               <  4:888447/125‑1 (MQ=255)
                                                                                                                             tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                       >  3:338538/1‑149 (MQ=255)
                                                                                                                              aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgg                      <  2:56645/149‑1 (MQ=255)
                                                                                                                               aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                           >  2:150202/1‑79 (MQ=255)
                                                                                                                               aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                           <  1:150202/79‑1 (MQ=255)
                                                                                                                                 gATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                   >  1:533394/1‑149 (MQ=255)
                                                                                                                                     attGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGATTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatc               <  3:736321/149‑1 (MQ=255)
                                                                                                                                      ttGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca              >  1:416718/1‑149 (MQ=255)
                                                                                                                                          ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGGTACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt          >  4:158018/1‑149 (MQ=255)
                                                                                                                                          ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                           <  1:499546/116‑1 (MQ=255)
                                                                                                                                          ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                           >  2:499546/1‑116 (MQ=255)
                                                                                                                                          ccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt          >  3:974998/1‑149 (MQ=255)
                                                                                                                                           cTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa         >  1:447797/1‑149 (MQ=255)
                                                                                                                                           cTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa         >  4:73552/1‑149 (MQ=255)
                                                                                                                                                ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc    <  4:338538/149‑1 (MQ=255)
                                                                                                                                                ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTATTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc    <  2:671496/149‑1 (MQ=255)
                                                                                                                                                  tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                        <  4:166365/63‑1 (MQ=255)
                                                                                                                                                  tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                        >  3:166365/1‑63 (MQ=255)
                                                                                                                                                  tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCa  <  4:1111931/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCA  >  NC_000913/120008‑120302

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