Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,160,919 C→A 6.3% C219* (TGC→TGA) ‡ 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,9190CA 6.3% 202.9 / 3.6 64C219* (TGC→TGA) ‡pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base A (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 = 2.72e-01

TGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACGTTAC  >  NC_002947/2160774‑2161067
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                             <  5:4750/149‑3 (MQ=255)
      aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac                                                                                                                                             <  1:243899/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                                                                                                            <  8:235118/129‑1 (MQ=255)
                                                           gACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAttt                                                                                                            >  7:235118/1‑129 (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)
                                                                                                                                                 cGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACGTTAc  <  3:261127/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACGTTAC  >  NC_002947/2160774‑2161067

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: ‑