Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 6.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. 6.2% 452.0 / 21.5 146intergenic (+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 (66/70);  new base . (7/2);  total (74/72)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 148 total aligned reads displayed.

TCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCA  >  NC_000913/120005‑120299
                                                                                                                                                     |                                                                                                                                                  
tctctcTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCt                                                                                                                                                    >  2:1132435/1‑150 (MQ=255)
 ctctctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTg                                                                                                                                                   >  1:1673465/1‑150 (MQ=255)
       ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACggg                                                                                                                                            <  1:932837/150‑1 (MQ=255)
       ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCAATACTGACGATTTACGGGGCCTTGGTGTAAAGATTATTTCCCACACACTGTaatgg                                                                                                                                             >  2:113911/1‑146 (MQ=255)
             tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                                                                      <  1:557476/150‑1 (MQ=255)
             tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGAGAggg                                                                                                                                      >  1:1248602/1‑150 (MQ=255)
               aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCg                                                                                                                                    >  2:351339/1‑150 (MQ=255)
                 ccGCATGTCGATTAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTa                                                                                                                                  <  2:1638235/150‑1 (MQ=255)
                  ctcATGTCCATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGGGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                 <  1:13021/148‑1 (MQ=255)
                   gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTaga                                                                                                                                >  2:1058491/1‑150 (MQ=255)
                      tGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagag                                                                                                                             >  2:804481/1‑150 (MQ=255)
                      tGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACATCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACTGGTGAGGGCGTagagag                                                                                                                             >  2:1014827/1‑150 (MQ=255)
                       gTCGATAAGGTGGCGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAgagaga                                                                                                                            <  1:1228464/150‑1 (MQ=255)
                        tCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAt                                                                                                                           <  2:1233208/150‑1 (MQ=255)
                         cGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAtt                                                                                                                          <  2:412689/150‑1 (MQ=255)
                              aGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                     <  2:184559/150‑1 (MQ=255)
                              aGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGGGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                     >  1:350235/1‑150 (MQ=255)
                              aGGTGGAGACGCCGAAGGCAGCGGGTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                     <  1:497254/150‑1 (MQ=255)
                                      tcGCCGAAGGGAGCGGTTTCAATGGCGAAGAAATCCTCGTTAGTGGTCGTGGTGGTGATTGCACTTCTGTCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTTATGCGCTTTTAc                                                                                                             <  1:1145507/149‑1 (MQ=255)
                                      aCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGGCGGGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                                                                             <  1:1497368/150‑1 (MQ=255)
                                       cGCCGATGGCAGCGGTTTGATGGGCGACGAAAACCTCGTGCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTGGGTGTAAATATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTATTGCGCTTTTACg                                                                                                            <  1:158072/150‑1 (MQ=255)
                                       cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                            >  1:764023/1‑150 (MQ=255)
                                       cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                            <  2:1536284/150‑1 (MQ=255)
                                         ccGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                                                                          >  1:414735/1‑150 (MQ=255)
                                         ccGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                                                                          >  2:1273913/1‑150 (MQ=255)
                                            atggCAGCGGTTTCAATGGCGAAGAAATCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCttt                                                                                                       <  1:225646/148‑1 (MQ=255)
                                             aGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                                                                      <  1:380950/150‑1 (MQ=255)
                                             aGGCAGCGGTTTCAATGGCGAAGAAAACCTAGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                                                                      <  2:399849/150‑1 (MQ=255)
                                               gCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                    <  1:1008058/150‑1 (MQ=255)
                                               gCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                    <  2:1277226/150‑1 (MQ=255)
                                                cAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                   >  2:724971/1‑150 (MQ=255)
                                                 aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                                                                  <  2:253763/150‑1 (MQ=255)
                                                  gCGGTTTCAATGGCGCAGAATACCTCGTTCGTGGTCGTGGTGGTTATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGt                                                                                                 <  2:1051246/150‑1 (MQ=255)
                                                   cGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                <  2:1333710/150‑1 (MQ=255)
                                                   cGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGGGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCGTTGGCGGtt                                                                                                <  2:120486/150‑1 (MQ=255)
                                                    ggTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                                                                               >  2:1312180/1‑150 (MQ=255)
                                                      tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                             <  2:158610/150‑1 (MQ=255)
                                                           aTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                        >  1:1528783/1‑150 (MQ=255)
                                                              gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                     <  1:1132435/150‑1 (MQ=255)
                                                                gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                                                                   <  1:844648/150‑1 (MQ=255)
                                                                  agaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTTGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCTTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                                                                 <  1:1014827/150‑1 (MQ=255)
                                                                  agaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                                                                 >  1:1488545/1‑150 (MQ=255)
                                                                   gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAt                                                                                >  2:1470592/1‑150 (MQ=255)
                                                                   gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAt                                                                                >  2:534980/1‑150 (MQ=255)
                                                                   gaaAACCTCGTTCGTGGTCGTGGTGGGGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAt                                                                                <  2:901740/150‑1 (MQ=255)
                                                                       aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                            <  1:431436/150‑1 (MQ=255)
                                                                       aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                            <  1:379037/150‑1 (MQ=255)
                                                                       aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                            >  1:355320/1‑150 (MQ=255)
                                                                        ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                           >  1:503409/1‑150 (MQ=255)
                                                                        ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                           <  2:764023/150‑1 (MQ=255)
                                                                        ccTCGTTCGAGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                           >  2:410820/1‑150 (MQ=255)
                                                                          tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGGACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                                                                         <  2:1045551/150‑1 (MQ=255)
                                                                           cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACACTGTACGGGGGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                                                                        >  1:228987/1‑150 (MQ=255)
                                                                           cgttcgtGGGCGTGGTGGTGATTTCACTACTGACGATTTACGGGGCGTTGGTGTATAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                                                                        <  2:1673465/150‑1 (MQ=255)
                                                                            gttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                                                                       <  1:1058491/150‑1 (MQ=255)
                                                                             ttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                                                                      <  1:1603301/150‑1 (MQ=255)
                                                                              tcgtggtcgtggTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                                                                     <  1:680852/150‑1 (MQ=255)
                                                                               cgtgttcgttgtTGTTTTTGCACTACTTACGATGTACGGGGCGTTGGGGTAAAGATTATTGCACTCACCCTGTACGGGTGAGGGCGTAGATAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                                                                    <  2:97511/150‑1 (MQ=255)
                                                                                 tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                                                                  >  2:378579/1‑150 (MQ=255)
                                                                                          tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc                                                         <  1:378579/150‑1 (MQ=255)
                                                                                             tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCAt                                                      >  2:1001506/1‑150 (MQ=255)
                                                                                             tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGACGATACCTAACACGATCAGCCAt                                                      >  2:351799/1‑150 (MQ=255)
                                                                                                  gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCAATACCgg                                                 >  2:555213/1‑150 (MQ=255)
                                                                                                      tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATTCCTAACACGATCAGCCATACCGGGATc                                             <  1:804481/150‑1 (MQ=255)
                                                                                                          gACGAGTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGGTACGGGTGAGGGTGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGt                                         >  2:71005/1‑150 (MQ=255)
                                                                                                            cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                       <  1:1314715/150‑1 (MQ=255)
                                                                                                            cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTGTGGCGGTTTTCTCTTTTATCAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                       <  1:410820/150‑1 (MQ=255)
                                                                                                                ttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCg                                   <  1:637464/150‑1 (MQ=255)
                                                                                                                        cGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                           <  2:77070/150‑1 (MQ=255)
                                                                                                                             gtgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc                                                   >  2:1423582/1‑121 (MQ=255)
                                                                                                                             gtgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc                                                   <  1:1423582/121‑1 (MQ=255)
                                                                                                                             gtgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTcc                      >  1:1601597/1‑150 (MQ=255)
                                                                                                                               gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                    >  1:448567/1‑150 (MQ=255)
                                                                                                                                 aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATATCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                  >  1:1407472/1‑150 (MQ=255)
                                                                                                                                  aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgggg                 <  1:534980/150‑1 (MQ=255)
                                                                                                                                   aGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGTTTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                <  1:1143216/150‑1 (MQ=255)
                                                                                                                                    gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtc               >  2:683809/1‑150 (MQ=255)
                                                                                                                                      ttattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATGCCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcat             >  2:1097276/1‑150 (MQ=255)
                                                                                                                                       tattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatc            <  1:1660034/150‑1 (MQ=255)
                                                                                                                                        attGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca           <  1:553506/150‑1 (MQ=255)
                                                                                                                                         ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAg          >  2:212124/1‑150 (MQ=255)
                                                                                                                                          cgCCCTCACCCGTACAGCGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCATTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATTCCGAATTCGCGATTCCTGGGGTCTTCAGc         >  2:1159567/2‑150 (MQ=255)
                                                                                                                                           gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa        >  2:1247394/1‑150 (MQ=255)
                                                                                                                                           gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa        >  2:675833/1‑150 (MQ=255)
                                                                                                                                           gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa        >  1:912353/1‑150 (MQ=255)
                                                                                                                                            cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                        >  1:59158/1‑117 (MQ=255)
                                                                                                                                            cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                        <  2:59158/117‑1 (MQ=255)
                                                                                                                                            cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                    <  1:732260/137‑1 (MQ=255)
                                                                                                                                            cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt       >  2:1322172/1‑150 (MQ=255)
                                                                                                                                            cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt       >  2:612401/1‑150 (MQ=255)
                                                                                                                                            cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGTGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt       >  1:657425/1‑150 (MQ=255)
                                                                                                                                            cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt       >  1:679447/1‑150 (MQ=255)
                                                                                                                                            cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt       >  2:1306352/1‑150 (MQ=255)
                                                                                                                                              cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCTGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa     >  2:825220/1‑150 (MQ=255)
                                                                                                                                              cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa     <  2:527681/150‑1 (MQ=255)
                                                                                                                                              cTCACCCTGTACGGGTGACGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa     >  2:558488/1‑150 (MQ=255)
                                                                                                                                               ccACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTTGCGGATTTCTCTTTAATCAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAAt    <  1:113911/149‑1 (MQ=255)
                                                                                                                                                cACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATc   <  2:1076399/150‑1 (MQ=255)
                                                                                                                                                 aCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCa  <  2:1159237/150‑1 (MQ=255)
                                                                                                                                                 aCCCTGTACGGGAGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCa  <  2:1248602/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCA  >  NC_000913/120005‑120299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: