Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,629 T→C 6.5% intergenic (‑10/+72) PP_mr62 ← / ← trxA Pseudomon‑Rho regulatory region/thioredoxin I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,6290TC 6.5% 184.4 / 2.5 64intergenic (‑10/+72)PP_mr62/trxAPseudomon‑Rho regulatory region/thioredoxin I
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGC  >  NC_002947/5948481‑5948774
                                                                                                                                                    |                                                                                                                                                 
aTGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAAt                                                                                                                                                   <  5:54583/149‑1 (MQ=255)
 tGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATg                                                                                                                                                  <  1:177601/149‑1 (MQ=255)
   aCGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGt                                                                                                                                                 >  4:175975/1‑148 (MQ=255)
    cGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCACGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTa                                                                                                                                                >  7:32734/1‑148 (MQ=255)
     gTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGTGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGAAac                                                                                                                                              <  3:186024/149‑1 (MQ=255)
           aTGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTg                                                                                                                                         <  6:52379/148‑1 (MQ=255)
            tGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCt                                                                                                                                       <  3:79720/149‑1 (MQ=255)
                 gaagGTAGAATTAAGCTTCTGGCTAATTGAGTCCGCGCCGCTGGATGGCGACAGGATCGCGTACTTATTCTAATTAGGGATGCTTCGGCTACGGCATGTAGAGGCCACTGGAAAAGAAGTGCGAGGCCGAATGTAACACTTGCtttttt                                                                                                                                  <  7:195111/149‑1 (MQ=255)
                    gggAGAATTTAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGGACTGGAAAAGGAGTGCGAGGCCGAATGTACCAGTTGGTTTTTTCTg                                                                                                                               >  4:60706/1‑149 (MQ=255)
                    gggAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTg                                                                                                                               >  2:142292/1‑149 (MQ=255)
                        gaATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTc                                                                                                                           <  3:175975/149‑1 (MQ=255)
                           ttAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAg                                                                                                                        >  1:47891/1‑149 (MQ=255)
                              aGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGc                                                                                                                     >  5:4224/1‑149 (MQ=255)
                               gCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGct                                                                                                                    >  5:312587/1‑149 (MQ=255)
                                 ttCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGctct                                                                                                                  >  6:23798/1‑149 (MQ=255)
                                      gCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaa                                                                                                              >  1:43097/1‑148 (MQ=255)
                                       cTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaaaa                                                                                                            <  1:236416/149‑1 (MQ=255)
                                       cTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaaaa                                                                                                            <  3:202208/149‑1 (MQ=255)
                                       cTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaa                                                                                                              >  5:94782/1‑147 (MQ=255)
                                       cTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCACGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaaaa                                                                                                            <  8:32734/149‑1 (MQ=255)
                                        ttAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaa                                                                                                              >  6:218766/1‑146 (MQ=255)
                                        ttAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGaaa                                                                                                              >  4:175080/1‑146 (MQ=255)
                                          aaTGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAaaaaga                                                                                                          >  7:31039/1‑148 (MQ=255)
                                          aaTGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAaaaag                                                                                                           <  2:43097/147‑1 (MQ=255)
                                          aaTGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAaaaag                                                                                                           <  5:23798/147‑1 (MQ=255)
                                           aTGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAaaaagaa                                                                                                         <  6:312587/148‑1 (MQ=255)
                                            tGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAaaaagaa                                                                                                         >  4:213758/1‑147 (MQ=255)
                                               ggCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGcc                                                                                                    <  6:4224/149‑1 (MQ=255)
                                               ggCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGc                                                                                                     <  1:176004/148‑1 (MQ=255)
                                                 ccgcgccgcTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGcccc                                                                                                  <  6:211707/149‑1 (MQ=255)
                                                  cgcgccgcTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCg                                                                                                 >  4:174206/1‑149 (MQ=255)
                                                                       gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcaca                                                                            >  3:271285/1‑149 (MQ=255)
                                                                              tACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAgg                                                                     <  3:175080/149‑1 (MQ=255)
                                                                                cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTg                                                                   >  4:30849/1‑149 (MQ=255)
                                                                                   aTTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGc                                                                >  2:69902/1‑149 (MQ=255)
                                                                                    ttCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGc                                                                <  2:47891/148‑1 (MQ=255)
                                                                                             gggATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                  >  2:8320/1‑137 (MQ=255)
                                                                                             gggATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTCGCGGGGCCTTTTCACATCACAGGTgg                                                                  <  1:8320/137‑1 (MQ=255)
                                                                                                     tCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAACGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                              >  4:86837/1‑149 (MQ=255)
                                                                                                     tCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAACGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                              >  4:91204/1‑149 (MQ=255)
                                                                                                       ggCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCt                                            >  7:111328/1‑149 (MQ=255)
                                                                                                       ggCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCt                                            >  2:2341/1‑149 (MQ=255)
                                                                                                              gCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACTGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAttt                                     >  3:55131/1‑149 (MQ=255)
                                                                                                               cATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTg                                    <  3:30849/149‑1 (MQ=255)
                                                                                                                aTGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTg                                    >  7:96781/1‑148 (MQ=255)
                                                                                                                    agagGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAg                               >  7:249726/1‑149 (MQ=255)
                                                                                                                          cACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTGGACAGTGCGcc                         >  6:200582/1‑148 (MQ=255)
                                                                                                                           aCTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGcc                         >  4:94571/1‑148 (MQ=255)
                                                                                                                               gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCt                    >  8:100424/1‑149 (MQ=255)
                                                                                                                                   aGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                <  2:173311/149‑1 (MQ=255)
                                                                                                                                   aGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                <  6:94782/149‑1 (MQ=255)
                                                                                                                                   aGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                <  3:213758/149‑1 (MQ=255)
                                                                                                                                       gTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGt                                                                    >  8:110727/1‑93 (MQ=255)
                                                                                                                                       gTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGt                                                                    <  7:110727/93‑1 (MQ=255)
                                                                                                                                       gTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGGCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCt            <  5:218766/149‑1 (MQ=255)
                                                                                                                                       gTGCGAGGCCGAACGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCATCCTTGGTGGCt            <  3:91204/149‑1 (MQ=255)
                                                                                                                                       gTGCGAGGCCGAACGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCATCCTTGGTGGCt            <  3:86837/149‑1 (MQ=255)
                                                                                                                                        tGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCtt           <  3:174206/149‑1 (MQ=255)
                                                                                                                                          cGAGGCCGAAACTAACACTTGCATTTATCCGCCTCTACTGCTCTGAAAAAGCAAAGCCCCGCATTTGCGGGGCCGTATCACATCACAGGTGCGCGGCGAGGAATGCGGCCAGCTGCGATTTCGACAGTGCGCCAACCTTGGTGGCtt           <  2:167426/147‑1 (MQ=255)
                                                                                                                                               ccGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                >  4:239328/1‑137 (MQ=255)
                                                                                                                                               ccGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                <  3:239328/137‑1 (MQ=255)
                                                                                                                                               ccGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGtt    <  4:271285/149‑1 (MQ=255)
                                                                                                                                               ccGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGtt    >  2:149503/1‑149 (MQ=255)
                                                                                                                                                cGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg   <  1:69902/149‑1 (MQ=255)
                                                                                                                                                 gAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCAGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGc  >  4:230023/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGC  >  NC_002947/5948481‑5948774

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: