Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,981,526 (C)5→6 coding (87/1425 nt) otsA ← trehalose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,981,5211.C99.6% 948.2 / ‑3.7 233coding (92/1425 nt)otsAtrehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (117/115);  total (118/115)
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.97e-01

Only 100 of 262 total aligned reads displayed.

GGCCCACGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACGACTCATAGTC  >  NC_000913/1981421‑1981616
                                                                                                     |                                                                                               
gacccACGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTg                                                                                                  >  2:2581154/3‑101 (MQ=255)
 acccACGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGc                                                                                                >  1:2581157/2‑101 (MQ=255)
 acccACGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGc                                                                                                >  1:2581155/2‑101 (MQ=255)
     aCGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccc                                                                                            >  2:2581161/1‑101 (MQ=255)
     aCGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccc                                                                                            <  1:2580858/101‑1 (MQ=255)
       gTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGcccccca                                                                                          >  2:2581165/1‑99 (MQ=255)
        tAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccc                                                                                            >  1:2581168/1‑98 (MQ=255)
         aaTGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccccagt                                                                                        >  1:2581170/1‑97 (MQ=255)
          aTGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccccagta                                                                                       >  2:2581172/1‑96 (MQ=255)
          aTGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGccccccagta                                                                                       <  1:2580825/101‑6 (MQ=255)
           tGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTAt                                                                                       <  1:2580970/101‑1 (MQ=255)
           tGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTAt                                                                                       >  2:2581177/1‑101 (MQ=255)
            gTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATg                                                                                      <  1:2580629/101‑1 (MQ=255)
            gTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATg                                                                                      >  2:2581179/1‑101 (MQ=255)
             ttACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGc                                                                                     >  1:2581180/1‑101 (MQ=255)
               aCCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCa                                                                                   >  2:2581187/1‑101 (MQ=255)
                ccTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATg                                                                                      <  2:2580865/97‑1 (MQ=255)
                ccTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCaa                                                                                  <  1:2580897/101‑1 (MQ=255)
                ccTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCaa                                                                                  >  1:2581188/1‑101 (MQ=255)
                  ttttttCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACg                                                                                >  2:2581192/1‑101 (MQ=255)
                   tttttCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACgg                                                                               <  2:2580781/101‑1 (MQ=255)
                     tttCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCa                                                                             >  2:2581195/1‑101 (MQ=255)
                        cACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAgg                                                                          >  1:2581196/1‑101 (MQ=255)
                          ccTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAg                                                                           <  2:2580966/98‑1 (MQ=255)
                          ccTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGcc                                                                        >  1:2581202/1‑101 (MQ=255)
                          ccTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGcc                                                                        >  2:2581203/1‑101 (MQ=255)
                            ttttttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCAc                                                                      <  1:2580787/101‑1 (MQ=255)
                             tttttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCAcc                                                                     >  2:2581208/1‑101 (MQ=255)
                              ttttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCg                                                                    >  1:2581211/1‑101 (MQ=255)
                               tttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCgg                                                                   >  2:2581213/1‑101 (MQ=255)
                                ttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGc                                                                  <  1:2580701/101‑1 (MQ=255)
                                ttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGc                                                                  >  1:2581218/1‑101 (MQ=255)
                                ttAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGc                                                                  >  2:2581217/1‑101 (MQ=255)
                                 tAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCAc                                                                      >  2:2581220/1‑96 (MQ=255)
                                     ggCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTgg                                                             >  2:2581224/1‑101 (MQ=255)
                                        tGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTggcgg                                                          >  1:2581227/1‑101 (MQ=255)
                                        tGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTggcgg                                                          >  1:2581229/1‑101 (MQ=255)
                                           tCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGt                                                       >  1:2581232/1‑101 (MQ=255)
                                             cTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGc                                                     <  1:2580692/101‑1 (MQ=255)
                                              tCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCt                                                    >  2:2581237/1‑101 (MQ=255)
                                              tCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCt                                                    >  1:2581236/1‑101 (MQ=255)
                                               cATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTc                                                   <  1:2581016/101‑1 (MQ=255)
                                                aTTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCg                                                  >  1:2581240/1‑101 (MQ=255)
                                                aTTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCg                                                  >  2:2581241/1‑101 (MQ=255)
                                                aTTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCg                                                  <  1:2581060/101‑1 (MQ=255)
                                                   cccTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTggcg                                                           >  2:2581247/1‑89 (MQ=255)
                                                   cccTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCt                                               <  1:2580867/101‑1 (MQ=255)
                                                   cccTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCt                                               >  1:2581246/1‑101 (MQ=255)
                                                      tGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCtggt                                            <  2:2580793/101‑1 (MQ=255)
                                                       gTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCtgg                                             >  2:2581251/1‑99 (MQ=255)
                                                        tttCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCtggtgg                                          >  2:2581254/1‑101 (MQ=255)
                                                           caccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGc                                       <  2:2580908/101‑1 (MQ=255)
                                                           caccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGc                                       <  1:2580800/101‑1 (MQ=255)
                                                           caccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGcccccca                                                                                          <  1:2580823/49‑3 (MQ=255)
                                                             ccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGc                                                                              <  2:2580751/60‑1 (MQ=255)
                                                             ccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTggcgg                                                          <  2:2580625/80‑1 (MQ=255)
                                                             ccacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGc                                       >  1:2581257/1‑99 (MQ=255)
                                                              cacTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAAt                                    >  2:2581258/1‑101 (MQ=255)
                                                                cTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAg                                                                           <  1:2580481/60‑1 (MQ=255)
                                                                cTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCa                                                                 <  2:2580884/70‑1 (MQ=255)
                                                                cTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATcc                                  <  1:2580953/101‑1 (MQ=255)
                                                                cTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGcccccca                                                                                          <  2:2580835/44‑3 (MQ=255)
                                                                 tCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCg                                 <  1:2580921/101‑1 (MQ=255)
                                                                  ccAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCAc                                                                <  2:2580702/69‑1 (MQ=255)
                                                                  ccAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCtgg                                             <  2:2580947/88‑1 (MQ=255)
                                                                  ccAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCg                                 <  2:2580923/100‑1 (MQ=255)
                                                                    aGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGtt                              >  1:2581261/1‑101 (MQ=255)
                                                                     gCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCgg                                <  2:2581020/98‑1 (MQ=255)
                                                                      ccAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTggc                                                            <  2:2581040/69‑1 (MQ=255)
                                                                      ccAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAg                            >  2:2581264/1‑101 (MQ=255)
                                                                       cAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCgg                                <  1:2580818/96‑1 (MQ=255)
                                                                       cAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGa                           >  1:2581267/1‑101 (MQ=255)
                                                                       cAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGa                           <  2:2580930/101‑1 (MQ=255)
                                                                       cAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGa                           >  2:2581268/1‑101 (MQ=255)
                                                                       cAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGa                           <  1:2580813/101‑1 (MQ=255)
                                                                        aaaCCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGAt                          <  2:2580980/101‑1 (MQ=255)
                                                                        aaaCCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGAt                          >  1:2581270/1‑101 (MQ=255)
                                                                         aaCCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACg                                                                                <  2:2580673/46‑1 (MQ=255)
                                                                         aaCCACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGAta                         <  1:2580935/101‑1 (MQ=255)
                                                                           ccACAGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGAtact                       <  1:2580730/101‑1 (MQ=255)
                                                                              caGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGAtactact                    <  1:2580909/101‑2 (MQ=255)
                                                                               aGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGa                   <  1:2580861/101‑1 (MQ=255)
                                                                               aGTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGa                   >  1:2581275/1‑101 (MQ=255)
                                                                                gTCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACg                                                                                <  1:2580843/39‑1 (MQ=255)
                                                                                 tCCGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACt                 >  2:2581278/1‑101 (MQ=255)
                                                                                  ccGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTgg                                                             <  2:2580740/56‑1 (MQ=255)
                                                                                  ccGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTc                                                   <  1:2580994/66‑1 (MQ=255)
                                                                                  ccGCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTa                <  2:2580955/101‑1 (MQ=255)
                                                                                    gCCTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTaaa              >  1:2581280/1‑101 (MQ=255)
                                                                                     ccTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTaaa              >  1:2581282/1‑100 (MQ=255)
                                                                                     ccTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAAc             <  2:2580937/101‑1 (MQ=255)
                                                                                     ccTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAAc             <  2:2580775/101‑1 (MQ=255)
                                                                                     ccTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAAc             <  1:2580886/101‑1 (MQ=255)
                                                                                     ccTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAAc             <  1:2580795/101‑1 (MQ=255)
                                                                                      cTGCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACg            >  2:2581283/1‑101 (MQ=255)
                                                                                        gCGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACt                 <  2:2580839/94‑1 (MQ=255)
                                                                                         cGGCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTgg                                                             >  1:2581287/1‑49 (MQ=39)
                                                                                          ggCTTTCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACGACTc        >  2:2581290/1‑101 (MQ=255)
                                                                                             tttCAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACGACTCATa     >  1:2581292/1‑101 (MQ=255)
                                                                                                cAGTGCCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACGACTCATAGtc  >  1:2581296/1‑101 (MQ=255)
                                                                                                     |                                                                                               
GGCCCACGTAATGTTACCTTTTTTCACCTTTTTTAGCGGCTGATCCTCATTCCCTGTTTCACCACTCCAGCCAAACCACAGTCCGCCTGCGGCTTTCAGTGCCCCCAGTATGCCAACGGCAAGGCCACCGGCACTGGCGGCGTGCTCGTCTGGTGGTGCAATCCGGTTAGATACTACGACTAAACGACTCATAGTC  >  NC_000913/1981421‑1981616

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