Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,147,209 A→G intergenic (+461/‑560) RPPX_18425 → / → RPPX_18430 enoyl‑CoA hydratase/conjugal transfer protein TraR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,147,2090AG96.2% 75.9 / ‑2.4 53intergenic (+461/‑560)RPPX_18425/RPPX_18430enoyl‑CoA hydratase/conjugal transfer protein TraR
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (26/25);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC  >  CP009974/4147068‑4147348
                                                                                                                                                                  |                                                                                                                                           
ggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGt                                                                                                                                           >  1:4452375/1‑151 (MQ=11)
ggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGt                                                                                                                                           <  1:3509786/151‑1 (MQ=11)
 gCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTg                                                                                                                                          >  1:5355779/1‑151 (MQ=11)
  cGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTgg                                                                                                                                         <  2:989232/151‑1 (MQ=11)
  cGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTgg                                                                                                                                         >  1:2088857/1‑151 (MQ=11)
   gAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgc                                                                                                                                        <  2:5688595/151‑1 (MQ=11)
   gAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgc                                                                                                                                        >  2:3016620/1‑151 (MQ=11)
     aaaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgcgc                                                                                                                                      >  2:6432614/1‑151 (MQ=11)
      aaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgcgct                                                                                                                                     >  1:5793897/1‑150 (MQ=11)
      aaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCa                                                                                                                                     <  2:6498490/151‑1 (MQ=11)
       aGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCAt                                                                                                                                    <  1:6171686/151‑1 (MQ=11)
           gAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTc                                                                                                                                >  2:2679186/1‑151 (MQ=11)
             gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  1:3852058/1‑151 (MQ=11)
             gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  1:3894096/1‑151 (MQ=11)
             gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  2:2054799/1‑151 (MQ=11)
              tCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAcc                                                                                                                             <  1:1906979/151‑1 (MQ=11)
               cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  1:4402424/1‑151 (MQ=11)
               cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  1:4595830/1‑151 (MQ=11)
               cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  2:6497106/1‑151 (MQ=11)
               cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            <  2:1978736/151‑1 (MQ=11)
                 gCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCttt                                                                                                                          >  1:127402/1‑151 (MQ=11)
                  cGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTc                                                                                                                         <  2:2024451/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:6064884/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:2885875/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:3717757/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:3629621/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:6485117/151‑1 (MQ=11)
                   gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  1:1235597/151‑1 (MQ=11)
                    tCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCg                                                                                                                       <  1:1806109/151‑1 (MQ=11)
                     cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:552540/1‑151 (MQ=11)
                     cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:2039499/1‑151 (MQ=11)
                     cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:552404/1‑151 (MQ=11)
                      gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGTGCATAGTCACCTTTCCGcc                                                                                                                     >  1:2302800/1‑151 (MQ=11)
                      gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGcc                                                                                                                     >  1:524568/1‑151 (MQ=11)
                      gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGcc                                                                                                                     >  2:5843614/1‑151 (MQ=11)
                       ccATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGccc                                                                                                                    <  2:4803367/151‑1 (MQ=11)
                        cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   >  2:6481485/1‑151 (MQ=11)
                        cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   <  1:2997509/151‑1 (MQ=11)
                        cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   >  1:1537166/1‑151 (MQ=11)
                         aTCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCtt                                                                                                                  <  1:6419024/151‑1 (MQ=11)
                         aTCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTTCTTGCTCCTTGCTTCTTGGATGTGCGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTCGCGCATAGTAACCTTTCCCCCCtt                                                                                                                  >  2:3816966/1‑151 (MQ=11)
                            gCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcg                                                                                                               <  2:1119856/151‑1 (MQ=11)
                             cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              <  1:333655/151‑1 (MQ=11)
                             cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              <  1:4477192/151‑1 (MQ=11)
                             cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              >  1:689441/1‑151 (MQ=11)
                               aaTGAACTGTTCACTGAGCCAAAACAGACGACTTC‑‑‑‑‑‑‑GCTTCTTGCGTCTTGCTTCATACTTATTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCAATCGCGCAGGTGGCGCATAGTCCCCTTTCCGCCCTTAcggcg                                                                                                            <  2:2304970/151‑1 (MQ=11)
                                  ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         <  1:6091095/151‑1 (MQ=11)
                                  ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         <  2:2720154/151‑1 (MQ=11)
                                  ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         >  1:1116213/1‑151 (MQ=11)
                                   gACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTc                                                                                                        <  1:1056206/151‑1 (MQ=1)
                                   gACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTc                                                                                                        >  1:5323102/1‑151 (MQ=11)
                                                                                                                       cttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACTAAAAAACGCTGGCTCCTATCATCCGGCTCCtgcgctgcgc                                 >  1:831187/1‑151 (MQ=255)
                                                                                                                                                      ttcattcaCCCAAGTGGCGCATAGTCACCTTTCCGCCCGTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTTCGCTCCGGGGTCCCCTCGCCACGGGCTGGCTcac  >  1:3914216/1‑149 (MQ=11)
                                                                                                                                                                  |                                                                                                                                           
GGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC  >  CP009974/4147068‑4147348

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