Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 19.6% 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.19.6% 606.6 / 142.9 214intergenic (+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 (84/88);  new base . (19/23);  total (103/111)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 220 total aligned reads displayed.

GCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGG  >  NC_000913/120055‑120253
                                                                                                   |                                                                                                    
gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGCTGGTGTAAAGATTATTGCCCTCACCCt                                                                                                      >  2:4210733/1‑100 (MQ=255)
     ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                <  2:2310238/100‑1 (MQ=255)
     ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                <  1:1840980/100‑1 (MQ=255)
      tCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGATGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gg                                                                                               >  2:3115988/1‑100 (MQ=255)
       cAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑ggg                                                                                              <  1:3115970/100‑1 (MQ=255)
          tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                           <  2:1952603/100‑1 (MQ=255)
           ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                          <  1:1365979/100‑1 (MQ=255)
           ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                          <  1:58890/100‑1 (MQ=255)
           ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAg                                                                                          <  2:3583174/100‑1 (MQ=255)
            gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAgg                                                                                         <  1:31755/100‑1 (MQ=255)
             cGAAGAAAGCGTCGTGCGTGGTCGTTGTGGTGATTGCAATACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                        <  1:1577702/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        >  2:3772812/1‑100 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                        <  1:1864757/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                        <  2:3448400/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTCCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        <  2:3224548/100‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCATTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                          <  2:566749/97‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  2:205161/100‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  2:3051588/100‑1 (MQ=255)
               aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCg                                                                                      <  2:8850/100‑1 (MQ=255)
                 gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGcaat                                                                                    <  1:3324880/100‑4 (MQ=38)
                    aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagag                                                                                 >  1:2829299/1‑100 (MQ=255)
                     aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                <  1:1322323/100‑1 (MQ=255)
                     aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                >  1:3479179/1‑100 (MQ=255)
                      ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                               >  2:1589022/1‑100 (MQ=255)
                      ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                               >  1:1996931/1‑100 (MQ=255)
                      ccTCGTTCGGGGTCGTGGTGGTGATTGCACTGCTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                               <  1:531477/100‑1 (MQ=255)
                       cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAgagaga                                                                              <  1:925265/100‑1 (MQ=255)
                         cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTCCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                            >  1:655586/1‑100 (MQ=255)
                         cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                            <  2:146030/100‑1 (MQ=255)
                          nttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTa                                                                           >  2:2240419/2‑100 (MQ=255)
                             cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                        <  2:3479202/100‑1 (MQ=255)
                             cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                        <  2:3303536/100‑1 (MQ=255)
                             cgtggtcgtggtGGTGATTGCACTACTGACGAGTTACGGTGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                        <  1:2991105/100‑1 (MQ=255)
                                ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcgc                                                                     <  1:1069604/100‑1 (MQ=255)
                                ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcgc                                                                     >  1:1457742/1‑100 (MQ=255)
                                     tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCttttn                                                                <  1:2240409/100‑2 (MQ=255)
                                       gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                              >  2:1833333/1‑100 (MQ=255)
                                        tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                             <  2:2513830/100‑1 (MQ=255)
                                         ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                            <  2:2393075/100‑1 (MQ=255)
                                         ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                            >  1:1018772/1‑100 (MQ=255)
                                         ggtgATTGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                            <  1:3061893/100‑1 (MQ=255)
                                            gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCttt                                                         >  1:2067263/1‑100 (MQ=255)
                                               tGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                      >  1:3646690/1‑100 (MQ=255)
                                                gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                     <  2:319432/100‑1 (MQ=255)
                                                gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                     <  2:234968/100‑1 (MQ=255)
                                                gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                     >  1:1165473/1‑100 (MQ=255)
                                                 cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGACCTCACCCTGTAC‑GGGTGAGGGCGTGGAGAGATTAATGCGCTTTTACGGCGTTGGCgg                                                    >  2:3804542/1‑100 (MQ=255)
                                                   ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                  <  2:1969327/100‑1 (MQ=255)
                                                   ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                  <  2:2849111/100‑1 (MQ=255)
                                                    ttctGACGATTTACGGGGTGTTGGTTTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGAGCTTTTACGGCTTTGGCGGt                                                   <  1:2162880/96‑1 (MQ=255)
                                                      ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                               <  2:2590792/100‑1 (MQ=255)
                                                      ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                               <  2:2962743/100‑1 (MQ=255)
                                                        gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                             >  2:767980/1‑100 (MQ=255)
                                                         aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                            <  1:4376758/100‑1 (MQ=255)
                                                           gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                          <  2:2884890/100‑1 (MQ=255)
                                                            aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                         >  1:3906921/1‑100 (MQ=255)
                                                            aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                         >  1:1196122/1‑100 (MQ=255)
                                                             tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                        >  1:1234967/1‑100 (MQ=255)
                                                               tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                      <  1:3389470/100‑1 (MQ=255)
                                                                aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                     <  1:2918770/100‑1 (MQ=255)
                                                                 cGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAg                                    >  2:4088210/1‑100 (MQ=255)
                                                                  ggggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                   <  2:41712/100‑1 (MQ=255)
                                                                  ggggCGTTGGGGTAAAGATTATTGCCCTCGCCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                   <  1:3925319/100‑1 (MQ=255)
                                                                   gggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACACAt                                  >  1:2274761/1‑100 (MQ=255)
                                                                         tGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGa                            <  1:2125744/100‑1 (MQ=255)
                                                                         tGGTGTAAAGATTATTGCCCTCACCCTGGAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGa                            <  2:3733390/100‑1 (MQ=255)
                                                                          ggTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                           >  1:3538640/1‑100 (MQ=255)
                                                                          ggTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                           >  1:3994045/1‑100 (MQ=255)
                                                                            tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                         <  1:3782655/100‑1 (MQ=255)
                                                                            tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                         >  2:2103205/1‑100 (MQ=255)
                                                                            tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTATCTCTTTAAACAGATAGCCGATAc                         >  2:369819/1‑100 (MQ=255)
                                                                             gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                        <  1:1952619/100‑1 (MQ=255)
                                                                             gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                        >  2:2890634/1‑100 (MQ=255)
                                                                             gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                        >  2:4098123/1‑100 (MQ=255)
                                                                              tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                       <  2:3400173/100‑1 (MQ=255)
                                                                              tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                       <  2:966870/100‑1 (MQ=255)
                                                                              tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                       >  1:4279196/1‑100 (MQ=255)
                                                                               aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAgagatcggaag                      >  2:917186/1‑90 (MQ=255)
                                                                               aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                      <  1:124880/100‑1 (MQ=255)
                                                                               aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                      >  2:1430077/1‑100 (MQ=255)
                                                                               aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGCTGTCTCTTTAAACAGATAGCCGATACCTa                      >  2:2561635/1‑100 (MQ=255)
                                                                                 aGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                    >  2:1367733/1‑100 (MQ=255)
                                                                                  gATTATTGGCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGGTTAATGCGCTTTTGCGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                    <  1:3146372/99‑1 (MQ=255)
                                                                                   attattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAacac                  <  2:2919343/100‑1 (MQ=255)
                                                                                        tGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             <  2:2770576/100‑1 (MQ=255)
                                                                                        cgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             >  1:997772/2‑100 (MQ=39)
                                                                                         gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg            >  1:2465118/1‑100 (MQ=255)
                                                                                         gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg            >  1:3084554/1‑100 (MQ=255)
                                                                                          cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  2:349326/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  2:3644635/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  2:1068910/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  1:2761107/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  2:2181607/1‑100 (MQ=255)
                                                                                           ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          >  1:406744/1‑100 (MQ=255)
                                                                                             tCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCAt        >  2:3898737/1‑100 (MQ=255)
                                                                                             tCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCAt        >  1:1342194/1‑100 (MQ=255)
                                                                                              cACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa       <  2:1165477/100‑1 (MQ=255)
                                                                                                cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc     >  2:2582020/1‑100 (MQ=255)
                                                                                                  cTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCgg   >  2:1263163/1‑100 (MQ=255)
                                                                                                  cTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCgg   <  2:1656207/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
GCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGG  >  NC_000913/120055‑120253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑