Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,160,918 G→C 6.2% C219S (TGC→TCC) ‡ pabC → 4‑amino‑4‑deoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,160,9180GC 6.2% 207.5 / 3.1 64C219S (TGC→TCC) ‡pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

TGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGC  >  NC_002947/2160774‑2161056
                                                                                                                                                |                                                                                                                                          
tGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGAc                                                                                                                                        <  6:169600/149‑1 (MQ=255)
    cGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGt                                                                                                                                    >  2:8504/1‑149 (MQ=255)
      aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACCGGTCTGCGACCTGTCg                                                                                                                                  >  4:261127/1‑149 (MQ=255)
      aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCg                                                                                                                                  >  1:180902/1‑149 (MQ=255)
      aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac                                                                                                                                  <  1:243899/149‑3 (MQ=255)
      aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac                                                                                                                                  <  5:4750/149‑3 (MQ=255)
           gTATACAGAAATCTGTTTCTGGTGCGCGATTGAGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTc                                                                                                                               >  8:130680/1‑147 (MQ=255)
            tataCAGCAATCTGTTTCTTGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGATGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAg                                                                                                                            <  3:71776/149‑1 (MQ=255)
                 aGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCt                                                                                                                       >  5:52893/1‑149 (MQ=255)
                    aaTCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGa                                                                                                                    <  3:284660/149‑1 (MQ=255)
                      tCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGTCGTGGAGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAg                                                                                                                   >  5:245585/1‑148 (MQ=255)
                      tCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAg                                                                                                                   >  3:180882/1‑148 (MQ=255)
                         gTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGTGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAg                                                                                                               >  7:225075/1‑149 (MQ=255)
                                gTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTTTTTCGGTAGAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAAt                                                                                                          <  8:33355/147‑1 (MQ=255)
                                gTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGAc                                                                                                                                        >  3:135749/1‑117 (MQ=255)
                                gTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGAc                                                                                                                                        <  4:135749/117‑1 (MQ=255)
                                 tGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGaa                                                                                                       <  3:7730/149‑1 (MQ=255)
                                   cgcgATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGTGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGt                                                                                                     <  8:225075/149‑1 (MQ=255)
                                            gTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgtttg                                                                                            <  6:249562/149‑1 (MQ=255)
                                                           gACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAttt                                                                                                 >  7:235118/1‑129 (MQ=255)
                                                           gACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAttt                                                                                                 <  8:235118/129‑1 (MQ=255)
                                                             ccTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCa                                                                           >  6:258202/1‑149 (MQ=255)
                                                                 aGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTg                                                                       >  6:317790/1‑149 (MQ=255)
                                                                  gCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTg                                                                       <  5:317790/148‑1 (MQ=255)
                                                                   ccGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgtttg                                                                                            >  6:205554/1‑126 (MQ=255)
                                                                   ccGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgtttg                                                                                            <  5:205554/126‑1 (MQ=255)
                                                                     gTTGTGGCGTGGCGGGCGTGATGCGTGCCGCCCTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGccc                                                                   >  6:247110/1‑149 (MQ=255)
                                                                        gtggcgtggcGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGTGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcgt                                                                 <  6:351105/148‑1 (MQ=255)
                                                                           gcgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTg                                                                                                                                     >  4:22656/1‑77 (MQ=255)
                                                                           gcgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTg                                                                                                                                     <  3:22656/77‑1 (MQ=255)
                                                                           gcgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGCTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcgtgcgt                                                             >  6:225438/1‑149 (MQ=255)
                                                                            cgtgtctGGCGTGATGCGTGCCTCACTGCTGGATAAGGCCCCGCTGGTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGTCGCTGGAGCAAGCGAATTTAGTATTTGTTTTCATTAGCGTTTACGGCATTTTGGCgttgcgtg                                                            <  7:130680/149‑1 (MQ=255)
                                                                            cgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcgtgcgtg                                                            <  4:180882/149‑1 (MQ=255)
                                                                             gtggcgGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcggcgac                                                             >  6:311483/1‑142 (MQ=255)
                                                                                   ggcgTGATGCGTGCCGCACTGCTGGAACAGGCTCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGc                                                     <  6:245585/149‑1 (MQ=255)
                                                                                   ggcgTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGc                                                     <  1:8504/149‑1 (MQ=255)
                                                                                     cgTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCcg                                                   <  3:270441/149‑1 (MQ=255)
                                                                                           gCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgtttg                                                                                            <  4:30778/102‑1 (MQ=255)
                                                                                           gCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgtttg                                                                                            >  3:30778/1‑102 (MQ=255)
                                                                                           gCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTa                                              >  7:254267/1‑148 (MQ=255)
                                                                                                  gCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTc                                      >  7:46465/1‑149 (MQ=255)
                                                                                                      tgctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCcgc                                                  <  7:203569/133‑1 (MQ=255)
                                                                                                      tgctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCcgc                                                  >  8:203569/1‑133 (MQ=255)
                                                                                                        ctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCAtt                                                                         <  2:192347/108‑1 (MQ=255)
                                                                                                        ctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCgg                                 <  5:262596/148‑1 (MQ=255)
                                                                                                               aGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTc                         <  5:120791/149‑1 (MQ=255)
                                                                                                                 gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTCGCCGCGCTAAAc                                           >  6:240082/1‑129 (MQ=255)
                                                                                                                 gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCAc                       >  6:331999/1‑149 (MQ=255)
                                                                                                                 gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCAc                       >  6:65801/1‑149 (MQ=255)
                                                                                                                 gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGAGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAAc                                           <  5:240082/129‑1 (MQ=255)
                                                                                                                   cccGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTc                     <  8:254267/149‑1 (MQ=255)
                                                                                                                   cccGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTc                     <  8:46465/149‑1 (MQ=255)
                                                                                                                      gctgctGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCgg                                 >  4:303323/1‑134 (MQ=255)
                                                                                                                      gctgctGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCgg                                 <  3:303323/134‑1 (MQ=255)
                                                                                                                         gctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAAc               <  6:272151/149‑1 (MQ=255)
                                                                                                                          ctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGt                                       <  2:209512/124‑1 (MQ=255)
                                                                                                                          ctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGt                                       >  1:209512/1‑124 (MQ=255)
                                                                                                                           tgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTg             >  7:5772/1‑149 (MQ=255)
                                                                                                                           tgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTg             >  1:184235/1‑149 (MQ=255)
                                                                                                                             ggTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGc            <  5:65801/148‑1 (MQ=255)
                                                                                                                                    cGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGtt    >  1:58519/1‑149 (MQ=255)
                                                                                                                                    aacgaCAGGTCGCAGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGtt    >  1:266332/5‑149 (MQ=255)
                                                                                                                                      gTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGc  <  2:184235/149‑1 (MQ=255)
                                                                                                                                      cgaCAGGTCGCAGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGc  >  6:118944/3‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
TGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGC  >  NC_002947/2160774‑2161056

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