Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,453,162 A→T 6.6% C161S (TGC→AGC)  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,1620AT 6.6% 205.2 / 5.2 61C161S (TGC→AGC) 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 = 5.55e-01

ATTGCGCCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCT  >  NC_000913/1453020‑1453309
                                                                                                                                              |                                                                                                                                                   
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:163947/1‑149 (MQ=255)
                   gACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGc                                                                                                                            >  3:4674/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                                                                                                                    >  3:497098/1‑149 (MQ=255)
                           tCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCAc                                                                                                                    >  2:510337/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                                                                                                                                         >  4:223745/1‑102 (MQ=255)
                                                     aGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTaaa                                                                                                                                         <  3:223745/102‑1 (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                                                                              <  3:255552/140‑1 (MQ=255)
                                                                          aGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATc                                                                              >  4:255552/1‑140 (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:510337/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                            <  3:74806/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                            <  1:194159/149‑1 (MQ=255)
                                                                                   gatgatGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAgg                                                            <  1:448047/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                       >  3:156208/1‑149 (MQ=255)
                                                                                                                        aGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAgg                       >  4:517415/1‑149 (MQ=255)
                                                                                                                        aGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAgg                       >  4:274830/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)
                                                                                                                                            gcagGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCgg             <  1:550879/139‑1 (MQ=255)
                                                                                                                                            gcagGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCgg             >  2:550879/1‑139 (MQ=255)
                                                                                                                                             cagGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCt  <  1:101697/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ATTGCGCCTTCGGGAACAAGACCACTATCGACAATTGATTTCACCATCGCCTGAGTCAGTTGGGCCGTCGCGGTAGCTGGTTTGATGATGGCTGGCATTCCGCCCAGCCACGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCT  >  NC_000913/1453020‑1453309

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