Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,481,4020GT72.2% 89.0 / 44.4 54intergenic (+53/‑956)unknown/unknownZona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base G (8/7);  new base T (18/21);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGAAGTGTGGCGTATTTTCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGAGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GGAAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCGCTGATGATTT  >  1/1481155‑1481637
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
gatgaAGTGTGGCGTATTTTCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGt                                                                                                                                                                                                                                            >  1:453597/1‑251 (MQ=255)
      gtgtGGCGTATTTTCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGa                                                                                                                                                                                                                                         >  1:422482/1‑248 (MQ=255)
               aTTTTCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAAc                                                                                                                                                                                                                             <  1:744671/251‑1 (MQ=255)
                ttttCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTATAAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAAc                                                                                                                                                                                                                             <  1:342888/251‑1 (MQ=255)
                       gACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGt                                                                                                                                                                                                                                            <  1:122617/228‑1 (MQ=255)
                                                  gTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGtttcacgttgta                                                                                                                                                                                          >  1:492024/1‑242 (MQ=255)
                                                      tttATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGtttcacgttgtaagct                                                                                                                                                                                      <  1:796360/251‑14 (MQ=255)
                                                      tttATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGtttcacgttgtaagct                                                                                                                                                                                      <  1:245385/251‑14 (MQ=255)
                                                                                           cGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCAAAAACAAGATAACTg                                                                                                                                                                                                                   >  1:439398/1‑185 (MQ=255)
                                                                                               cGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGtttcacgttgtaagct                                                                                                                                                                                      <  1:466489/210‑14 (MQ=255)
                                                                                                                       taATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAAc                                                                                                                                                                                                                             <  1:896693/147‑1 (MQ=255)
                                                                                                                                 gTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCTGAAAATAAGATAACTGAAGGCGAAACAGGTTTTACGTTGTTTGCTATAGGATATTTCACGTGAAAACCGAAAAATGATACTAAGTTAGAGATGGTTATTTTATTCCAAGTgtaagttggg                                                                                                           >  1:769522/1‑242 (MQ=255)
                                                                                                                                 gTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGtttcacgttgtaag                                                                                                                                                                                        >  1:211439/1‑163 (MQ=255)
                                                                                                                                              atctaattttgttaagacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTc                                                                                                           <  1:735103/218‑1 (MQ=255)
                                                                                                                                              atctaattttgttaagacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttga                                                                                             >  1:897988/21‑248 (MQ=255)
                                                                                                                                                   attttgttaagacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaag                                                                                        >  1:47361/16‑243 (MQ=255)
                                                                                                                                                             gacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAAc                                                                                                                                                                        <  1:598161/157‑1 (MQ=255)
                                                                                                                                                             gacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                             <  1:356024/184‑1 (MQ=255)
                                                                                                                                                             gacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTa                                                                                                         <  1:645917/220‑1 (MQ=255)
                                                                                                                                                             gacggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatg                                                                              >  1:631333/6‑233 (MQ=255)
                                                                                                                                                                ggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAATAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatgtag                                                                           >  1:92644/3‑230 (MQ=255)
                                                                                                                                                                ggataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                             >  1:126598/3‑186 (MQ=255)
                                                                                                                                                                gTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGAt                                                                                                                                                                                                                        >  1:469852/1‑111 (MQ=255)
                                                                                                                                                                 gataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                   >  1:710266/2‑195 (MQ=255)
                                                                                                                                                                 gataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatgtagg                                                                          >  1:60513/2‑229 (MQ=255)
                                                                                                                                                                  ataCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCAt                                                                                                                                                                                              >  1:216121/1‑135 (MQ=255)
                                                                                                                                                                   taCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                   <  1:573317/193‑1 (MQ=255)
                                                                                                                                                                   taCAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTAc                                                                                                                            >  1:740583/1‑200 (MQ=255)
                                                                                                                                                                        gtgtTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTg                                                                                                                                                                                                                   >  1:501948/1‑108 (MQ=255)
                                                                                                                                                                                  atcTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCaaaa                                                                                                                                                          <  1:858884/155‑1 (MQ=255)
                                                                                                                                                                                         acacTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTaaac                                                   <  1:828510/251‑2 (MQ=255)
                                                                                                                                                                                         acacTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTCTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTaaac                                                   <  1:639394/251‑2 (MQ=255)
                                                                                                                                                                                          cacTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTaaacc                                                  <  1:119271/251‑3 (MQ=255)
                                                                                                                                                                                           acTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTaaacca                                                 <  1:283417/251‑4 (MQ=255)
                                                                                                                                                                                            cTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                             <  1:700016/158‑1 (MQ=255)
                                                                                                                                                                                                t‑aAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTc                                            <  1:21991/251‑1 (MQ=255)
                                                                                                                                                                                                   aaaCGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTaaac                                                   >  1:882242/1‑241 (MQ=255)
                                                                                                                                                                                                     aCGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATg                                        <  1:871823/251‑1 (MQ=255)
                                                                                                                                                                                                      cgcgACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGt                                       <  1:796521/251‑1 (MQ=255)
                                                                                                                                                                                                        cgACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTa                                     <  1:624089/251‑1 (MQ=255)
                                                                                                                                                                                                               gTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCggg                               >  1:315815/1‑250 (MQ=255)
                                                                                                                                                                                                                tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGatat                                                                                                                                       >  1:119379/1‑145 (MQ=255)
                                                                                                                                                                                                                tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGt                                       <  1:92661/241‑1 (MQ=255)
                                                                                                                                                                                                                tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGc                              <  1:251811/250‑1 (MQ=255)
                                                                                                                                                                                                                tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCa                             >  1:321878/1‑251 (MQ=255)
                                                                                                                                                                                                                     gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCaaaaa                                                                                                                                                         >  1:435792/1‑122 (MQ=255)
                                                                                                                                                                                                                     gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGa                                                                                                                                  <  1:472060/145‑1 (MQ=255)
                                                                                                                                                                                                                     gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCAtttt                         >  1:791560/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACt                                                                                                               <  1:595573/163‑1 (MQ=255)
                                                                                                                                                                                                                      aTGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTa                       >  1:18069/1‑251 (MQ=255)
                                                                                                                                                                                                                         tAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAGCCATGCTCGATGTTATTCGGGCATTTTTAttt                    <  1:721810/251‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                   >  1:645560/1‑137 (MQ=255)
                                                                                                                                                                                                                              gtAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAg                                                                                                                                     <  1:602292/133‑1 (MQ=255)
                                                                                                                                                                                                                                           aTGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTgcgcgtacctaaatt  >  1:212832/1‑241 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
GATGAAGTGTGGCGTATTTTCCCGACTAATGAAATTAAATCCAATCATCGGTCTTTTATTGCCGAACATCGTCATTTTACCGATGAAGAAACGGGCGTTTGTTGTGATTTAGTTGTTATTAATCAGTCTGTTGCAGGCATTCCTCGGTTTGTCAAAGACCGTATACAGGTGTTAATCAATCTCCGACACTTAT‑AAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGAGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GGAAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCGCTGATGATTT  >  1/1481155‑1481637

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