Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 10.9% 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.10.9% 162.2 / 13.7 55intergenic (+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 (24/24);  new base . (4/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGAT  >  NC_000913/119862‑120451
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                          
tGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                                                                                                                                                                                                                   >  2:138892/1‑301 (MQ=255)
 ttACTGGATTGGGTGCCGTTTCGTCTTGCGGTTTTGGTTTTTGCCTTGTTCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                                                                                                                                                                                                                                  <  2:292168/300‑1 (MQ=255)
   attGGATTGGGTGCCGGTTCGTCTTGGGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACGTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                                                                                                                                                                                                                                <  1:471886/299‑1 (MQ=255)
     tGGTTTGGGTGCCGTTTCTTTTTGCTGGGTTGTTATATGCCTTGATCGGTCATGGTGAGACAGCGTTACCGGGCTTGTTTTCTTCGCTTGGTGATTTCCCTAATTCGCAGGATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCTCATGTCGATAAGGTGGCTACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTTGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                                                                                                                                                                                                                              <  1:365935/301‑1 (MQ=255)
       gtttggtggcgggTTCGTGTTGGGGGTGTGGTTTATGCCTGGATCGGTCATGGTGTGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACGTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCg                                                                                                                                                                                                                                                                                            <  1:464159/290‑1 (MQ=255)
          tggggtcgggTTCGTTGTGTGTGTGTGGTTTTTGCTTTTATCGGTCATGGTGAGATAGCGTTACCGGCCTGGTTTGCTTCGCTGGTTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGGAGTTCTCTCTGGCGCGTGAACCTCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCGCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                                                                                                                                                                                                                         <  1:382398/293‑1 (MQ=255)
            ggTGGCGGTTGGTGTTGTGGTGTTGGTTTTTGCCTTGTTCGGTCATGGTTGGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagag                                                                                                                                                                                                                                                                                       <  1:138767/301‑1 (MQ=255)
              tGCCGGTTCTTCTTGCGGGTGTGGTATTTGTCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGcaataatct                                                                                                                                                                                                                                                                                     <  2:270768/301‑9 (MQ=255)
               gCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGAAGATTTACGGGGCGTAGGTGTAAAGATACATACCCCCACCC‑GTACAGGGGGAGGGccaacatcaa                                                                                                                                                                                                                                                                                    >  1:270470/1‑292 (MQ=255)
                       gTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTTGTCGTGTTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGAAAAGATTATTTCCCTCACCCTGCAC‑TGGTGAGGGCGTACAAAGATTAAtact                                                                                                                                                                                                                                                                            >  2:47607/1‑298 (MQ=255)
                                gttctgTATATGCCTTGATCGGTCATGTTGGGAAAGCGTTACCGTCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                                                                                                                                                                                                                   <  1:534557/296‑1 (MQ=255)
                                    ggTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGCTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTAACCCTTTAA‑GGGAGAGGGCGTAGAGAGAATAATGCCCTTTCACCGCatt                                                                                                                                                                                                                                                               >  1:596285/1‑301 (MQ=255)
                                     gTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACggagatc                                                                                                                                                                                                                                                              >  2:535230/1‑296 (MQ=255)
                                          tGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCggg                                                                                                                                                                                                                                                         >  2:370345/1‑300 (MQ=255)
                                          tGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCggg                                                                                                                                                                                                                                                         >  2:218024/1‑300 (MQ=255)
                                                       gttggtggggaggggTTACCGGGGTGGGTTGTGTTGGTGGGGGCTTTTCGGACTTGGCCGGAGCAGGGGTTAACTCTGCGGGCGCAGTTATTTCTGGGGGGTGAACCGGATGTGGAGAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAATCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                                            <  1:440287/287‑1 (MQ=255)
                                                       aaTGGTGGGAAATGGTTTCCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                                            <  1:406695/300‑1 (MQ=255)
                                                               gaAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGATTTTACGGCTTTGGCGGGTTTCTCTTTAAACAGATAGcc                                                                                                                                                                                                                                    >  2:539630/1‑301 (MQ=255)
                                                                    ggTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                                                                                                                                                                                                                               <  2:577661/300‑1 (MQ=255)
                                                                                          gCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTAACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGCGGGCGTAGAGAGCTTAATGCGCTTTTACGCCTTTCGCGGTTTTCTCTTTAAACAACAAGCCCATACCTAACACCATCAGCCAtcacgg                                                                                                                                                                                                         >  1:626031/1‑296 (MQ=255)
                                                                                                         ccgatttcAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGa                                                                                                                                                                                          <  2:73240/294‑1 (MQ=255)
                                                                                                                cAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGaagaacgg                                                                                                                                                                                   >  1:73150/1‑295 (MQ=255)
                                                                                                                  gTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAACAGATACCTAACACGATCAGCCATACCGGGAACAAGTATACCCAACGCGCCAt                                                                                                                                                                                 >  1:491390/1‑301 (MQ=255)
                                                                                                                       aGGTGTTATCGCTTCTGTCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                                                                                                                                                                            <  2:539465/301‑1 (MQ=255)
                                                                                                                             tAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATAGCCATTCCAGGGGtca                                                                                                                                                                      >  2:338926/1‑301 (MQ=255)
                                                                                                                              aaCGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGCTCCATAACACGATCAGCCCTAACGGGATCAGATATAACGACAACGCCATTCCAGTGGacat                                                                                                                                                                     >  1:587509/1‑301 (MQ=255)
                                                                                                                                  gtgttggggtatgtctgggtggtgtgTGAACCGAATGTCGGTAAGGTGTAGTCTTCGAAGGGTGCGGGTTCGATGGCGAAGTATACCTCTTGGGTGGTCGTGGTGGTGATTGCACTACTGACGATTGACGGGGTGTTGTTGTAAAGAGTATTGCCCTCACCTTGTAC‑GGGTGATGGCGTAGAGAGAGTAATGCGCTTTTACGGCGTGGGCGGTTTTCTCTGTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc                                                                                                                                                                 <  1:242313/276‑1 (MQ=255)
                                                                                                                                      tGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCACTTCCTGGGACCATCAACATaa                                                                                                                                                             >  2:183909/1‑301 (MQ=255)
                                                                                                                                      gggtgcagtgtttgtgggggcgTGAACCGCATGGCGTTAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa                                                                                                                                                             <  1:369849/282‑1 (MQ=255)
                                                                                                                                                             agtacggatgtggtttaGGTGGTGACGCGTATGGCAGTGTTTTCAATGGCGAATAAAAACTCGTTCGTGGTCGGGGTGGTGATTGCACTACTGAGGATTGACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATGAATGCGGTTGTACGGCTTAGGCGGTTGTAAGTTTAAATAGATAGCCGATACCTAACACGATCAGGCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAc                                                                                                                                      <  1:538957/285‑1 (MQ=255)
                                                                                                                                                                                 ggAGACGCCGAAGGCAGCTGTTTCAATTGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTAcc                                                                                                                  <  2:344924/301‑1 (MQ=255)
                                                                                                                                                                                      cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAAAGCAGGCAGATCCAATTACCAAGCg                                                                                                             >  2:580709/1‑301 (MQ=255)
                                                                                                                                                                                                            ggggaggtaatctgtgttttgggtggtggtTGTTGTTGCATTTCTGATGTGTTACGGGGCGTTGGTGTAAAGATTATTGGCCTCACCCGGTAC‑GGGTGAGGGGGTAGAGAGAGTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCg                                                                                       <  1:457864/280‑1 (MQ=255)
                                                                                                                                                                                                            ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAACTCGCCATTCCTGGGGTCATCAGCATAATCCCCAGTACCGCCCCCATAAACAGCAGGCAGATCCAGTTACCCCACGGATAAAGCAGAGAAGGGAAACg                                                                                       >  1:339314/1‑301 (MQ=255)
                                                                                                                                                                                                                   ttaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTAcc                                                                                <  2:571047/299‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTAca                                                                                >  2:92085/1‑300 (MQ=255)
                                                                                                                                                                                                                     aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAACCAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCCCCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGAAGAGAAGGGAAGCGAGtacacag                                                                              >  1:570333/1‑295 (MQ=255)
                                                                                                                                                                                                                                     ttgttggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTg                                                              <  2:492036/297‑1 (MQ=255)
                                                                                                                                                                                                                                             gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAATTACCCAGCAGATAAAGCACAGCAAGGACACGAGTTACCACCCCCTGTTCTAGCTTGGcgcggcc                                                      >  1:5551/1‑298 (MQ=255)
                                                                                                                                                                                                                                               ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAACCAGAGCAGGGAAGCGAGTTAACACGCCTTGTTCCTACTAGCCCCGACggg                                                    >  1:237896/1‑300 (MQ=255)
                                                                                                                                                                                                                                                ggtccttctGAGGTTTTGCGGGGCGTTGGTGTAAAGATTGTTGCTCTCACCCTGTAC‑GGGTGAGGGCGTAGGGAGATTAATGCGCTTTTACGGCTTTGGCGGTTGTCTCTTTAAACAGATAGCCGATACCTAACACTATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGaa                                                   <  1:419710/294‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGaaa                                                  >  2:274588/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    ttactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTc                                               <  2:339777/300‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAAAGCCGGCAGATCCAGTTACCCAGCGGATAAAACAGAGCAGGGAAGCGAGATAACAAGGCTTGTTCCAGCTTGGCACGAAGGAAttac                                               >  2:141350/1‑299 (MQ=255)
                                                                                                                                                                                                                                                     tactGACTATTAACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCa                                              <  2:626908/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      attgatgaTTTAGGGGGCGTTGGTGTAAAGATTATTGCCCGCACCCTGTAC‑GGGGGAGGGCGTAGGGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCat                                             <  1:217781/299‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTACATATg                                          >  2:396017/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCAAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATACGCCaaa                                   >  2:115790/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                    gggCGTTGGTGTAATGATTATTGTCTGCACCTTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCATCCATACCGGGATCAGGTATACCGAAATCGCCATTCCGGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGAGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTaa                               <  1:91961/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                         <  2:207167/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTg                 <  2:238163/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAATTAATTACCaa         >  2:54628/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTCCCACGCCTTGTTCCTGCTTGGCTCGACGGAATTTCATAATCCCCAAGCCAACCATCCCCCAGTTGAATACCact        >  1:477229/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                             cTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGTAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATCTGCCCCAGGCTAATCATCGCCAAGTTGATTCCCAggtg      >  1:470680/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                 cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCACAGTTGATTACCAGTACAGaa  >  2:474281/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                          
TGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGAT  >  NC_000913/119862‑120451

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