Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,453,159 A→T 6.6% W162R (TGG→AGG)  paaZ ← oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,453,1590AT 6.6% 203.2 / 5.3 61W162R (TGG→AGG) paaZoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (33/24);  new base T (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATTGCGCCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGG  >  NC_000913/1453020‑1453307
                                                                                                                                           |                                                                                                                                                    
aTTGCGCCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGaa                                                                                                                                             <  2:381731/149‑1 (MQ=255)
   gcgcCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGtt                                                                                                                                          <  4:127307/149‑1 (MQ=255)
       cTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAg                                                                                                                                      <  4:438715/149‑1 (MQ=255)
       cTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGACTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCTGCTGGGGaatgc                                                                                                                                          >  3:514741/1‑142 (MQ=255)
       cTTCGGGAACAAGACCACTATCGACAATTCATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCTGCTGGGGaatgc                                                                                                                                          <  4:514741/145‑4 (MQ=255)
           gggANCAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGt                                                                                                                                  >  1:200716/1‑149 (MQ=255)
              aaCAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTaa                                                                                                                               <  1:341582/149‑1 (MQ=255)
               aCAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAat                                                                                                                              <  1:71560/149‑1 (MQ=255)
                   gACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGc                                                                                                                          >  3:4674/1‑149 (MQ=255)
                   gACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGc                                                                                                                          >  3:163947/1‑149 (MQ=255)
                       aCTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTg                                                                                                                      >  2:149327/1‑149 (MQ=255)
                         tATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGcc                                                                                                                    >  4:491232/1‑149 (MQ=255)
                           tCGACAATTGATTTCACCATTGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGAAGGGGAAGTTAAAGGCGTTAATATGCACTGCCAc                                                                                                                  >  2:448047/1‑149 (MQ=255)
                           tCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCAc                                                                                                                  >  2:510337/1‑149 (MQ=255)
                           tCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCAc                                                                                                                  >  3:497098/1‑149 (MQ=255)
                                      tttCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTg                                                                                                       >  3:256502/1‑149 (MQ=255)
                                       ttCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGa                                                                                                      <  4:163947/149‑1 (MQ=255)
                                          aCCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGt                                                                                                   <  1:149327/149‑1 (MQ=255)
                                            cttCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCa                                                                                                 >  2:474904/3‑149 (MQ=255)
                                             aTCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAg                                                                                                >  1:353942/1‑149 (MQ=255)
                                                 ccTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTaaa                                                                                            >  3:207886/1‑149 (MQ=255)
                                                  cTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAAt                                                                                           <  4:398363/149‑1 (MQ=255)
                                                     aGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTaaa                                                                                                                                       <  3:223745/102‑1 (MQ=255)
                                                     aGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTaaa                                                                                                                                       >  4:223745/1‑102 (MQ=255)
                                                       tCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCg                                                                                                                                   <  3:573665/104‑1 (MQ=255)
                                                       tCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCg                                                                                                                                   >  4:573665/1‑104 (MQ=255)
                                                                ccGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAAt                                                                             >  3:362329/1‑149 (MQ=255)
                                                                    cgcgGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCAc                                                                         <  4:129099/149‑1 (MQ=255)
                                                                     gcgGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCTGCTGGGGaatgct                                                                                                                                         <  1:562174/84‑5 (MQ=255)
                                                                     gcgGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCTGCTGGGGaatgct                                                                                                                                         >  2:562174/1‑80 (MQ=255)
                                                                        gTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACcttc                                                                     <  3:139202/149‑1 (MQ=255)
                                                                         tAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACcttct                                                                    <  2:400758/149‑1 (MQ=255)
                                                                          aGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATc                                                                            >  4:255552/1‑140 (MQ=255)
                                                                          aGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATc                                                                            <  3:255552/140‑1 (MQ=255)
                                                                           gCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCttt                                                                  <  2:140357/149‑1 (MQ=255)
                                                                             tGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCg                                                                >  2:588470/1‑149 (MQ=255)
                                                                               gTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGAt                                                              >  3:326019/1‑149 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                          <  1:194159/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                          <  3:74806/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                          <  1:448047/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                          <  1:510337/149‑1 (MQ=255)
                                                                                    atgatgGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAggg                                                         >  3:64280/1‑149 (MQ=255)
                                                                                    atgatgGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAggg                                                         >  2:93854/1‑149 (MQ=255)
                                                                                            tggcATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAAtt                                                 >  1:502404/1‑149 (MQ=255)
                                                                                                  tCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCtt                                           >  2:231795/1‑149 (MQ=255)
                                                                                                    cGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTcc                                         <  3:395265/149‑1 (MQ=255)
                                                                                                    cGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTcc                                         <  1:442279/149‑1 (MQ=255)
                                                                                                     tgccAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCg                                        <  3:456822/147‑1 (MQ=255)
                                                                                                           ccACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTc                                          >  4:552379/1‑141 (MQ=255)
                                                                                                           ccACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTc                                          <  3:552379/141‑1 (MQ=255)
                                                                                                           ccACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCcaca                                  >  3:28026/1‑149 (MQ=255)
                                                                                                               gTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGt                              <  4:362329/149‑1 (MQ=255)
                                                                                                                      ccAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCa                       >  3:306684/1‑149 (MQ=255)
                                                                                                                       cAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAg                      >  3:511850/1‑149 (MQ=255)
                                                                                                                        aGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAgg                     >  4:274830/1‑149 (MQ=255)
                                                                                                                        aGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAgg                     >  3:156208/1‑149 (MQ=255)
                                                                                                                        aGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAgg                     >  4:517415/1‑149 (MQ=255)
                                                                                                                                cAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTccc             >  3:224220/1‑149 (MQ=255)
                                                                                                                                  gCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCgg           >  4:472510/1‑149 (MQ=255)
                                                                                                                                       ccccAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTAcc      >  3:459281/1‑149 (MQ=255)
                                                                                                                                         ccAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGa    >  1:344453/1‑149 (MQ=255)
                                                                                                                                           agcagGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAgg  <  4:4674/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
ATTGCGCCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGG  >  NC_000913/1453020‑1453307

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