Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158:1 +A 15.9% intergenic (+23/+20) 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,1581.A15.9% 537.9 / 113.8 170intergenic (+23/+20)ampE/aroPampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (64/79);  new base A (14/13);  total (78/92)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 178 total aligned reads displayed.

TTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCA  >  NC_000913/120060‑120257
                                                                                                   |                                                                                                   
   aaTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACgg                                                                                                   <  2:3562817/98‑1 (MQ=255)
   aaTGGCGAAGAAAACCTCGTTCGTGGTCGGGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGt                                                                                                 <  2:3314691/100‑1 (MQ=255)
    aTGGCGAAGAAAACCTCGTTTGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACggg                                                                                                  <  2:1994815/98‑1 (MQ=255)
      ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAg                                                                                              <  2:546192/100‑1 (MQ=255)
      ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACGACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAg                                                                                              <  1:2076865/100‑1 (MQ=255)
        cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                            >  1:1236150/1‑100 (MQ=255)
        cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                            <  2:3740431/100‑1 (MQ=255)
        cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                            <  1:3621760/100‑1 (MQ=255)
        cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                            <  1:1963357/100‑1 (MQ=255)
         gaagaaAACCTCGTTCGTGGTCGTGGTGGTGGTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGCTTATTGCCCTCACCCGTACAGGGTGAggg                                                                                            <  2:1317484/99‑1 (MQ=255)
         gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGc                                                                                           >  2:3011147/1‑100 (MQ=255)
         gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                           <  1:2919582/100‑1 (MQ=255)
         gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                           <  1:2728926/100‑1 (MQ=255)
          aagaaAACCTCGTTCGTGGTCGGGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGTggg                                                                                            <  2:2505403/98‑1 (MQ=255)
            gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTa                                                                                        <  2:3315072/100‑1 (MQ=255)
            gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGcaat                                                                                        <  2:813576/100‑4 (MQ=38)
              aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATCGCCCTCACCCGTACAGGGTGAGGGcaataa                                                                                      <  1:243652/100‑6 (MQ=37)
               aaCCTCGTTCGTGGTCGTGATGGAGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATAGCCATCACACTGTACGGGTGAGGGCGTagag                                                                                     >  2:1234287/1‑100 (MQ=255)
                aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTACAGATTATTGCCCTCACCCTGTACGGGCGAGGGCGTagaga                                                                                    >  2:1201372/1‑100 (MQ=255)
                aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagaga                                                                                    <  1:2992156/100‑1 (MQ=255)
                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                   <  2:141975/100‑1 (MQ=255)
                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                   >  2:1803541/1‑100 (MQ=255)
                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                   <  1:2182239/100‑1 (MQ=255)
                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                   <  1:630688/100‑1 (MQ=255)
                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                   >  2:3220847/1‑100 (MQ=255)
                  cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagag                                                                                     <  2:105712/97‑1 (MQ=255)
                  cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAgagaga                                                                                  >  1:1562327/1‑100 (MQ=255)
                  cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAgagaga                                                                                  <  1:3759042/100‑1 (MQ=255)
                   tcgctcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAt                                                                                 <  1:2499780/100‑1 (MQ=255)
                     gttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGTCGATGGCGTAAAGTTCCTTGCCCTCACCCTGTACGGGTGGGGGCGTAGAGCGAtt                                                                                >  2:629158/1‑99 (MQ=255)
                        cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                            >  2:3470744/1‑100 (MQ=255)
                          tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgcg                                                                          <  1:1161185/100‑1 (MQ=255)
                                  gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                  >  1:3939615/1‑100 (MQ=255)
                                   tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                 >  2:3168794/1‑100 (MQ=255)
                                        aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                            <  2:2670069/100‑1 (MQ=255)
                                         ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTgg                                                           <  1:1173221/100‑1 (MQ=255)
                                          tGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                          <  2:1696095/100‑1 (MQ=255)
                                            cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGCGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                        >  1:3833034/1‑100 (MQ=255)
                                            cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                        <  1:1803550/100‑1 (MQ=255)
                                              ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                      >  1:2895136/1‑100 (MQ=255)
                                              ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                      >  1:2138221/1‑100 (MQ=255)
                                              ctactGACGATTTACGGGGCGTTGCTGTAAAGATTATTTCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                      <  2:3740740/100‑1 (MQ=255)
                                               tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                                     <  1:308195/100‑1 (MQ=255)
                                               tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                                     <  2:1352010/100‑1 (MQ=255)
                                                 ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                   >  1:1069956/1‑100 (MQ=255)
                                                 ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                   <  2:3774761/100‑1 (MQ=255)
                                                   gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                 <  1:1736537/100‑1 (MQ=255)
                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                               >  1:1086334/1‑100 (MQ=255)
                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                               >  1:432174/1‑100 (MQ=255)
                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                               <  1:1741851/100‑1 (MQ=255)
                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                               <  1:3462545/100‑1 (MQ=255)
                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                               <  1:3476495/100‑1 (MQ=255)
                                                       aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                             <  1:1658593/100‑1 (MQ=255)
                                                        tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                            <  1:398828/100‑1 (MQ=255)
                                                        tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                            >  1:1681838/1‑100 (MQ=255)
                                                          tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                          >  1:3390699/1‑100 (MQ=255)
                                                           aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGGACGGGTGAGGGCGTAGAGAGATTAAAGCGCTTTTACGGCATTGGCGGTTTTCTCTATAAAc                                          <  2:559651/99‑1 (MQ=255)
                                                            cGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAg                                        <  2:3908338/100‑1 (MQ=255)
                                                              gggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAt                                      <  1:2794792/100‑1 (MQ=255)
                                                               ggCGTTGGTGTAAAGATTATTGCCCTCACCCTTTACGGGTGAGTGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGTCGGTTTTCTCTTTAAACAGATa                                     <  2:2138209/100‑1 (MQ=255)
                                                               ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                     <  1:792648/100‑1 (MQ=255)
                                                                gCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAg                                    >  1:2032361/1‑100 (MQ=255)
                                                                  gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                  <  1:3251695/100‑1 (MQ=255)
                                                                  gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                  >  1:1269815/1‑100 (MQ=255)
                                                                   ttGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                 <  2:2204193/100‑1 (MQ=255)
                                                                     ggTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                               >  1:1831899/1‑100 (MQ=255)
                                                                      gtgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                              <  1:3124636/100‑1 (MQ=255)
                                                                      gtgtAAAGATTATTGCCCTCACCCTGTACGGGCGAGGGGGTAGAGAGATTAATGCGCTTTTACGGCTCTCGCGGTTGTCTCTTTGAACAGATAGGCGATa                              >  2:1812530/1‑100 (MQ=255)
                                                                       tgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                               <  2:1942081/98‑1 (MQ=255)
                                                                       tgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                             <  2:228908/100‑1 (MQ=255)
                                                                       tgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                             >  1:1241759/1‑100 (MQ=255)
                                                                        gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAGCAGATAGCCGATAcc                            >  1:2353636/1‑100 (MQ=255)
                                                                        gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                            >  1:2587491/1‑100 (MQ=255)
                                                                        gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                            <  2:2872725/100‑1 (MQ=255)
                                                                         tAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                           >  1:993087/1‑100 (MQ=255)
                                                                          aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                          >  1:960473/1‑100 (MQ=255)
                                                                           aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTaa                         >  1:1878154/1‑100 (MQ=255)
                                                                             gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAaca                       >  1:191607/1‑100 (MQ=255)
                                                                                 tacgCCCTCACCCGTCCAGGGGGAGGGCGTAGTGGGCTCAATGCGCTTCTACGTCTGTGGCGGTATTCTCTTTAAACAGATAGCCACTACCTaagc                       >  2:946771/4‑94 (MQ=255)
                                                                                 attGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGAt                   >  2:2838527/1‑100 (MQ=255)
                                                                                  ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATc                  <  2:1031436/100‑1 (MQ=255)
                                                                                  ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATc                  <  2:621456/100‑1 (MQ=255)
                                                                                    gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                >  2:2552964/1‑100 (MQ=255)
                                                                                    gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                >  1:2049403/1‑100 (MQ=255)
                                                                                    gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                >  2:653632/1‑100 (MQ=255)
                                                                                     cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc               >  1:1857140/1‑100 (MQ=255)
                                                                                     cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGCTACCTACCACGATCAGc               >  1:1710981/1‑100 (MQ=255)
                                                                                     cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc               >  1:1325776/1‑100 (MQ=255)
                                                                                     cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGGTGTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc               >  2:3948238/1‑100 (MQ=255)
                                                                                     accTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                >  2:1894990/2‑99 (MQ=255)
                                                                                       cTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa             >  2:3685559/1‑100 (MQ=255)
                                                                                           cccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGCTCAGCCATAcc         >  1:1876797/1‑100 (MQ=255)
                                                                                            ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCg        >  1:3721140/1‑100 (MQ=255)
                                                                                            ccTGTACGGGTGAGGGCGTAGAGAGATTAATCCGCGTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa           <  2:472738/97‑1 (MQ=255)
                                                                                             cTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCgg       >  1:475147/1‑100 (MQ=255)
                                                                                              tGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg      >  2:2626957/1‑100 (MQ=255)
                                                                                               gTACGGGTGTGGGCGTAGAGAGTTTAAGGAGCTTGCACGGCTTGCGTGGTTTTCTCTTAGAACAGATAGCCGATACCTAACACGATCAGCCATACCgg       <  2:1567176/98‑1 (MQ=255)
                                                                                                tACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGAt    >  2:3802713/1‑100 (MQ=255)
                                                                                                 aCGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATc   <  2:2451804/100‑1 (MQ=255)
                                                                                                  cGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa  <  1:1609594/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCA  >  NC_000913/120060‑120257

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