Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,851,206 G→T 6.1% G116V (GGG→GTG) ‡ ansA → L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,851,2060GT 6.1% 239.0 / 3.4 72G116V (GGG→GTG) ‡ansAL‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base T (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCC  >  NC_000913/1851060‑1851351
                                                                                                                                                  |                                                                                                                                                 
cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                >  2:1242198/1‑150 (MQ=255)
cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                <  1:242744/150‑1 (MQ=255)
  gccagcaTACTCCCGAAGATAATAAAGCCCACTATGACGACTATGAGGGGTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcac                                                                                                                                              <  1:1902556/148‑1 (MQ=255)
   gcagcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                                                                                                                                             >  2:616171/1‑150 (MQ=255)
   gcagcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                                                                                                                                             >  2:328477/1‑150 (MQ=255)
     agcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaat                                                                                                                                           <  1:1868100/150‑2 (MQ=255)
      gcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAAt                                                                                                                                          >  1:1700037/1‑150 (MQ=255)
         tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                       <  2:1441593/150‑6 (MQ=255)
         tatTGCTGAAGATATTAAAGCGAACTATGACGACTATGATGGTTTTGTCATTCTGCACGCCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTAGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATccc                                                                                                                                       >  1:1650301/1‑150 (MQ=255)
          atTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCg                                                                                                                                      <  2:177448/150‑1 (MQ=255)
           ttGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgc                                                                                                                                     >  2:1090510/1‑150 (MQ=255)
             gCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                   >  1:1118157/1‑150 (MQ=255)
             gCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                   <  2:454701/150‑1 (MQ=255)
              cTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                   >  1:1339973/1‑149 (MQ=255)
              cTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                   <  2:1339973/149‑1 (MQ=255)
               tGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                                 >  2:1207259/1‑150 (MQ=255)
               tGAAGATATTAAAGCGCACTATCAAGACTATGATGGTATTGTAATTATGCACCGCACAGACACGATGGAGTATAACGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                                 >  2:927720/1‑150 (MQ=255)
                        tAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACg                                                                                                                        >  2:131171/1‑150 (MQ=255)
                        tAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACg                                                                                                                        >  2:138358/1‑150 (MQ=255)
                            gcgcACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGctct                                                                                                                    <  1:877388/150‑1 (MQ=255)
                             cgcACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTg                                                                                                                   <  2:85350/150‑1 (MQ=255)
                               cACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGAc                                                                                                                 >  2:137958/1‑150 (MQ=255)
                                  tATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGa                                                                                                              <  2:625655/150‑1 (MQ=255)
                                   aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGAc                                                                                                             <  1:208065/150‑1 (MQ=255)
                                    tGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACa                                                                                                            >  1:1090534/1‑150 (MQ=255)
                                    tGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACa                                                                                                            >  1:116719/1‑150 (MQ=255)
                                       cgacTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAAt                                                                                                         <  1:795778/150‑1 (MQ=255)
                                         acTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTa                                                                                                       <  1:616171/150‑1 (MQ=255)
                                               atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgc                                                                                                 <  2:516934/150‑1 (MQ=255)
                                                         cATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGtt                                                                                       <  1:328477/150‑1 (MQ=255)
                                                           ttCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGt                                                                                     >  1:457319/1‑150 (MQ=255)
                                                               gCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGt                                                                                 >  1:495499/1‑150 (MQ=255)
                                                               gCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGt                                                                                 >  1:879790/1‑150 (MQ=255)
                                                                 ccGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCTTGCTCGAGAATCTCGGTACACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAAATAATCTGCTGACTGCGTTGTACGTTg                                                                               <  2:1847269/149‑1 (MQ=255)
                                                                 aCGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTg                                                                               >  1:1091179/1‑150 (MQ=255)
                                                                    gCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTTTTTTTGACAGTGTCAGAAATCCGGCTGGTTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGAGTTGTAGGTTgagg                                                                            <  2:189199/150‑1 (MQ=255)
                                                                    gCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcgg                                                                            >  1:1599931/1‑150 (MQ=255)
                                                                     cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                           <  1:1207259/150‑1 (MQ=255)
                                                                         gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                       <  2:952635/150‑1 (MQ=255)
                                                                         gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                       <  2:1118157/150‑1 (MQ=255)
                                                                         gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                       <  1:1090510/150‑1 (MQ=255)
                                                                             cGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATc                                                                   >  2:1117538/1‑150 (MQ=255)
                                                                             cGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGAGTTGTACGTTGCGGCGAATTATc                                                                   >  2:1621612/1‑150 (MQ=255)
                                                                                  gCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                              <  1:131171/150‑1 (MQ=255)
                                                                                  gCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                              <  1:138358/150‑1 (MQ=255)
                                                                                         ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                           <  2:1733235/130‑1 (MQ=255)
                                                                                         ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                           >  1:1733235/1‑130 (MQ=255)
                                                                                         ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                       >  2:1102775/1‑150 (MQ=255)
                                                                                         ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                       <  1:137958/150‑1 (MQ=255)
                                                                                           gCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAACCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTa                                                     >  2:1127062/1‑150 (MQ=255)
                                                                                               ctGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGc                                                 <  2:1700037/150‑1 (MQ=255)
                                                                                                tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGct                                                >  1:924999/1‑150 (MQ=255)
                                                                                                       tCGTTCATGCTCGAGAATCTCGGTAACCCGGTGATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGAGAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCtttttg                                         <  1:152405/150‑2 (MQ=255)
                                                                                                       tCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                         <  2:1090534/150‑1 (MQ=255)
                                                                                                       tCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                         <  2:631565/150‑1 (MQ=255)
                                                                                                       gcGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                         <  2:1259029/149‑1 (MQ=255)
                                                                                                        cGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCttttttt                                        <  1:927720/150‑3 (MQ=255)
                                                                                                         gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                       <  1:1600986/150‑1 (MQ=255)
                                                                                                           tCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGTTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                     <  2:116719/150‑1 (MQ=255)
                                                                                                               gCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCg                                 >  1:1864241/1‑150 (MQ=255)
                                                                                                                         ctcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgc                       >  2:1517494/1‑150 (MQ=255)
                                                                                                                              tAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCaa                  >  1:1526862/1‑150 (MQ=255)
                                                                                                                                aaCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAAcc                >  1:139040/1‑150 (MQ=255)
                                                                                                                                 aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg               >  2:1972246/1‑150 (MQ=255)
                                                                                                                                 aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg               >  1:1359387/1‑150 (MQ=255)
                                                                                                                                    ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCgg                     <  2:1506461/141‑1 (MQ=255)
                                                                                                                                    ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCgg                     >  1:1506461/1‑141 (MQ=255)
                                                                                                                                     ggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCCGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt           >  2:408036/5‑150 (MQ=255)
                                                                                                                                       cATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAc         >  1:1000995/1‑150 (MQ=255)
                                                                                                                                       atttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAc         >  2:1307359/3‑150 (MQ=255)
                                                                                                                                        aTTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAcc        >  2:741487/1‑150 (MQ=255)
                                                                                                                                              aCAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc  >  1:771413/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCC  >  NC_000913/1851060‑1851351

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