Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158:1 +A 7.1% 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.A 7.1% 746.3 / 55.0 183intergenic (+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 . (85/85);  new base A (8/5);  total (93/90)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 184 total aligned reads displayed.

GGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCA  >  NC_000913/120012‑120302
                                                                                                                                                   |                                                                                                                                                
ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACgg                                                                                                                                                <  3:125054/149‑1 (MQ=255)
   gcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTg                                                                                                                                             >  4:837923/1‑149 (MQ=255)
    cgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGTTGTAAAGATTATTGCCCTCTCCCTGTACGGGTGa                                                                                                                                            <  4:446708/149‑1 (MQ=255)
      tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                                                                          <  2:846846/149‑1 (MQ=255)
      tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                                                                          <  3:1178641/149‑1 (MQ=255)
      tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAgg                                                                                                                                          <  1:295999/149‑1 (MQ=255)
      tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAgg                                                                                                                                          <  3:1356843/149‑1 (MQ=255)
       gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                                                                         <  2:83185/149‑1 (MQ=255)
       gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                                                                         <  1:1414669/149‑1 (MQ=255)
          ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGt                                                                                                                                      >  1:1158864/1‑149 (MQ=255)
            gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                    >  2:1104345/1‑149 (MQ=255)
            gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                    >  2:759622/1‑149 (MQ=255)
             cATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTaga                                                                                                                                   <  1:937070/149‑1 (MQ=255)
               tGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagagc                                                                                                                                 <  1:757768/149‑2 (MQ=255)
               tGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagaga                                                                                                                                 <  2:1251543/149‑1 (MQ=255)
                gTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGt                                                                                                                                              >  1:226746/1‑135 (MQ=255)
                  cGATAAGGTGGAGACGCCGAACGCAGCGGTTTCAATGGCGACGAAAACCTCGTTCGTGGTCGAGGTGGTGCTCGCACTACTGACGACTTACGGGGCGTTGGTGTATAGAATATTGCCGTCACCCTGTACGGGTGAGGGCGTAGAGAGAt                                                                                                                              <  3:900254/149‑1 (MQ=255)
                   gATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGa                                                                                                                                            >  3:296574/1‑134 (MQ=255)
                   gATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAtt                                                                                                                             >  4:460694/1‑149 (MQ=255)
                        ggTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                        <  1:787914/149‑1 (MQ=255)
                           ggAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCt                                                                                                                     <  4:65572/149‑1 (MQ=255)
                            gagaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCtt                                                                                                                    <  4:930209/149‑1 (MQ=255)
                             agaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCttt                                                                                                                   >  2:844755/1‑149 (MQ=255)
                                cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                                                                                <  3:621332/149‑1 (MQ=255)
                                 gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                  <  2:12065/146‑1 (MQ=255)
                                 gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                               >  3:191072/1‑149 (MQ=255)
                                 gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                               >  3:645855/1‑149 (MQ=255)
                                 gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                               >  2:111715/1‑149 (MQ=255)
                                  ccGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                                                              >  3:1275247/1‑149 (MQ=255)
                                   cGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGc                                                                                                             <  4:191072/149‑1 (MQ=255)
                                      aGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATCAATGCGCTTTTACGGCttt                                                                                                          <  4:362919/149‑1 (MQ=255)
                                         cAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                       >  1:170926/1‑149 (MQ=255)
                                          aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                      >  1:582716/1‑149 (MQ=255)
                                          aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                      >  4:1029105/1‑149 (MQ=255)
                                           gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgcgc                                                                                                                      >  2:215162/1‑132 (MQ=255)
                                            cGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGt                                                                                                    <  1:904561/149‑1 (MQ=255)
                                               tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                 >  3:931561/1‑149 (MQ=255)
                                               tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                 >  1:510468/1‑149 (MQ=255)
                                                ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                                <  4:1443109/149‑1 (MQ=255)
                                                   aaTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                                                             <  2:420157/149‑1 (MQ=255)
                                                         gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGAGTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                       >  2:94748/1‑149 (MQ=255)
                                                                 ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                               >  4:1013624/1‑149 (MQ=255)
                                                                   tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                         >  1:1341235/1‑105 (MQ=255)
                                                                        cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                                                                        >  2:1013486/1‑149 (MQ=255)
                                                                        cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                                                                        <  3:837923/149‑1 (MQ=255)
                                                                         gtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGACGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                                                                       >  4:864165/1‑149 (MQ=255)
                                                                          tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTaa                                                                      <  1:913838/149‑1 (MQ=255)
                                                                           ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                                                                     >  4:471039/1‑149 (MQ=255)
                                                                               gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                 <  1:471650/149‑1 (MQ=255)
                                                                               gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                 <  2:1450836/149‑1 (MQ=255)
                                                                               gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                 <  3:239655/149‑1 (MQ=255)
                                                                                  gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa                                                              >  2:1359629/1‑149 (MQ=255)
                                                                                   tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                                                             >  2:406869/1‑149 (MQ=255)
                                                                                        aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa                                                        >  2:1400130/1‑149 (MQ=255)
                                                                                           gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCGCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCg                                                     <  2:224857/149‑1 (MQ=255)
                                                                                               tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCtt                                                                                            <  1:981641/106‑1 (MQ=255)
                                                                                                 ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                               <  4:509880/149‑1 (MQ=255)
                                                                                                 ctGACGATGTACGGGGAGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAGTAATGCGCTTTTACGGCTTTGGCGGTTTGCGCTTTAAACAGATTGCCGATACCTAACCCGATCAGCCATACCGGGATCa                                               <  2:50895/149‑1 (MQ=255)
                                                                                                  tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                              >  2:1458445/1‑149 (MQ=255)
                                                                                                  tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                              >  1:522471/1‑149 (MQ=255)
                                                                                                  tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                              <  4:612163/149‑1 (MQ=255)
                                                                                                  tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                              <  3:1013624/149‑1 (MQ=255)
                                                                                                  tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                              >  4:339751/1‑149 (MQ=255)
                                                                                                   gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAgg                                             <  4:110454/149‑1 (MQ=255)
                                                                                                     cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                           <  3:580951/149‑1 (MQ=255)
                                                                                                      gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                          <  2:645223/149‑1 (MQ=255)
                                                                                                       aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtata                                         <  1:1179029/149‑1 (MQ=255)
                                                                                                               ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                 <  4:594529/149‑1 (MQ=255)
                                                                                                               ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                 <  1:837048/149‑1 (MQ=255)
                                                                                                               ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                 <  1:406869/149‑1 (MQ=255)
                                                                                                                  gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATc                                                <  1:1250161/131‑1 (MQ=255)
                                                                                                                  gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                              <  2:997047/149‑1 (MQ=255)
                                                                                                                   ttGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCa                             >  2:288454/1‑149 (MQ=255)
                                                                                                                         tAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                       >  1:1007479/1‑149 (MQ=255)
                                                                                                                         nanaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                       >  2:180867/2‑149 (MQ=255)
                                                                                                                          aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                                                                               <  3:36027/76‑1 (MQ=255)
                                                                                                                           aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                     <  1:1104345/149‑1 (MQ=255)
                                                                                                                           aaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                     >  3:302639/1‑149 (MQ=255)
                                                                                                                            aGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgggg                    >  1:247010/1‑149 (MQ=255)
                                                                                                                             gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                   >  3:702134/1‑149 (MQ=255)
                                                                                                                             gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                   >  1:769880/1‑149 (MQ=255)
                                                                                                                              attattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtc                  <  3:864165/149‑1 (MQ=255)
                                                                                                                               ttattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtca                 >  4:796390/1‑149 (MQ=255)
                                                                                                                                 tacgCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                 >  3:759390/4‑99 (MQ=38)
                                                                                                                                 tacgCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                 <  4:759390/96‑1 (MQ=38)
                                                                                                                                    gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc            >  1:185021/1‑149 (MQ=255)
                                                                                                                                     cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa           >  1:1018510/1‑149 (MQ=255)
                                                                                                                                      ccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt          >  2:1368856/1‑149 (MQ=255)
                                                                                                                                      ccTCACCCTGTACGGATGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTCTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt          >  4:1075957/1‑149 (MQ=255)
                                                                                                                                      ccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                              <  4:1179684/129‑1 (MQ=255)
                                                                                                                                      ccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt          >  1:219813/1‑149 (MQ=255)
                                                                                                                                       cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                        >  4:299264/1‑70 (MQ=255)
                                                                                                                                       cTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa         >  4:531837/1‑149 (MQ=255)
                                                                                                                                        tCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                                <  2:322210/61‑1 (MQ=255)
                                                                                                                                          aCCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATc      >  3:452954/1‑149 (MQ=255)
                                                                                                                                            ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc    <  4:645855/149‑1 (MQ=255)
                                                                                                                                            ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc    <  1:64185/149‑1 (MQ=255)
                                                                                                                                            ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc    <  4:1276837/149‑1 (MQ=255)
                                                                                                                                              tGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCa  <  1:1359629/149‑1 (MQ=255)
                                                                                                                                              tGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCa  >  4:217393/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCA  >  NC_000913/120012‑120302

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