Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158:1 +A 5.9% 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 5.9% 914.4 / 53.8 221intergenic (+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 . (103/105);  new base A (8/5);  total (111/110)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 223 total aligned reads displayed.

GCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCC  >  NC_000913/119890‑120431
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
gggggTGGGGTATTGGCCTTGTTCGGTCATGGTGGGAAGGGGGTACCGGCCTGGTTGGCTTGGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGAAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                                                                                                                                                                                                              <  4:39302/272‑1 (MQ=255)
       tgGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGACAACCTCGTTCGTGGTCTTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGt                                                                                                                                                                                                                                                                       >  2:130355/1‑274 (MQ=255)
                    gATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGCAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACCGACGATTTACGGGGCGTTGGTGTAACGATTTTTGCCCTCACCCTGTAC‑GGTTGAGGGTGTAGAGAGACTACTgc                                                                                                                                                                                                                                                         >  4:125213/1‑275 (MQ=255)
                                                                                    gaGTTTCAGGTGTTAAGGCGGCTGGCGCAGTTGCTTCTGGTGCGTGAACCGCATGTCGATAAGGGGGGGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGGGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa                                                                                                                                                                                         <  4:130355/274‑1 (MQ=255)
                                                                                      gTATCAGGTGTTAACGGGTCTGGCGCAGGTCTCTCTGGCGCGTGAAACGAATGTCTATAAGGTGGAGACGCCGAAAGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGCGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAAATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACg                                                                                                                                                                                                   <  2:243244/263‑1 (MQ=255)
                                                                                         tCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCAAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg                                                                                                                                                                                    >  2:93758/1‑275 (MQ=255)
                                                                                                                          ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gg                                                                                                                                                                                                                                                                                 <  7:125054/149‑1 (MQ=255)
                                                                                                                            cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGACGGCAGCGGTTTCAATGGCGAAGCAAACCTCGTTAGTGGTCGTGGTGGTGATTGCACTACTGACGGTTTACGAGTCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                                                                                                                                                                                                  >  1:53827/1‑146 (MQ=255)
                                                                                                                            cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                                                                                                                                                                                               <  8:447883/149‑1 (MQ=255)
                                                                                                                            cgcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcat                                                                                                                                                 >  2:337292/1‑275 (MQ=255)
                                                                                                                                tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                                                                                                                                                                                                           <  6:201422/149‑1 (MQ=255)
                                                                                                                                tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                                                                                                                                                                                                           <  6:432277/149‑1 (MQ=255)
                                                                                                                                tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                                                                                                                                                                                                           <  7:1178641/149‑1 (MQ=255)
                                                                                                                                tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                                                                                                                                                                                                           <  8:796019/149‑1 (MQ=255)
                                                                                                                                tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAgg                                                                                                                                                                                                                                                                           <  5:295999/149‑1 (MQ=255)
                                                                                                                                 gAACCGCATGTCGATAATGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAG‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                                                                                                                                                                                                     >  2:369328/1‑202 (MQ=255)
                                                                                                                                 gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                                                                                                                                                                                                          <  6:83185/149‑1 (MQ=255)
                                                                                                                                 gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                                                                                                                                                                                                          <  5:1414669/149‑1 (MQ=255)
                                                                                                                                    ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGt                                                                                                                                                                                                                                                                       <  8:403411/149‑1 (MQ=255)
                                                                                                                                    ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGt                                                                                                                                                                                                                                                                       >  6:107325/1‑149 (MQ=255)
                                                                                                                                    ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGt                                                                                                                                                                                                                                                                       >  5:1158864/1‑149 (MQ=255)
                                                                                                                                    ccGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGAGt                                                                                                                                                                                                                                                                       >  5:50895/1‑149 (MQ=255)
                                                                                                                                     cGCATGTCGNAAAGGTGGAGACGCCNAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                                                                                                                                                                                                                 <  4:369327/202‑1 (MQ=255)
                                                                                                                                      gCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                                                                                                                                                                                                     >  6:759622/1‑149 (MQ=255)
                                                                                                                                       cATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTaga                                                                                                                                                                                                                                                                    <  6:919658/149‑1 (MQ=255)
                                                                                                                                          gTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                                                                                                                                                                                               <  6:226746/135‑1 (MQ=255)
                                                                                                                                            cGATAAGGTGGAGACGCCGAACGCAGCGGTTTCAATGGCGACGAAAACCTCGTTCGTGGTCGAGGTGGTGCTCGCACTACTGACGACTTACGGGGCGTTGGTGTATAGAATATTGCCGTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAt                                                                                                                                                                                                                                                               <  7:900254/149‑1 (MQ=255)
                                                                                                                                             gATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                                                                                                                                                                                                             >  7:296574/1‑134 (MQ=255)
                                                                                                                                             gATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                                                                                                                                                                                                              >  6:64185/1‑149 (MQ=255)
                                                                                                                                                aaGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAAt                                                                                                                                                                                                                                                           <  6:1366648/149‑1 (MQ=255)
                                                                                                                                                 acgTGGAGCCGCCGAAGTCAGCGGTTTCAATGGCGAAGAAAATCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTAAGTGCGTAGAGAGATTAATGCGCTTTTACGGCGTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTAccg                                                                                                                             <  2:142144/272‑1 (MQ=255)
                                                                                                                                                 aGGTGGGGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTAccgc                                                                                                                            <  4:93758/275‑1 (MQ=255)
                                                                                                                                                 aGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                                                                                                                                                                                                          >  8:809063/1‑149 (MQ=255)
                                                                                                                                                  ggTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgc                                                                                                                                                                                                                                                         >  5:784713/1‑149 (MQ=255)
                                                                                                                                                     ggAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCt                                                                                                                                                                                                                                                      <  8:65572/149‑1 (MQ=255)
                                                                                                                                                       agaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCttt                                                                                                                                                                                                                                                    >  6:844755/1‑149 (MQ=255)
                                                                                                                                                          cGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                                                                                                                                                                                                                 <  7:621332/149‑1 (MQ=255)
                                                                                                                                                           gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                                                                                                                                                                >  7:191072/1‑149 (MQ=255)
                                                                                                                                                           gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                                                                                                                                                                                                                >  7:645855/1‑149 (MQ=255)
                                                                                                                                                              gAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCt                                                                                                                                                                                                                                             >  5:830860/1‑149 (MQ=255)
                                                                                                                                                                   cAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                                                                                                                                                        >  5:170926/1‑149 (MQ=255)
                                                                                                                                                                    aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                                                                                                                                                       >  5:582716/1‑149 (MQ=255)
                                                                                                                                                                         tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                                                                                                                                                  >  5:841945/1‑149 (MQ=255)
                                                                                                                                                                         tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                                                                                                                                                                                                  >  7:931561/1‑149 (MQ=255)
                                                                                                                                                                          ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGc                                                                                                                                                                <  4:193138/214‑1 (MQ=255)
                                                                                                                                                                                 gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                                                                                                                                                                                                          >  6:222957/1‑149 (MQ=255)
                                                                                                                                                                                    aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                                                                                                                                                                                                       >  8:638321/1‑149 (MQ=255)
                                                                                                                                                                                        aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                                                                                                                                                                                                   <  8:1361975/149‑1 (MQ=255)
                                                                                                                                                                                           ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                                                                                                                                                                >  8:580951/1‑149 (MQ=255)
                                                                                                                                                                                            cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                                                                                                                                                               >  5:940412/1‑149 (MQ=255)
                                                                                                                                                                                             tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                                                                                                                                                                                                          <  6:1341235/105‑1 (MQ=255)
                                                                                                                                                                                                  cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                                                                                                                                                                                                         <  7:837923/149‑1 (MQ=255)
                                                                                                                                                                                                    tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTaa                                                                                                                                                                                                       <  8:1445465/149‑1 (MQ=255)
                                                                                                                                                                                                    tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTaa                                                                                                                                                                                                       <  5:913838/149‑1 (MQ=255)
                                                                                                                                                                                                     ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                                                                                                                                                                                                      >  8:471039/1‑149 (MQ=255)
                                                                                                                                                                                                         gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                                                                                                                                                  <  7:239655/149‑1 (MQ=255)
                                                                                                                                                                                                         gtggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                                                                                                                                                  <  6:1450836/149‑1 (MQ=255)
                                                                                                                                                                                                                           ctGACGATTTACGGGGGGTTTGTGTAAAGATTATTGCCCTCGCCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGTCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAAAACGATCAGCCATACCGGGATCa                                                                                                                                                                                <  8:608823/149‑1 (MQ=255)
                                                                                                                                                                                                                           ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa                                                                                                                                                                                <  8:509880/149‑1 (MQ=255)
                                                                                                                                                                                                                            tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                                                                                                                                                               <  7:1013624/149‑1 (MQ=255)
                                                                                                                                                                                                                             gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAgg                                                                                                                                                                              <  8:110454/149‑1 (MQ=255)
                                                                                                                                                                                                                              aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCg                                                                                                                                                                                                                                       >  3:70798/1‑91 (MQ=255)
                                                                                                                                                                                                                              aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCCGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCCGCATAATCCCCAGTACCGCCCCCATAAACCAACGGCAGATCCCGTTACCCAGCGGGCAAAGCCAAGACGGACAGCGGGTTACCACGCCTTGTTccc                                               >  4:302388/1‑274 (MQ=255)
                                                                                                                                                                                                                               cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGta                                                                                                                                                                            <  7:580951/149‑1 (MQ=255)
                                                                                                                                                                                                                                gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                                                                                                                                                           >  8:855504/1‑149 (MQ=255)
                                                                                                                                                                                                                                gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                                                                                                                                                           >  7:1276837/1‑149 (MQ=255)
                                                                                                                                                                                                                                 aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtata                                                                                                                                                                          <  5:1179029/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                                                                                                                                                  <  5:222957/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                                                                                                                                                  <  5:406869/149‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                                                                                                                                                               <  6:997047/149‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAAATAGCCGATACCTAACACGATCAGCCATACCGGGATc                                                                                                                                                                                 >  6:1250161/1‑131 (MQ=255)
                                                                                                                                                                                                                                                  gtAACGATTATTGCCCTCACCCTGAAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCt                                                                                                                                                         >  5:1257237/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCt                                                                                                                                                         >  5:987215/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                                                                                                                                                                                                                >  8:36027/1‑76 (MQ=255)
                                                                                                                                                                                                                                                    aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgg                                                                                                                                                       <  8:890229/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                                                                                                                                                      <  5:1104345/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                                                                                                                                                      >  7:302639/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                                                                                                                                                      <  5:759622/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                                                                                                                                                      <  6:1376163/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgggg                                                                                                                                                     >  5:247010/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       gATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGt                                                                                                                                                    >  6:1104906/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        attattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtc                                                                                                                                                   <  7:864165/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        attattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtc                                                                                                                                                   <  6:417349/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           tacgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                                                                                                                                                  >  7:759390/4‑99 (MQ=38)
                                                                                                                                                                                                                                                           attGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAagca                                                                                >  2:164119/1‑213 (MQ=255)
                                                                                                                                                                                                                                                              gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc                                                                                                                                             >  5:185021/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa                                                                                                                                            <  8:1025073/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa                                                                                                                                            >  5:1018510/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAagcagagc                                                                            <  4:164119/213‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccTCACCCTGTAC‑GGATGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTCTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt                                                                                                                                           >  8:1075957/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGcc                                                                                                                                                               <  8:1179684/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaaa                                                                                                                                                                                                                         >  8:299264/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                 cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca                                                                                                                                               <  6:1459483/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa                                                                                                                                          >  5:515702/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                                                                                                                                                                                                 <  6:322210/61‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aCCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATc                                                                                                                                       >  7:452954/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc                                                                                                                                     <  5:64185/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCa                                                                                                                                   <  5:1359629/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCa                                                                                                                                   >  8:217393/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           aC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGccc  <  2:1131/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
GCGGGTGTGGTATATGCCTTGATCGGTCATGGTGAGAAAGCGTTACCGGCCTGGTTTGCTTCGCTGGGTGATTTCCATACTTCGCAGTATCAGGTGTTAACGCGTCTGGCGCAGTTCTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGGGAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCC  >  NC_000913/119890‑120431

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