Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,361,714 C→A R199S (CGC→AGC)  puuD → gamma‑glutamyl‑gamma‑aminobutyrate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,361,7140CA98.6% 517.2 / ‑0.5 145R199S (CGC→AGC) puuDgamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (68/75);  total (69/76)
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.81e-01

Only 100 of 154 total aligned reads displayed.

TGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTT  >  NC_000913/1361578‑1361849
                                                                                                                                        |                                                                                                                                       
 gCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAgctagctc                                                                                                                                     <  1:1224588/140‑1 (MQ=255)
   aCCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCtagctcgc                                                                                                                                   >  2:1888381/1‑140 (MQ=255)
   aCCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCtagctcgc                                                                                                                                   >  1:2237485/1‑140 (MQ=255)
     cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTGGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGCTGCGTGTTGAAGCTAGCTCGCCg                                                                                                                                 >  2:291518/1‑140 (MQ=255)
     cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCg                                                                                                                                 <  1:2878615/140‑1 (MQ=255)
     cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCg                                                                                                                                 >  2:2353900/1‑140 (MQ=255)
     cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCg                                                                                                                                 >  1:1875825/1‑140 (MQ=255)
      gTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCgg                                                                                                                                <  2:376022/140‑1 (MQ=255)
       tcacATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGa                                                                                                                               <  2:2962110/137‑1 (MQ=255)
       tCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGa                                                                                                                               <  2:1010560/140‑1 (MQ=255)
       tCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGa                                                                                                                               <  2:2291439/140‑1 (MQ=255)
          cATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATgg                                                                                                                            >  1:859437/1‑140 (MQ=255)
          cATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATgg                                                                                                                            >  2:1378142/1‑140 (MQ=255)
                  tCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTg                                                                                                                    >  2:561544/1‑140 (MQ=255)
                   cAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGa                                                                                                                   <  1:1321630/140‑1 (MQ=255)
                    aGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAg                                                                                                                  <  2:642098/140‑1 (MQ=255)
                        tGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCgg                                                                                                              <  1:2853506/140‑1 (MQ=255)
                          aaGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGtt                                                                                                            >  1:1231675/1‑140 (MQ=255)
                           agagGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGttt                                                                                                           >  2:2989230/1‑139 (MQ=255)
                           agaaGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTa                                                                                                           >  2:1674514/1‑140 (MQ=255)
                            gagGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAg                                                                                                          <  1:423516/140‑1 (MQ=255)
                              ggggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAg                                                                                                                       <  2:1388705/125‑1 (MQ=255)
                              ggggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAg                                                                                                                       >  1:1388705/1‑125 (MQ=255)
                               gggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCgg                                                                                                       >  2:3060558/1‑139 (MQ=255)
                                gtggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtc                                                                                                      >  1:1352317/3‑140 (MQ=255)
                                ggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtc                                                                                                      <  2:938682/140‑1 (MQ=255)
                                ggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtc                                                                                                      >  1:1194169/1‑140 (MQ=255)
                                  ggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtcat                                                                                                    <  1:1888381/140‑1 (MQ=255)
                                    aTTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtcatca                                                                                                  >  1:2605083/1‑140 (MQ=255)
                                    aTTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGtcatc                                                                                                   <  1:1670375/139‑1 (MQ=255)
                                         tGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGTCCTCAaacat                                                                                             >  2:766629/1‑140 (MQ=255)
                                           tCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATcc                                                                                           >  2:153741/1‑140 (MQ=255)
                                           tCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATcc                                                                                           >  2:1591117/1‑140 (MQ=255)
                                             tGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCtt                                                                                         >  2:1758201/1‑140 (MQ=255)
                                             tGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCtt                                                                                         >  2:57398/1‑140 (MQ=255)
                                             tGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCtt                                                                                         >  2:1954595/1‑140 (MQ=255)
                                             tGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCtt                                                                                         >  2:2819342/1‑140 (MQ=255)
                                              gCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCttt                                                                                        >  2:2860931/1‑140 (MQ=255)
                                                gttgttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTg                                                                                      >  2:2081769/1‑140 (MQ=255)
                                                 ttgttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgc                                                                                     <  2:2605083/140‑1 (MQ=255)
                                                  tgttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGGCCTCAATCATCCTTTTtcg                                                                                    >  2:747409/1‑140 (MQ=255)
                                                  tgttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcg                                                                                    >  1:261343/1‑140 (MQ=255)
                                                   gttACCTGAGTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   <  2:2116655/140‑1 (MQ=255)
                                                   gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   <  2:1854984/140‑1 (MQ=255)
                                                   gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   <  2:261343/140‑1 (MQ=255)
                                                   gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   <  2:1181208/140‑1 (MQ=255)
                                                   gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   <  2:547593/140‑1 (MQ=255)
                                                     tACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTg                                                                                 >  2:343174/1‑140 (MQ=255)
                                                           aaTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTa                                                                           >  1:2556771/1‑140 (MQ=255)
                                                           aaTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTa                                                                           >  2:1933393/1‑140 (MQ=255)
                                                            aTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTAc                                                                          >  1:573645/1‑140 (MQ=255)
                                                             tGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTgcgc                                                                                   >  1:2932710/1‑130 (MQ=255)
                                                             tGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACa                                                                         <  2:159511/140‑1 (MQ=255)
                                                               tAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGt                                                                       <  1:1758201/140‑1 (MQ=255)
                                                                 gCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTgg                                                                     >  2:3049751/1‑140 (MQ=255)
                                                                  cAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGc                                                                    >  1:972720/1‑140 (MQ=255)
                                                                  cAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGc                                                                    >  1:1628514/1‑140 (MQ=255)
                                                                     cTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCAcc                                                                 <  1:2819342/140‑1 (MQ=255)
                                                                     cTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCAcc                                                                 >  1:1118244/1‑140 (MQ=255)
                                                                     cTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCAcc                                                                 <  1:1502293/140‑1 (MQ=255)
                                                                     cTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCAcc                                                                 <  1:29358/140‑1 (MQ=255)
                                                                     cTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCAcc                                                                 <  2:1314931/140‑1 (MQ=255)
                                                                          gggTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGaa                                                            <  1:561544/140‑1 (MQ=255)
                                                                          gggTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGaa                                                            <  1:739250/140‑1 (MQ=255)
                                                                          gggTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGaa                                                            <  1:679281/140‑1 (MQ=255)
                                                                              aaaCTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGa                                                        <  1:2563119/140‑1 (MQ=255)
                                                                               aaCTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGaa                                                       <  1:1151196/140‑1 (MQ=255)
                                                                               aaCTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGaa                                                       >  1:1965654/1‑140 (MQ=255)
                                                                                     ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGaa                                                       <  2:1931258/134‑1 (MQ=255)
                                                                                     ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGaa                                                       >  1:1931258/1‑134 (MQ=255)
                                                                                     ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAg                                                 >  2:2170767/1‑140 (MQ=255)
                                                                                      tACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGc                                                >  2:746072/1‑140 (MQ=255)
                                                                                           ggACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTa                                           >  2:2375677/1‑140 (MQ=255)
                                                                                            gACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGGCATCAATAACCCTTTTTTGCTAGGCGCACAGTGGGGCCCCGGAAGGGACAGTAACgcggac                                          >  2:1545309/1‑135 (MQ=255)
                                                                                              cAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTAcgc                                        >  2:1887127/1‑140 (MQ=255)
                                                                                              cAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTAcgc                                        >  1:1067446/1‑140 (MQ=255)
                                                                                                aGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACgcgc                                      >  1:1831623/1‑140 (MQ=255)
                                                                                                aGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGGACgcgc                                      >  1:2081382/1‑140 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTATTCCACGGTTGCGTGTTGAAGCTAGCTCGACGGATGGCTTAGGTGAGCCGGTTAGTGTCATTAACCACCGTTTTGCTCTGGGCGTACAGTTCCGCCCGGAATGCAACAGTAGCGGGTACGCGCt                                     <  2:2911305/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCTTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  1:2218798/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  1:762915/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  1:353954/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  1:397360/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  1:1010267/140‑1 (MQ=255)
                                                                                                 ggggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCt                                     <  2:1067446/140‑1 (MQ=255)
                                                                                                   ggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCttt                                   <  1:2866643/140‑1 (MQ=255)
                                                                                                     cGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTcg                                 <  2:1965654/140‑1 (MQ=255)
                                                                                                      gAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTcgc                                <  2:1252689/140‑1 (MQ=255)
                                                                                                       caGGTCGTTAGCCCACCGTTGCGTGTTGAAGCTAGCCCGCCGGATGGCTTAGTTGAGGCGGTTTATCTGATCAATCATTCTTTTTTGCTCGGAGTACACCGGCACCCGGGCGGGGACCGAATCGTGTAAGCACTTTcgcc                               >  2:218604/2‑139 (MQ=255)
                                                                                                       aaGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTcgcg                               >  1:308519/1‑140 (MQ=255)
                                                                                                       aaGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGGACAGGGGGACCCGGGATGGAAAAGTAGCGAGGACGCGCTTTcgcg                               >  2:1174227/1‑140 (MQ=255)
                                                                                                       aaGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCCTCATCAATCATCCTTTTTCGCCGGGCGTACAGGGGGACCCCGAATGGAACAGTAGCGAGGAAGCGCTTTcccg                               >  2:386553/1‑140 (MQ=255)
                                                                                                             gTTAGCCCACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTGATCCATCATCCTTTTTCGCCTGGCGTGAAGGGGGACCCGGAATAGAACAGAAGAGAGAACGCGCTTTCGCGGATAtt                         >  2:1977241/1‑140 (MQ=255)
                                                                                                                  cccACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCg                    <  2:8087/140‑1 (MQ=255)
                                                                                                                    cACGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAg                  <  1:1028823/140‑1 (MQ=255)
                                                                                                                      cGGTTGCGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGGATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAggg                >  1:2524561/1‑140 (MQ=255)
                                                                                                                            cGTGTTGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTAt          >  2:2041985/1‑140 (MQ=255)
                                                                                                                                 tGAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCg     <  1:3049751/140‑1 (MQ=255)
                                                                                                                                  gAAGCTAGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGc    <  1:2041985/140‑1 (MQ=255)
                                                                                                                                    agctagctCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCtt  <  1:1545309/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTT  >  NC_000913/1361578‑1361849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: