Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158 +A 10.7% 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.A10.7% 187.1 / 21.2 56intergenic (+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 . (24/26);  new base A (4/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACC  >  NC_000913/119862‑120454
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                        
tGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGt                                                                                                                                                                                                                                                                                                      >  2:138892/1‑301 (MQ=255)
 ttACTGGATTGGGTGCCGTTTCGTCTTGCGGTTTTGGTTTTTGCCTTGTTCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTg                                                                                                                                                                                                                                                                                                     <  2:292168/300‑1 (MQ=255)
   attGGATTGGGTGCCGGTTCGTCTTGGGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACGTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAg                                                                                                                                                                                                                                                                                                   <  1:471886/299‑1 (MQ=255)
     tGGTTTGGGTGCCGTTTCTTTTTGCTGGGTTGTTATATGCCTTGATCGGTCATGGTGAGACAGCGTTACCGGGCTTGTTTTCTTCGCTTGGTGATTTCCCTAATTCGCAGGATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCTCATGTCGATAAGGTGGCTACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTTGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                                                                                                                                                                                                                                 <  1:365935/301‑1 (MQ=255)
       gtttggtggcgggTTCGTGTTGGGGGTGTGGTTTATGCCTGGATCGGTCATGGTGTGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACGTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCg                                                                                                                                                                                                                                                                                               <  1:464159/290‑1 (MQ=255)
          tggggtcgggTTCGTTGTGTGTGTGTGGTTTTTGCTTTTATCGGTCATGGTGAGATAGCGTTACCGGCCTGGTTTGCTTCGCTGGTTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGGAGTTCTCTCTGGCGCGTGAACCTCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCGCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                                                                                                                                                                            <  1:382398/293‑1 (MQ=255)
            ggTGGCGGTTGGTGTTGTGGTGTTGGTTTTTGCCTTGTTCGGTCATGGTTGGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagag                                                                                                                                                                                                                                                                                          <  1:138767/301‑1 (MQ=255)
              tGCCGGTTCTTCTTGCGGGTGTGGTATTTGTCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGcaataatct                                                                                                                                                                                                                                                                                        <  2:270768/301‑9 (MQ=255)
               gCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGAAGATTTACGGGGCGTAGGTGTAAAGATACATACCCCCACCCGTACAGGGGGAGGGccaacatcaa                                                                                                                                                                                                                                                                                       >  1:270470/1‑292 (MQ=255)
                       gTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTTGTCGTGTTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGAAAAGATTATTTCCCTCACCCTGCACTGGTGAGGGCGTACAAAGATTAAtact                                                                                                                                                                                                                                                                               >  2:47607/1‑298 (MQ=255)
                                gttctgTATATGCCTTGATCGGTCATGTTGGGAAAGCGTTACCGTCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                                                                                                                                                                                                                      <  1:534557/296‑1 (MQ=255)
                                    ggTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGCTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTAACCCTTTAAGGGAGAGGGCGTAGAGAGAATAATGCCCTTTCACCGCatt                                                                                                                                                                                                                                                                  >  1:596285/1‑301 (MQ=255)
                                     gTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACggagatc                                                                                                                                                                                                                                                                 >  2:535230/1‑296 (MQ=255)
                                          tGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCggg                                                                                                                                                                                                                                                            >  2:218024/1‑300 (MQ=255)
                                          tGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCggg                                                                                                                                                                                                                                                            >  2:370345/1‑300 (MQ=255)
                                                       gttggtggggaggggTTACCGGGGTGGGTTGTGTTGGTGGGGGCTTTTCGGACTTGGCCGGAGCAGGGGTTAACTCTGCGGGCGCAGTTATTTCTGGGGGGTGAACCGGATGTGGAGAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAATCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                                               <  1:440287/287‑1 (MQ=255)
                                                       aaTGGTGGGAAATGGTTTCCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                                               <  1:406695/300‑1 (MQ=255)
                                                               gaAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGATTTTACGGCTTTGGCGGGTTTCTCTTTAAACAGATAGcc                                                                                                                                                                                                                                       >  2:539630/1‑301 (MQ=255)
                                                                    ggTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                                                                                                                                                                                                                                  <  2:577661/300‑1 (MQ=255)
                                                                                          gCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTAACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGCGGGCGTAGAGAGCTTAATGCGCTTTTACGCCTTTCGCGGTTTTCTCTTTAAACAACAAGCCCATACCTAACACCATCAGCCAtcacgg                                                                                                                                                                                                            >  1:626031/1‑296 (MQ=255)
                                                                                                         ccgatttcAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGa                                                                                                                                                                                             <  2:73240/294‑1 (MQ=255)
                                                                                                                cAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGaagaacgg                                                                                                                                                                                      >  1:73150/1‑295 (MQ=255)
                                                                                                                  gTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAACAGATACCTAACACGATCAGCCATACCGGGAACAAGTATACCCAACGCGCCAt                                                                                                                                                                                    >  1:491390/1‑301 (MQ=255)
                                                                                                                       aGGTGTTATCGCTTCTGTCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                                                                                                                                                                               <  2:539465/301‑1 (MQ=255)
                                                                                                                             tAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATAGCCATTCCAGGGGtca                                                                                                                                                                         >  2:338926/1‑301 (MQ=255)
                                                                                                                              aaCGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGCTCCATAACACGATCAGCCCTAACGGGATCAGATATAACGACAACGCCATTCCAGTGGacat                                                                                                                                                                        >  1:587509/1‑301 (MQ=255)
                                                                                                                                  gtgttggggtatgtctgggtggtgtgTGAACCGAATGTCGGTAAGGTGTAGTCTTCGAAGGGTGCGGGTTCGATGGCGAAGTATACCTCTTGGGTGGTCGTGGTGGTGATTGCACTACTGACGATTGACGGGGTGTTGTTGTAAAGAGTATTGCCCTCACCTTGTACGGGTGATGGCGTAGAGAGAGTAATGCGCTTTTACGGCGTGGGCGGTTTTCTCTGTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc                                                                                                                                                                    <  1:242313/276‑1 (MQ=255)
                                                                                                                                      tGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCACTTCCTGGGACCATCAACATaa                                                                                                                                                                >  2:183909/1‑301 (MQ=255)
                                                                                                                                      gggtgcagtgtttgtgggggcgTGAACCGCATGGCGTTAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa                                                                                                                                                                <  1:369849/282‑1 (MQ=255)
                                                                                                                                                             agtacggatgtggtttaGGTGGTGACGCGTATGGCAGTGTTTTCAATGGCGAATAAAAACTCGTTCGTGGTCGGGGTGGTGATTGCACTACTGAGGATTGACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATGAATGCGGTTGTACGGCTTAGGCGGTTGTAAGTTTAAATAGATAGCCGATACCTAACACGATCAGGCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAc                                                                                                                                         <  1:538957/285‑1 (MQ=255)
                                                                                                                                                                                 ggAGACGCCGAAGGCAGCTGTTTCAATTGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTAcc                                                                                                                     <  2:344924/301‑1 (MQ=255)
                                                                                                                                                                                      cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAAAGCAGGCAGATCCAATTACCAAGCg                                                                                                                >  2:580709/1‑301 (MQ=255)
                                                                                                                                                                                                            ggggaggtaatctgtgttttgggtggtggtTGTTGTTGCATTTCTGATGTGTTACGGGGCGTTGGTGTAAAGATTATTGGCCTCACCCGGTACGGGTGAGGGGGTAGAGAGAGTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCg                                                                                          <  1:457864/280‑1 (MQ=255)
                                                                                                                                                                                                            ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAACTCGCCATTCCTGGGGTCATCAGCATAATCCCCAGTACCGCCCCCATAAACAGCAGGCAGATCCAGTTACCCCACGGATAAAGCAGAGAAGGGAAACg                                                                                          >  1:339314/1‑301 (MQ=255)
                                                                                                                                                                                                                   ttaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTAcc                                                                                   <  2:571047/299‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTAca                                                                                   >  2:92085/1‑300 (MQ=255)
                                                                                                                                                                                                                     aaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAACCAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCCCCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGAAGAGAAGGGAAGCGAGtacacag                                                                                 >  1:570333/1‑295 (MQ=255)
                                                                                                                                                                                                                                     ttgttggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTg                                                                 <  2:492036/297‑1 (MQ=255)
                                                                                                                                                                                                                                             gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAATTACCCAGCAGATAAAGCACAGCAAGGACACGAGTTACCACCCCCTGTTCTAGCTTGGcgcggcc                                                         >  1:5551/1‑298 (MQ=255)
                                                                                                                                                                                                                                               ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAACCAGAGCAGGGAAGCGAGTTAACACGCCTTGTTCCTACTAGCCCCGACggg                                                       >  1:237896/1‑300 (MQ=255)
                                                                                                                                                                                                                                                ggtccttctGAGGTTTTGCGGGGCGTTGGTGTAAAGATTGTTGCTCTCACCCTGTACGGGTGAGGGCGTAGGGAGATTAATGCGCTTTTACGGCTTTGGCGGTTGTCTCTTTAAACAGATAGCCGATACCTAACACTATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGaa                                                      <  1:419710/294‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGaaa                                                     >  2:274588/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    ttactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTc                                                  <  2:339777/300‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAAAAGCCGGCAGATCCAGTTACCCAGCGGATAAAACAGAGCAGGGAAGCGAGATAACAAGGCTTGTTCCAGCTTGGCACGAAGGAAttac                                                  >  2:141350/1‑299 (MQ=255)
                                                                                                                                                                                                                                                     tactGACTATTAACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCa                                                 <  2:626908/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      attgatgaTTTAGGGGGCGTTGGTGTAAAGATTATTGCCCGCACCCTGTACGGGGGAGGGCGTAGGGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCat                                                <  1:217781/299‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTACATATg                                             >  2:396017/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCAAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATACGCCaaa                                      >  2:115790/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                    gggCGTTGGTGTAATGATTATTGTCTGCACCTTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCATCCATACCGGGATCAGGTATACCGAAATCGCCATTCCGGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGAGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTaa                                  <  1:91961/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCg                            <  2:207167/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTg                    <  2:238163/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAATTAATTACCaa            >  2:54628/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTCCCACGCCTTGTTCCTGCTTGGCTCGACGGAATTTCATAATCCCCAAGCCAACCATCCCCCAGTTGAATACCact           >  1:477229/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                             cTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGTAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATCTGCCCCAGGCTAATCATCGCCAAGTTGATTCCCAggtg         >  1:470680/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                 cccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCACAGTTGATTACCAGTACAGaa     >  2:474281/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    tttACGGGTGAGGGCGTAGAGAGTTTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGAGCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATacc  <  1:274285/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                        
TGTACTGGATTGGGTGCCGGTTCGTCTTGCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACC  >  NC_000913/119862‑120454

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: