Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,560,016 T→A 6.1% intergenic (‑53/+90) mcpA ← / ← yhaJ_1 Methyl‑accepting chemotaxis protein McpA/HTH‑type transcriptional regulator YhaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,560,0160TA 6.1% 190.6 / 3.8 66intergenic (‑53/+90)mcpA/yhaJ_1Methyl‑accepting chemotaxis protein McpA/HTH‑type transcriptional regulator YhaJ
Reads supporting (aligned to +/- strand):  ref base T (38/24);  new base A (2/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-02

TAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTG  >  ACB122/2559884‑2560150
                                                                                                                                    |                                                                                                                                      
tccagCGGGACCAGGCCAAATCCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCg                                                                                                                                     <  3:176179/133‑1 (MQ=255)
        gAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGATACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTa                                                                                                                              >  1:136342/1‑135 (MQ=255)
        gAACAGGGCAAATGCCAAAACCACGAAGCCGCAGGCGCTGAGCACGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAAAGAGGGGAGGGCGCAGCCCGCTCATGACGCGGTGCGGGGCCTa                                                                                                                              >  6:323840/1‑135 (MQ=255)
          aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGCCCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATc                                                                                                                            >  2:394820/1‑135 (MQ=255)
          aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATc                                                                                                                            >  7:300174/1‑135 (MQ=255)
          aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATc                                                                                                                            >  5:136086/1‑135 (MQ=255)
          aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCg                                                                                                                           >  7:256065/1‑136 (MQ=255)
           cAGGGCAAATCCCAAACCCCCGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGa                                                                                                                          <  7:72684/136‑1 (MQ=255)
            aGGGCAAATCCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGa                                                                                                                          <  2:66195/135‑1 (MQ=255)
             gggCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAcc                                                                                                                        >  8:69903/1‑136 (MQ=255)
                 aaaTGCCAACACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCGAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAc                                                                                                                         >  2:377002/1‑131 (MQ=255)
                 aaaTGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAcc                                                                                                                        >  4:85123/1‑132 (MQ=255)
                 aaaTGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAcc                                                                                                                        >  4:370117/1‑132 (MQ=255)
                 aaaTGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAc                                                                                                                         >  6:101416/1‑131 (MQ=255)
                         aaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGc                                                                                                             <  1:14715/135‑1 (MQ=255)
                         aaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCt                                                                                                            <  4:43427/136‑1 (MQ=255)
                           acaacgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGCTTCATGTTTCGAGAACTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGAAGCGTTGCGGGGCCAATCGACCGAGATATCGGCtt                                                                                                           >  7:284580/4‑135 (MQ=255)
                           aCCACTACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCt                                                                                                            >  5:120154/1‑134 (MQ=255)
                           aCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCtt                                                                                                           <  6:166266/135‑1 (MQ=255)
                               cgacgCCGCAGGCGCTGAGCAGGAGCTTGTTGCGAAATTTCATGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                      <  3:13407/136‑1 (MQ=255)
                                gacgCCGCAGGCGCTGATCAGGAGCTTGTGGCCAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGAGCTCATGACTAGTTGCTTTGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                      <  2:136342/135‑1 (MQ=255)
                                   gccgcAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGt                                                                                                    >  5:47183/1‑134 (MQ=255)
                                        aGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGCCGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAACAGa                                                                                              >  6:53787/1‑135 (MQ=255)
                                        aGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGa                                                                                              <  5:286915/135‑1 (MQ=255)
                                        aGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGAt                                                                                             >  8:304152/1‑136 (MQ=255)
                                                  cAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGc                                                                                     >  1:387306/1‑134 (MQ=255)
                                                  cAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCt                                                                                    >  4:173623/1‑135 (MQ=255)
                                                  cAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCt                                                                                    >  8:370063/1‑135 (MQ=255)
                                                   aGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGc                                                                                     >  2:12011/1‑133 (MQ=255)
                                                      aTCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAgg                                                                                <  1:9576/135‑1 (MQ=255)
                                                             ggCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGc                                                                                                             <  5:90291/99‑1 (MQ=255)
                                                             ggCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGc                                                                                                             >  6:90291/1‑99 (MQ=255)
                                                             ggCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAAtt                                                                         >  2:293755/1‑135 (MQ=255)
                                                             ggCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAAtt                                                                         >  2:9548/1‑135 (MQ=255)
                                                             ggCGAAATTTCAGGTTCATGTTTCGAGACATTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATTACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTc                                                                        >  8:33662/1‑136 (MQ=255)
                                                               cGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTc                                                                        >  5:272065/1‑134 (MQ=255)
                                                                gAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCtt                                                                      <  1:161072/135‑1 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTa                                                                     >  4:83836/1‑134 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGAGCGGGGCCTATCGACCgag                                                                                                                     >  2:224213/1‑86 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGAGCGGGGCCTATCGACCgag                                                                                                                     <  1:224213/86‑1 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGAGCGGGGCCTATCGACCgag                                                                                                                     <  1:217933/86‑1 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGAGCGGGGCCTATCGACCgag                                                                                                                     >  2:217933/1‑86 (MQ=255)
                                                                  aaTTTCAGGTTCATGTTTCGAGACCGTGAGCCAATGCGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGGCCGAGATATCGGCGTGAGGGGGAAAAGATTGAAGTGCTTCGCGAAAATTCTTag                                                                    >  6:292961/1‑134 (MQ=255)
                                                                         ggTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGttt                                                            <  8:233779/136‑1 (MQ=255)
                                                                                gTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCgacgg                                                      >  2:24855/1‑135 (MQ=255)
                                                                                          cTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGTCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCAc                                           <  8:300174/136‑1 (MQ=255)
                                                                                               gCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGAt                                       >  4:27451/1‑135 (MQ=255)
                                                                                                ccAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGAt                                       >  8:240634/1‑134 (MQ=255)
                                                                                                      aggggagggCGCAGCGCGCTCATGACGCGGTGCGGTGCCTATCGAGCGAGTTATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGTCAGTCACGGATcagccag                                <  8:284580/135‑1 (MQ=255)
                                                                                                      aggggagggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGc                                                                                                             >  8:396614/1‑58 (MQ=255)
                                                                                                      aggggagggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGc                                                                                                             <  7:396614/58‑1 (MQ=255)
                                                                                                       ggggagggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAgct                              <  3:370117/136‑1 (MQ=255)
                                                                                                              gcgcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGa                        >  1:321188/1‑135 (MQ=255)
                                                                                                              gcgcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGAt                       >  8:3372/1‑136 (MQ=255)
                                                                                                              gcgcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGGAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGAt                       <  7:7321/136‑1 (MQ=255)
                                                                                                              gcgcAGCGCGCGCAGGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATGGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGAt                       <  7:13016/136‑1 (MQ=255)
                                                                                                                gcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGa                        >  1:191158/1‑133 (MQ=255)
                                                                                                                gcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGa                        >  6:283030/1‑133 (MQ=255)
                                                                                                                     gcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAgg                >  7:90310/1‑136 (MQ=255)
                                                                                                                      cgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGc               <  7:34199/136‑1 (MQ=255)
                                                                                                                       gcgcATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGc               <  6:180494/132‑1 (MQ=255)
                                                                                                                         tCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGACGTTAGGGAAAATTCTTTCAGACGATTCGACGGACGGCAGTCCCGGATCAGCCAGCTGCTCGATCAGCAGGCGt             >  1:185322/1‑135 (MQ=255)
                                                                                                                           aTGACGCGGTGCGGGGCCAATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGc           <  5:389801/135‑1 (MQ=255)
                                                                                                                                 cgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgc    >  4:94898/1‑136 (MQ=255)
                                                                                                                                  ggTGCGGGGCCTATCGACCGAGATATCGGCGTGAGGTGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTCCAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgc    <  5:310536/135‑1 (MQ=255)
                                                                                                                                   gTGCGGGGCCTATCGACCGATATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGAAGGCGTGCAGCctgctg  <  7:263322/136‑1 (MQ=255)
                                                                                                                                    tGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgctg  <  1:293755/135‑1 (MQ=255)
                                                                                                                                    tGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgctg  >  4:306128/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTG  >  ACB122/2559884‑2560150

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