Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,345 T→A 15.7% T232S (ACG→TCG) ‡ ddpF ← D,D‑dipeptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,3450TA15.7% 314.0 / 41.1 127T232S (ACG→TCG) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (46/60);  new base A (9/11);  total (55/72)
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.68e-01

Only 100 of 129 total aligned reads displayed.

GCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAA  >  NC_000913/1557246‑1557442
                                                                                                   |                                                                                                 
gCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGt                                                                                                   <  2:1420892/100‑1 (MQ=255)
  aGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCa                                                                                                 <  2:950255/100‑1 (MQ=255)
  aGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCa                                                                                                 <  2:1041947/100‑1 (MQ=255)
    aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGc                                                                                               <  1:1093385/100‑1 (MQ=255)
    aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGc                                                                                               <  1:2791395/100‑1 (MQ=255)
      ccGTTTTACGTAATGCCCATTCTTCCTCCAGCCGTTTATCAATGGCGGGGAGGGAATCCCGCAATAATCGGGTGTATGGCTGTGCAGGTGGCGTCAGc                                                                                               >  2:3272913/1‑98 (MQ=255)
       cGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCAcc                                                                                            <  1:1447908/100‑1 (MQ=255)
       cGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAACGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCAcc                                                                                            <  1:3605225/100‑1 (MQ=255)
        gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACct                                                                                           >  2:2183492/1‑100 (MQ=255)
        gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACct                                                                                           <  2:2998912/100‑1 (MQ=255)
         ttttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                          <  2:2371122/100‑1 (MQ=255)
         ttttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                          <  1:3777192/100‑1 (MQ=255)
          tttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGACGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgc                                                                                         <  2:1174797/100‑1 (MQ=255)
           ttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgca                                                                                        <  1:3162402/100‑2 (MQ=255)
            tACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcac                                                                                       <  2:444405/100‑3 (MQ=255)
            tACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgctg                                                                                       >  1:3515325/1‑100 (MQ=255)
             aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcaca                                                                                      <  2:2561686/100‑4 (MQ=255)
             aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcaca                                                                                      <  1:2182858/100‑4 (MQ=255)
             aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcaca                                                                                      <  1:2159109/100‑4 (MQ=255)
             aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcaca                                                                                      <  1:3748445/100‑4 (MQ=255)
             aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCtgctgc                                                                                      <  2:2609808/100‑1 (MQ=255)
              cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcg                                                                                     <  2:2772802/100‑1 (MQ=255)
              cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcg                                                                                     <  1:2647272/100‑1 (MQ=255)
              cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcg                                                                                     <  2:3657175/100‑1 (MQ=255)
              cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcg                                                                                     >  1:1047285/1‑100 (MQ=255)
                tAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGcgcg                                                                                   <  1:193937/100‑1 (MQ=255)
                 aaTGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcacatcca                                                                                  <  2:1457929/100‑8 (MQ=38)
                   tGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCg                                                                                >  1:2360719/1‑100 (MQ=255)
                    gCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGGTGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGc                                                                               <  1:2726515/100‑1 (MQ=255)
                    gCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGc                                                                               <  2:23123/100‑1 (MQ=255)
                      ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGGAGGTGCCGTCAGCACCTGCTGCGCGTCGc                                                                               >  2:864072/1‑98 (MQ=255)
                      ccATTCTTCCTCCAGCGGTGTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGGGTATGGAGGTGCAGGGGGCGACAGCACCTGCTGCGTGTCGc                                                                               >  2:2320478/1‑98 (MQ=255)
                      ccATTCTTCCTCCAGCGGTGTATCAATGGCGCGGAGGGCATCCCGCAATAATTGGGTGTCTGGGTGTGCACGTGGCGTCAGCACCTGCTGGGCGGCGcc                                                                              >  2:373418/1‑99 (MQ=255)
                       cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCa                                                                            >  1:479796/1‑100 (MQ=255)
                         ttcttcCGCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGt                                                                          <  2:1440135/100‑1 (MQ=255)
                           cttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTc                                                                        <  1:234124/100‑1 (MQ=255)
                           cttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTc                                                                        <  1:204634/100‑1 (MQ=255)
                              cctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                     >  1:2406468/1‑100 (MQ=255)
                                tccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACaa                                                                   <  2:3214993/100‑1 (MQ=255)
                                  cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATc                                                                 <  1:340942/100‑1 (MQ=255)
                                  cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATc                                                                 <  2:639125/100‑1 (MQ=255)
                                    gCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAGGGATGTGCAGGTGCCGTCAGCACCTTCTGCGCGTCGCTCAGTTCTACAATCTg                                                               >  1:858680/1‑100 (MQ=255)
                                      ggTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                             <  1:3483249/100‑1 (MQ=255)
                                      ggTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                             >  1:3460469/1‑100 (MQ=255)
                                      ggTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                             <  1:2732077/100‑1 (MQ=255)
                                      ggTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                             >  2:837012/1‑100 (MQ=255)
                                         ttATCAAGGGCGGGGAGGGAATCCAGCAATAATCGTGTGTATGGATGTGCAGGTGCTGTCAGCACCTGCTGGGCGTCGCCCAGTTCTACAATCTGCCCg                                                           >  2:2435004/1‑99 (MQ=255)
                                          tATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgag                                                         <  2:2602286/100‑1 (MQ=255)
                                            tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                       <  1:473860/100‑1 (MQ=255)
                                            tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                       >  1:2768682/1‑100 (MQ=255)
                                              aaTGGCGGGGCGGGAATCCAGCAATAATTGGGTGTTAGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGCCGCCCAGTTCTACAAACTGCCCGAGATa                                                      >  2:2054423/1‑99 (MQ=255)
                                               aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACa                                                    <  1:903748/100‑1 (MQ=255)
                                               aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACa                                                    <  2:1036529/100‑1 (MQ=255)
                                               aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACa                                                    <  1:3498431/100‑1 (MQ=255)
                                               aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACa                                                    >  1:1902262/1‑100 (MQ=255)
                                               aTGGCGGGGAGGGAATCCAGCAATAAGCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACa                                                    <  1:779200/100‑1 (MQ=255)
                                                  gCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCCGTTCTACAATCTGCCCGAGATACATCa                                                 <  1:713054/100‑1 (MQ=255)
                                                     gggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                              >  1:3145894/1‑100 (MQ=255)
                                                         gggaATCCAGCAATAATCGGGTGTATGGATGTTCAGGTGTCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                          <  2:2684588/100‑1 (MQ=255)
                                                         gggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                          <  2:3022774/100‑1 (MQ=255)
                                                          ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAcc                                         <  1:3685211/100‑1 (MQ=255)
                                                          ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAcc                                         >  1:2804877/1‑100 (MQ=255)
                                                           gaATCCAGCACTAATCGGGCGCCTGGATTTGCAGGTGCCGTCAGCACCTGCTGGGCGTCGCCCAGTTCTACAATCTGCCCCAGCTACATCACCGCCAccc                                        >  2:1295137/1‑100 (MQ=255)
                                                              tCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGAt                                     <  1:76669/100‑1 (MQ=255)
                                                               ccAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                    <  1:3558554/100‑1 (MQ=255)
                                                                cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                   >  1:1500113/1‑100 (MQ=255)
                                                                cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                   <  1:78838/100‑1 (MQ=255)
                                                                cAGCAATAATCGGGTGTATGGATGTGCAGGCGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                   <  2:1489799/100‑1 (MQ=255)
                                                                cAGCAATAATCAGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGa                                      >  2:785381/1‑97 (MQ=255)
                                                                   cAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTc                                <  2:2804331/100‑1 (MQ=255)
                                                                     ataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCat                              >  2:1760396/1‑100 (MQ=255)
                                                                       aatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                            <  2:2793836/100‑1 (MQ=255)
                                                                        atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCCCTTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATg                           <  2:3813638/100‑1 (MQ=255)
                                                                           gggTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACg                        <  1:1971800/100‑1 (MQ=255)
                                                                            ggTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGt                       >  2:3699414/1‑100 (MQ=255)
                                                                            ggTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGt                       <  1:1669281/100‑1 (MQ=255)
                                                                             gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTa                      >  1:228330/1‑100 (MQ=255)
                                                                             gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTa                      >  1:472317/1‑100 (MQ=255)
                                                                               gtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTGCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGt                       >  2:3308148/1‑97 (MQ=255)
                                                                                   ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                <  1:2571740/100‑1 (MQ=255)
                                                                                   cgcgcagCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                >  1:2381220/7‑100 (MQ=38)
                                                                                     cgcagCAGGTGCTGACGGCACCTGTTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCACATGACGGATCAc                  >  2:1281064/5‑96 (MQ=37)
                                                                                     cgcagCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgag              >  1:586996/5‑100 (MQ=38)
                                                                                     cgcagCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATAGCTCATACGACGTATGACCga               >  2:3666026/5‑99 (MQ=37)
                                                                                     aTGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgag              >  1:1156127/1‑100 (MQ=255)
                                                                                      gcagCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgaga             >  1:909330/4‑100 (MQ=38)
                                                                                       gtgCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGAc            >  1:1497299/1‑100 (MQ=255)
                                                                                        tgCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACa           >  1:1113160/1‑100 (MQ=255)
                                                                                         gCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTGTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt          >  2:1301658/1‑100 (MQ=39)
                                                                                         gCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt          <  2:2073681/100‑1 (MQ=255)
                                                                                         gCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt          <  2:3793169/100‑1 (MQ=255)
                                                                                          cAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCTAGACAtt         >  2:556677/1‑100 (MQ=255)
                                                                                           aGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCAGCGAGACATtg        >  1:3138349/1‑100 (MQ=38)
                                                                                           aGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg        <  1:2885892/100‑1 (MQ=255)
                                                                                            ggTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt       >  1:1753467/1‑100 (MQ=255)
                                                                                               gCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgtg    <  2:3069930/100‑1 (MQ=255)
                                                                                                ccGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa   <  2:3421388/100‑1 (MQ=255)
                                                                                                ccGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa   >  1:2016805/1‑100 (MQ=255)
                                                                                                 cGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaa  <  2:157383/100‑1 (MQ=255)
                                                                                                 cGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa   >  2:1602280/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
GCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAA  >  NC_000913/1557246‑1557442

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