Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 5.8% 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. 5.8% 796.6 / 37.3 225intergenic (+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 (107/104);  new base . (8/5);  total (115/110)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-02

Only 100 of 228 total aligned reads displayed.

GGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATC  >  NC_000913/119892‑120424
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                
                  gATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATTGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGAGTTGGTGTAAAGATTTTTGCCCTCACCCTGTACGGGGGAGGGCGTAGGGAGATTAATgg                                                                                                                                                                                                                                                   >  4:99521/1‑274 (MQ=255)
                         cATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGCACCGCCTGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAACACCCCGTTCGTTGTCGTGGTGGTGTGTGCACTACTGACGCTTTCCGGGGCGTTGGTGTAAAGCTTATGTCCCTCACCCTGTACAGGTGTAGGCGTAGGTAGATACATGATCTTTTa                                                                                                                                                                                                                                            >  4:1131/1‑275 (MQ=255)
                            ggTGGGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCt                                                                                                                                                                                                                                                                                  <  4:175000/235‑1 (MQ=255)
                                                      tgattCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGGCTGGGGCAGTGCTCTCTGGCGAGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAg                                                                                                                                                                                                               <  2:125213/271‑1 (MQ=255)
                                                                                cGCCGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTCGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAACGCAGCGGTTTCAATGGCAAAGACAACCTCGTTCGTGGTCGTGGTGGTGATCGCACTACTGACGATTTACGGGGCGTTGGCGTAAAGATTATTGCCCTCACCCTGTCCGGGTGAGGGCGTAGAGTGATTAATCCCATTATACGGCTTTGGCGGTTTTCGTCTTCAACCAGTCGGCGATACCCAaaacaaactac                                                                                                                                                                                      >  4:243244/1‑266 (MQ=255)
                                                                                cGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCATCACCCTGTACGGGTGAGGGCGGAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTGTCTCTTTAAAACGATTGCCGATACCTAACACGAACAGCCa                                                                                                                                                                                    >  2:291454/1‑275 (MQ=255)
                                                                                   aGTGTCCGGTGTTAACGCGTGTGGCGCAGTTCTCTCTGGCGCGGGAACCGCAGGTCGCTAAGGGGGAGACGCCGAAGGAAGCGAGTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAc                                                                                                                                                                                 <  4:52011/275‑1 (MQ=255)
                                                                                                                   tctctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTg                                                                                                                                                                                                                                                                                 >  5:224857/1‑149 (MQ=255)
                                                                                                                     tctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTa                                                                                                                                                                                                                                                                               >  5:417349/1‑149 (MQ=255)
                                                                                                                          cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGACGGCAGCGGTTTCAATGGCGAAGCAAACCTCGTTAGTGGTCGTGGTGGTGATTGCACTACTGACGGTTTACGAGTCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                                                                                                                                                                                                            >  1:53827/1‑146 (MQ=255)
                                                                                                                          cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGt                                                                                                                                                                                                                                                                         <  5:936194/149‑1 (MQ=255)
                                                                                                                            cgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGTTGTAAAGATTATTGCCCTCTCCCTGTACGGGTGa                                                                                                                                                                                                                                                                       <  8:446708/149‑1 (MQ=255)
                                                                                                                              tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                                                                                                                                                                                                     <  7:1178641/149‑1 (MQ=255)
                                                                                                                              tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAgg                                                                                                                                                                                                                                                                     <  7:1356843/149‑1 (MQ=255)
                                                                                                                              ggAACCGNATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTg                                                                                                                                                                                                                                                                        <  3:53827/145‑1 (MQ=255)
                                                                                                                               gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                                                                                                                                                                                                    <  6:83185/149‑1 (MQ=255)
                                                                                                                                  ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGAGt                                                                                                                                                                                                                                                                 >  5:50895/1‑149 (MQ=255)
                                                                                                                                    gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                                                                                                                                               >  6:759622/1‑149 (MQ=255)
                                                                                                                                    gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                                                                                                                                               >  6:1104345/1‑149 (MQ=255)
                                                                                                                                               aGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                                                                                                                                                    >  8:809063/1‑149 (MQ=255)
                                                                                                                                               aGGTGGAGAAGCCGAAGGCAGCGGTTTCAATGGCGAAGAAACCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTAccg                                                                                                                      <  2:148635/275‑1 (MQ=255)
                                                                                                                                                ggTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                                                                                                                                                   <  5:787914/149‑1 (MQ=255)
                                                                                                                                                 gTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgcg                                                                                                                                                                                                                                                  >  8:1167511/1‑149 (MQ=255)
                                                                                                                                                   ggAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCt                                                                                                                                                                                                                                                <  8:65572/149‑1 (MQ=255)
                                                                                                                                                     agaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCttt                                                                                                                                                                                                                                              >  6:844755/1‑149 (MQ=255)
                                                                                                                                                         gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                                                                                                                                             >  5:12065/1‑146 (MQ=255)
                                                                                                                                                         gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                                                                                                                                                          >  6:111715/1‑149 (MQ=255)
                                                                                                                                                         gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                                                                                                                                                          >  7:191072/1‑149 (MQ=255)
                                                                                                                                                          ccGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                                                                                                                                                                                         >  7:1275247/1‑149 (MQ=255)
                                                                                                                                                           cGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                                                                                                                                                                                                        <  8:191072/149‑1 (MQ=255)
                                                                                                                                                            gAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCt                                                                                                                                                                                                                                       >  5:830860/1‑149 (MQ=255)
                                                                                                                                                              aGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATCAATGCGCTTTTACGGCttt                                                                                                                                                                                                                                     <  8:362919/149‑1 (MQ=255)
                                                                                                                                                                  aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                                                                                                                                                 >  5:1156021/1‑149 (MQ=255)
                                                                                                                                                                  aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                                                                                                                                                 >  8:1029105/1‑149 (MQ=255)
                                                                                                                                                                   gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgcgc                                                                                                                                                                                                                                                 >  6:215162/1‑132 (MQ=255)
                                                                                                                                                                       tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                                                                                                                                            >  5:841945/1‑149 (MQ=255)
                                                                                                                                                                       tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                                                                                                                                            >  7:585148/1‑149 (MQ=255)
                                                                                                                                                                        ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                                                                                                                                                           <  8:1443109/149‑1 (MQ=255)
                                                                                                                                                                           aaTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                                                                                                                                                                                       <  6:420157/149‑1 (MQ=255)
                                                                                                                                                                             tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                                                                                                                                                     <  7:879906/149‑1 (MQ=255)
                                                                                                                                                                                 gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                 <  7:673542/149‑1 (MQ=255)
                                                                                                                                                                                 gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGAGTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                                                                                                                                                 >  6:94748/1‑149 (MQ=255)
                                                                                                                                                                                  aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                                                                                                                                                                                                <  5:67267/149‑1 (MQ=255)
                                                                                                                                                                                  aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                                                                                                                                                                                                >  8:638321/1‑149 (MQ=255)
                                                                                                                                                                                         ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                                                                                                                                                         >  8:1013624/1‑149 (MQ=255)
                                                                                                                                                                                          cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                                                                                                                                                        >  5:940412/1‑149 (MQ=255)
                                                                                                                                                                                           tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                                                                                                                                                    <  6:1341235/105‑1 (MQ=255)
                                                                                                                                                                                           tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                                                                                                                                                    >  5:1341235/1‑105 (MQ=255)
                                                                                                                                                                                                   ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                                                                                                                                                                                               >  8:471039/1‑149 (MQ=255)
                                                                                                                                                                                                       gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                                                                                                                                           <  6:1450836/149‑1 (MQ=255)
                                                                                                                                                                                                           tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                                                                                                                                                                                       >  6:406869/1‑149 (MQ=255)
                                                                                                                                                                                                                aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa                                                                                                                                                                                  >  6:1400130/1‑149 (MQ=255)
                                                                                                                                                                                                                   gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCGCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCg                                                                                                                                                                               <  6:224857/149‑1 (MQ=255)
                                                                                                                                                                                                                     actactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg                                                                                                                                                                             >  5:63488/1‑149 (MQ=255)
                                                                                                                                                                                                                       tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGAt                                                                                                                                                                           >  8:1079146/1‑149 (MQ=255)
                                                                                                                                                                                                                       tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGAt                                                                                                                                                                           <  5:429805/149‑1 (MQ=255)
                                                                                                                                                                                                                         ctGACGATTTACGGGGGGTTTGTGTAAAGATTATTGCCCTCGCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGTCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAAAACGATCAGCCATACCGGGATCa                                                                                                                                                                         <  8:608823/149‑1 (MQ=255)
                                                                                                                                                                                                                         ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGCTAGCCGATACCTAACACGATCAGCCATACCGGGAACACGTATACCGAAATCGCCATTCCTGGGTTCATCAGCATAATACCCAGTACCGGCTCCATAACCATCAGGAAGATCCAGTTACCCCGGGGATAAACCACAGCAGGGCACCTAGTCACCACGCCtgtt                                            >  4:125240/1‑271 (MQ=255)
                                                                                                                                                                                                                         ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                                                                                                                                                         <  8:509880/149‑1 (MQ=255)
                                                                                                                                                                                                                         ctGACGATGTACGGGGAGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAGTAATGCGCTTTTACGGCTTTGGCGGTTTGCGCTTTAAACAGATTGCCGATACCTAACCCGATCAGCCATACCGGGATCa                                                                                                                                                                         <  6:50895/149‑1 (MQ=255)
                                                                                                                                                                                                                          tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                                                                                                                                                        >  8:339751/1‑149 (MQ=255)
                                                                                                                                                                                                                          tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                                                                                                                                                        >  5:522471/1‑149 (MQ=255)
                                                                                                                                                                                                                          tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                                                                                                                                                        <  7:1013624/149‑1 (MQ=255)
                                                                                                                                                                                                                            aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCCGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCCGCATAATCCCCAGTACCGCCCCCATAAACCAACGGCAGATCCCGTTACCCAGCGGGCAAAGCCAAGACGGACAGCGGGTTACCACGCCTTGTTccc                                        >  4:302388/1‑274 (MQ=255)
                                                                                                                                                                                                                             cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                                                                                                                                                     <  5:1400130/149‑1 (MQ=255)
                                                                                                                                                                                                                              gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                                                                                                                                                    >  8:855504/1‑149 (MQ=255)
                                                                                                                                                                                                                              gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                                                                                                                                                    <  6:645223/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                                                                                                                                                        <  6:997047/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAAATAGCCGATACCTAACACGATCAGCCATACCGGGATc                                                                                                                                                                          >  6:1250161/1‑131 (MQ=255)
                                                                                                                                                                                                                                           ttGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCa                                                                                                                                                       >  6:288454/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCAt                                                                                                                                                      <  5:107325/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCt                                                                                                                                                  >  5:987215/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 tAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGCTACCTAAAACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                                                                                                                                                 >  7:80226/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 tAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                                                                                                                                                 <  7:930415/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 nanaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                                                                                                                                                 >  6:180867/2‑149 (MQ=255)
                                                                                                                                                                                                                                                  aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                                                                                                                                                                                                         <  7:36027/76‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgg                                                                                                                                                <  8:890229/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                                                                                                                                               <  7:855504/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGAAGAGCAGGGAAGCGAGTTACCACGCCATGTTCCTCCTTGGCCCCACGGAATTTCata                 >  2:232855/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    aGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgggg                                                                                                                                              >  5:247010/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                                                                                                                                             >  5:769880/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        tcttGCCCTCACCCCGTACGGGTTAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGTGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATgcgc            <  2:343034/273‑1 (MQ=255)
                                                                                                                                                                                                                                                         tacgCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                                                                                                                                           <  8:759390/96‑1 (MQ=38)
                                                                                                                                                                                                                                                         attGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc                                                                                                                                                                                >  4:376774/1‑110 (MQ=255)
                                                                                                                                                                                                                                                            gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc                                                                                                                                      >  5:185021/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGa                                                                                                                                                                            <  2:376775/110‑1 (MQ=255)
                                                                                                                                                                                                                                                             cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAagcagagc                                                                     <  4:164119/213‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt                                                                                                                                    >  7:644673/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              ccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt                                                                                                                                    >  6:1368856/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              ccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt                                                                                                                                    >  5:219813/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                                                                                                                                                  <  7:299264/70‑1 (MQ=255)
                                                                                                                                                                                                                                                               cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                                                                                                                                                  >  8:299264/1‑70 (MQ=255)
                                                                                                                                                                                                                                                               cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca                                                                                                                                        >  5:1459483/1‑144 (MQ=255)
                                                                                                                                                                                                                                                               cTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa                                                                                                                                   >  8:531837/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                tCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                                                                                                                                                           >  5:322210/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                  aCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATc                                                                                                                                >  7:452954/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGAAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAAGCTAatc  >  2:13884/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    cccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCa                                                                                                                               <  5:1161203/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCa                                                                                                                               >  7:290254/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc                                                                                                                              <  8:1276837/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                
GGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATC  >  NC_000913/119892‑120424

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