Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,314 Δ1 bp 5.5% coding (129/1269 nt) umuC → translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,231,3140G. 5.5% 248.0 / 7.6 73coding (129/1269 nt)umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base . (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GTCGTTAAGGCGATGCGCTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGC  >  NC_000913/1231166‑1231456
                                                                                                                                                    |                                                                                                                                              
gtcgtTAAGGCGATGCGCTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAg                                                                                                                                                >  3:58640/1‑149 (MQ=255)
       aGGGGATGCGCTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAgg                                                                                                                                         <  7:288685/149‑1 (MQ=255)
                gcTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCg                                                                                                                                >  7:329857/1‑149 (MQ=255)
                gcTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCg                                                                                                                                >  3:202901/1‑149 (MQ=255)
                  tGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTAGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGt                                                                                                                               >  3:231791/1‑148 (MQ=255)
                    aTGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTa                                                                                                                             <  8:44649/148‑1 (MQ=255)
                          gCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTCTGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAACGCGCTTGGCGTTAAAATgg                                                                                                                       <  5:105317/148‑1 (MQ=255)
                          gCCCTCTGTGATGTAAACGCGGTTTATGCCAGCTGGGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGGTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATgg                                                                                                                       <  7:150675/148‑1 (MQ=255)
                          gCCCTCTGTGATATAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATgg                                                                                                                       <  1:293945/148‑1 (MQ=255)
                             ctctGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGa                                                                                                                   >  7:225849/1‑149 (MQ=255)
                                tgtgATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATcc                                                                                                                >  1:278634/1‑149 (MQ=255)
                                  tgATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCt                                                                                                              >  7:235852/1‑149 (MQ=255)
                                              gTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTcaaacaa                                                                                                   >  2:150476/1‑148 (MQ=255)
                                                tttATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTcaaacaaa                                                                                                  >  5:149831/1‑147 (MQ=255)
                                                      ctaGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCg                                                                                                                                >  1:267051/3‑111 (MQ=255)
                                                      ctaGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCg                                                                                                                                <  2:267051/109‑1 (MQ=255)
                                                      ccAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCt                                                                                           <  1:265880/148‑1 (MQ=255)
                                                       cAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATgg                                                                                                                       >  6:258378/1‑119 (MQ=255)
                                                       cAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATgg                                                                                                                       <  5:258378/119‑1 (MQ=255)
                                                        aGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGtt                                                                                        >  3:135653/1‑149 (MQ=255)
                                                        aGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGtt                                                                                        <  4:135653/149‑1 (MQ=255)
                                                          cTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTtc                                                                                      <  6:56970/149‑1 (MQ=255)
                                                               agaCGTTGGTCCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTTCGTTATCGCCCGAAACGCTGAGGCAAACGCGCTTGGCGTTAAAATTGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGc                                                                                 <  8:156060/149‑1 (MQ=255)
                                                                  cGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCtgt                                                                              <  4:58640/149‑1 (MQ=255)
                                                                       tttCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGt                                                                         >  6:151/1‑149 (MQ=255)
                                                                         tCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGGGGCGTg                                                                        <  6:149831/148‑1 (MQ=255)
                                                                           gCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTcaaa                                                                                                      >  6:122987/1‑116 (MQ=255)
                                                                           gCCCTGATTTATGGGGTAAACCGGTGGTTGGGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTcaaa                                                                                                      <  5:122987/116‑1 (MQ=255)
                                                                             ccTGATTTATGGGGTAATCCGGTGGTTGTGCTATCGAATAATGACGGTTGCGCTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTg                                                                   >  1:52204/1‑149 (MQ=255)
                                                                             ccTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGAt                                                                                             >  5:28070/1‑123 (MQ=255)
                                                                             ccTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGAt                                                                                             <  6:28070/123‑1 (MQ=255)
                                                                              cTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGc                                                                  <  5:151/149‑1 (MQ=255)
                                                                               tGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCt                                                                 >  6:185248/1‑149 (MQ=255)
                                                                                   ttATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTag                                                             >  5:371659/1‑149 (MQ=255)
                                                                                   ttATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGATGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTag                                                             >  8:346791/1‑149 (MQ=255)
                                                                                     aTGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTtc                                                                                      >  1:321084/1‑122 (MQ=255)
                                                                                     aTGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTtc                                                                                      <  2:321084/122‑1 (MQ=255)
                                                                                      tGGGGTAAACCGGTGGTTGTGCTATCGAATCATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTagcag                                                          <  5:186660/149‑1 (MQ=255)
                                                                                      tGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTagcag                                                          >  8:312576/1‑149 (MQ=255)
                                                                                        gggTAAACCGGTGGTTGTGCTATCGAAAAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTagcagc                                                         >  5:174644/1‑148 (MQ=255)
                                                                                          gTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGTCGTGGTTTGCTTTAGCAGCAAt                                                      >  3:352712/1‑149 (MQ=255)
                                                                                             aaCCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTAt                                                   >  7:87295/1‑149 (MQ=255)
                                                                                             aaCCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTAt                                                   <  4:171504/149‑1 (MQ=255)
                                                                                               ccGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGa                                                 <  4:231791/149‑1 (MQ=255)
                                                                                                cGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAg                                                <  4:286312/149‑1 (MQ=255)
                                                                                                cGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAg                                                <  7:312576/149‑1 (MQ=255)
                                                                                                 ggtggtTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGc                                               <  2:278634/149‑1 (MQ=255)
                                                                                                     gtTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGtt                                                                                        >  5:144457/1‑104 (MQ=255)
                                                                                                     gtTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGtt                                                                                        <  6:144457/104‑1 (MQ=255)
                                                                                                      ttGTGCTATCGAATAATGACGGGTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCGTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTAc                                          <  1:319572/149‑1 (MQ=255)
                                                                                                            tATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCtt                                                                >  3:339315/1‑121 (MQ=255)
                                                                                                            tATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCtt                                                                <  4:339315/121‑1 (MQ=255)
                                                                                                            tATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGa                                     >  8:223102/1‑148 (MQ=255)
                                                                                                              tCGAATAATGACGGTTGCGTTATCTCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGCTCTGTTTCGTCTCTGTGGCGTGGTTTGCTTTAGCAGCAATTATCAGCTTTACTCAGACAt                                  >  8:182585/1‑149 (MQ=255)
                                                                                                                gAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAAGAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGa                                <  8:225849/149‑1 (MQ=255)
                                                                                                                gAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGa                                <  2:347045/149‑1 (MQ=255)
                                                                                                                     atGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAAt                           <  1:150476/149‑1 (MQ=255)
                                                                                                                          ggTTGCGTTATCGCCCGAAACGCTGAGTCAAAGGCGATTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTGTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGt                      <  7:182585/149‑1 (MQ=255)
                                                                                                                          ggTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTgg                                                                            >  3:8024/1‑95 (MQ=255)
                                                                                                                          ggTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTgg                                                                            <  4:8024/95‑1 (MQ=255)
                                                                                                                          ggTTGCGGTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTGCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGt                      <  4:202901/149‑1 (MQ=255)
                                                                                                                           gTTGCGCTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTg                     <  2:52204/149‑1 (MQ=255)
                                                                                                                            ttGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGa                    <  8:329857/149‑1 (MQ=255)
                                                                                                                              gCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATg                  >  3:89947/1‑149 (MQ=255)
                                                                                                                                   aTCGCCCGAAACGCTGAGGCAAACGCGCTTGGCGTAAAAATGGGCGATCCCAGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGAc             >  2:286745/1‑149 (MQ=255)
                                                                                                                                      gCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCt          >  4:340546/1‑149 (MQ=255)
                                                                                                                                      gCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGGGATGTCGACGc           <  8:235852/148‑1 (MQ=255)
                                                                                                                                          gAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGaa      <  6:371659/149‑1 (MQ=255)
                                                                                                                                          gAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGaa      <  4:352712/149‑1 (MQ=255)
                                                                                                                                               gCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGc  >  3:294390/1‑148 (MQ=255)
                                                                                                                                               gCTGAGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCgg                        <  6:246097/125‑1 (MQ=255)
                                                                                                                                               gCTGAGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCgg                        <  6:240827/125‑1 (MQ=255)
                                                                                                                                               gCTGAGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCgg                        >  5:246097/1‑125 (MQ=255)
                                                                                                                                               gCTGAGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCgg                        >  5:240827/1‑125 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTCGTTAAGGCGATGCGCTGATGTTTGCCCTCTGTGATGTAAACGCGTTTTATGCCAGCTGTGAGACGGTGTTTCGCCCTGATTTATGGGGTAAACCGGTGGTTGTGCTATCGAATAATGACGGTTGCGTTATCGCCCGAAACGCTGAGGCAAAGGCGCTTGGCGTTAAAATGGGCGATCCCTGGTTCAAACAAAAAGATCTGTTTCGTCGCTGTGGCGTGGTTTGCTTTAGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGC  >  NC_000913/1231166‑1231456

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